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Summary 



This thesis concerns the compatibility of inhomogeneous cosmologies with our present understanding of the universe. It is a 
problem of some interest to find the class of all relativistic cosmological models which are capable of providing a reasonable 
'fit' to the universe. One can imagine building up an (infinite-dimensional) parameter space containing all cosmological 
models. At any time the understanding of the universe would represent a blob in parameter space in which, presumably, 
the real universe would sit. This thesis, in some respects, is part of this process. We consider Stephani models, which are 
a generalisation of the standard Friedmann-LemaTtre-Robertson-Walker (FLRW) models, which can be thought of as FLRW 
models with curvature which changes over time. This changing curvature reflects the existence of spatial pressure gradients 
which leads to an acceleration of the fundamental observers. Thus these models generalise the 'dust' assumption of standard 
cosmology. 

Models normally considered in the 'classification scheme' approach are usually homogeneous dust or barotropic perfect 
fluid models. They are anisotropic however, and thus generalise the FLRW models. Most importantly, because these models 
are homogeneous, they satisfy the Copernican principle. The crucial aspect of this work is the retention of the Copernican 
principle - an assumption regarded by many as crucial to cosmology. It states that we are not at a special location in the 
universe. This is a vital aspect of the original work in this thesis: consideration of an inhomogeneous model, while retaining 
the Copernican principle has, as far as the author is aware, not been considered in detail before. 

One may formulate the Copernican principle in many ways, from assuming we are not at a special location, to assuming 
that all (or most) locations are equivalent, which more or less forces homogeneity. Because the models considered here 
are inhomogeneous, they cannot satisfy the stronger version of the Copernican principle entirely - all locations will not be 
equivalent. However, we may demand that they are observationally indistinguishable. This is the tactic we use here. En route 
to this goal we must therefore calculate all observable quantities at any location in the spacetime. Certain properties of the 
Stephani models we consider allow us to do this exactly; consequently, many results of this thesis present, for the first time, 
observational relations for a class of inhomogeneous cosmological models which are exact, and valid for any observer position 
in the spacetime. 

It may reasonably be claimed that the standard model is perfectly acceptable. However, a number of the properties of the 
models considered here do make them rather appealing. For example it is shown in §|]^ these models do not suffer from the 
horizon problem which is prevalent in the standard model. Also, the current conflict between cosmologists measuring a small 
but non-zero cosmological constant, and particle physics requiring it to be either zero or one hundred and twenty orders of 
magnitude higher, may be motivation in itself for considering cosmologies with a somewhat non-standard matter content. 

In chapter |l| a brief review of the Copernican and cosmological principles in in homogeneous cosmologies is presented. A 
discussion is given of the present understanding of the important parameters of the standard model, some methods used to 
find these parameters, and some of the problems encountered. 

In chapter|| an overview of relativistic cosmology is presented. This is necessarily incomplete, with an emphasis on deriving 
observable relations, and on the FLRW models. The 1+3 formalism is presented. Some new results are presented concerning 
conformally related spacetimes (which turn out to be useful but not necessary for deriving the observational quantities in 
the Stephani models). FLRW models are then reviewed. We give a discussion of the generalised definitions of the Hubble 
constant and deceleration parameter in non-standard cosmologies. 

In chapter || we discuss a theorem of fundamental importance to cosmology - the Ehlers-Geren-Sachs theorem (1968). 
This condition singles out spacetimes which will allow an isotropic cosmic microwave background (CMB). When applied to 
a dust universe it says that the existence of an isotropic CMB for every observer in the spacetime implies that the universe 
is FLRW. We extend this theorem to include the case of non-geodesic observers (in a perfect fluid model), which singles out 
a subclass of the Stephani models with symmetry. These models form the basis of the rest of the thesis. This chapter is a 
slightly extended version of Clarkson and Barrett (1999). 

Chapter takes the models which were singled out in chapter |3| and derives the observational relations for these models. 
Consideration is given to non-central observers, and the observational relations are then derived for any location in the 
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spacetime. 

The following two chapters examine the models in some detail from every observer position. The 'worst case' - ie, the 
observer position most restrictive on the parameter space - is singled out for each constraint, and it is shown that there is 
a large area of parameter space which is allowed by the tests we consider. In addition it is shown that some of the allowed 
models are distinctly inhomogeneous. Chapter^ is the most thorough, and deals exclusively with the spherically symmetric 
subclass of the models derived in chapter^. 
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Chapter 1 

Introduction and Review 



Since Hubble discovered the expansion of the universe and the homogeneous and isotropic expanding models of Fried- 
man, Lemaitre, Robertson, and Walker (FLRW) were accepted as the correct model of the universe, there has been 
relatively little consideration of alternative cosmological models. It is natural that 'mainstream' cosmology focuses 
on the understanding of the simplest acceptable models. However, it is also important to consider other possibilities; 
seventy years concentrating on one class of models is likely to lead to undue conviction in these highly special solutions. 
It is important to examine the assumptions on which cosmology is based, in the hope of improving our understanding 
of the universe. Indeed, it is essential that the assumptions can be tested wherever possible. This thesis is an attempt 
to do just that. 

The assumption I will investigate in this thesis is that we are geodesic, or freely falling, observers. There are a 
number of reasons this assumption has been used: firstly it simplifies things enormously; secondly, if one 'imagines' 
galaxies floating about in space, then it seems 'obvious' that they must be freely falling; one assumes that galaxies are 
like particles of dust (the classic billiard ball approach to a physical system) which, in itself, implies geodesic observers. 

As convincing as the geodesic assumption is, a standard FLRW dust universe cannot, on its own, satisfy the latest 
supernovae la (SNIa) results (Riess et ai, 1998, Schmidt et ai, 1998 and Perlmutter et ai, 1999), which imply that the 
rate of expansion of the universe is increasing; galaxies are moving apart faster and faster as time goes on. The simplest 
generalisation to dust, which solves this problem is the cosmological constant, or vacuum energy density: a concept 



which has been invoked and rejected as each new crisis is faced by cosmologists (figure 1.1). It was originally invoked 
by Einstein because of the belief at the time that the universe was static. It was only when Hubble discovered that the 
universe is expanding and the FLRW models were generally accepted that an alternative route was available. This led 
to the standard big bang cosmology. 

It is possible to achieve an accelerated expansion of the universe without invoking the cosmological constant, but 
this requires a large negative pressure: note that gravity effectively becomes repulsive whenever the pressure is large 
and negative enough, that is when 

H + 3p<0; (1.1) 
(whe re ji, p are the energy density and pressure respectively) which is precisely when the strong energy condition fails 



(see §5.1.1 and §6.2). This happens because /i + 3p is the effective gravitational mass. This is most elegantly shown 



in the Raychaudhuri equation (2.33), which is the fundamental equation of gravitational attraction. It says that the 



expansion rate 6 changes with time according to 

9 ^ -^{fi + 3p) + A (1.2) 



where A is the cosmological constant. In fact, it is easy to see (§2.5.1) that a cosmological constant is indistinguishable 



from a constant pressure; we may acknowledge this and identify the effective active gravitational mass with 

fi + 3p~2A (1.3) 

which shows why a cosmological constant may lead to an increasing expansion rate. 

Central to this thesis is an assumption we shall keep: the Copernican principle. There is no precise definition of the 
Copernican principle. In its weakest form it states that we are not at the center of the universe. This is not particularly 
useful for this thesis, so we will take a slightly stronger version: we are in a 'typical' location as observers in the 
universe. (This is sometimes known as the weak cosmological principle; Ellis, 1975.) Stronger still would be to say that 
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Figure 1.1: The variation in the accepted value of the cosmological constant over time, from Freedman (1999). 



all observers are equivalent, which would then force homogeneity. This is the cosmological principle (CP), and cannot 
be satisfied for inhomogeneous models. The CP follows from (either version of) the Copernican principle if we assume 
perfect isotropy about ourselves. With these assumptions we are lead to the standard homogeneous and isotropic model 
of the universe. We aim to show that the homogeneity of the universe does not follow from the Copernican principle, 
given that observations about ourselves are not perfectly isotropic. This allows the CP, and thus the standard model to 
be questioned. 

Although it essentially a philosophical assumption, the Copernican principle must be taken seriously (Ellis 1975): 
in order to reject an inhomogeneous cosmological model on the basis of its conflicting with the observed isotropy of 
the universe it is necessary to consider all observer positions and to show that for most observers in that spacetime 
the anisotropy observed is too large to be compatible with observations. We could adopt this view of the Copernican 
principle in this thesis, but, for simplicity, we will require consistency with observations for all observers - our results 
will thus be rather stronger than is strictly required by the Copernican Principle. In fact, it is only really possible to 
consider the weaker version of the Copernican principle when the spatial sections of the universe have finite volume (or, 
at least, when the number of observers, as measured by the integral over the spatial sections of the number density of 
particles - cf. §4.2 - is finite), otherwise the expression 'most observers' has very little meaning. For infinite universes 
the stronger version of the Copernican principle must be adopted. With our stronger definition we are thus prepared 
for all eventualities (although it will turn out that the models we consider are all 'finite'). 

Non-central locations are rarely considered in the literature in the analysis of inhomogeneous cosmologies however, 
owing to the mathematical difficulties that this usually entails. Having said that, Humphreys et al. (1997) made a study 
of Tolman-Bondi models from a non-central location and applied their results to a 'Great Attractor' model. Most other 
non-central analyses, however, look only at perturbations of standard FLRW models. 

Having adopted (a strong version of) the Copernican principle, the question then arises as to whether the observed 
isotropy of the universe, when required to hold at every point, forces homogeneity, thus validating the CP. Well, the 
nearby universe is distinctly lumpy, so it would be difficult to claim there is isotropy on that basis. However, the CMB 
is isotropic to one part in 10^. Together with the ECS theorem (Ehlers, Geren and Sachs (1968)), or rather, the almost 
ECS theorem of Stoeger, Maartens, and Ellis (1995) this allows us to say that within our past lightconc the universe 
is almost FLRW (ie, almost homogeneous and isotropic), provided the fundamental observers in the universe follow 
geodesies (that is, as long as the fundamental fluid is dust). What the almost EGS theorem achieves is to provide 
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support for the CP without the need to blithely assume the (near) isotropy of every observable: isotropy of the CMB 
alone is enough to ensure the validity of the CP (for geodesic observers). 

The ECS theorem is crucial to this thesis. It states that if all observers in an expanding dust universe see an 
isotropic radiation field, then that universe is FLRW. As the observed high isotropy of the CMB so well established, we 
can combine it with the Copernican principle as the starting point of this work. We generalise the ECS theorem to the 
case of an irrotational perfect fluid (ie, allowing for acceleration), to find a class of models which generalise the FLRW 
model. These models are a subclass of the inhomogeneous Stephani spacetimes with symmetry. We are left with a class 
of spacetimes which allow an isotropic CMB for all observers. We then study these models from all observer locations 
to show that these inhomogeneous models are acceptable given present observational constraints. They thus satisfy the 
Copernican principle while being inhomogeneous. 

A number of inhomogeneous or anisotropic cosmological models have been studied in relation to the CP. The 
homogeneous but anisotropic Bianchi and Kantowki-Sachs models (Kantowski and Sachs 1966; see Ellis 1998, §6 and 
references therein) have been investigated with regard to the time evolution of the anisotropy. It can be shown, for 
example, that there exist Bianchi models for which a significant phase of their evolution is spent in a near-FLRW 
state, even though at early and late times they may be highly anisotropic (again, see §6 of Ellis 1998 and references 
therein). Of the inhomogeneous models that arise in cosmological applications, by far the most common are the 
(Lemaitre-)Tolman-Bondi dust spacetimes (Tolman 1934; Bondi 1947). These are used both as global inhomogeneous 
cosmologies - probably the most important papers being Hellaby and Lake (1984, 1985), studying geometrical aspects, 
and Rindler and Suson (1989), Goicoechea and Martin-Mirones (1987), Schneider and Celerier (1999), Celerier (1999), 
and Maartens et al. (1996) investigating observational aspects - and also as models of local, nonlinear perturbations 
(over- or under-densities) in an FLRW background (Tomita 1995, 1996; Moffat and Tatarski 1995; Krasihski 1998; Nakao 
et al. 1995). See also Krasihski (1998) for a review. 

There has been some consideration of Stephani solutions (Stephani 1967a, b; see also Kramer et al. 1980 and Krasihski 
1983, 1997). These are the most general conformally flat perfect fluid solutions - and obviously therefore contain the 
FLRW models. They differ from FLRW models in general because they have inhomogeneous pressure, which leads to 
acceleration of the fundamental observers. D§,browski and Hendry (1998) fitted a certain subclass of these models to 
the first SNIa data of Perlmutter et al. (1997) using a low-order series expansion of the magnitude-redshift relation for 
central observers derived in Dsibrowski (1995), and found that they were significantly older than the FLRW models that 
fit that data. 

One unusual feature of the Stephani models is their matter content. The usual perfect-fluid interpretation precludes 
the existence of a barotropic equation of state in general (because the density is homogeneous but the pressure is not), 
although they can be provided with a strict thermodynamic scheme (Bona and Coll 1988) - it has recently been shown 
explicitly that they may be given a physically reasonable interpretation (Sussman 1999). Moreover, even individual fluid 
elements can behave in a rather exotic manner, having negative pressure, for example (c/. § ^.1.1| ). For these reasons, 
amongst others, Lorenz-Petzold (1986) has claimed that Stephani models are not a viable description of the universe, 
but Krasihski (1997, p. 170) argues rather vigorously that this conclusion is incorrect, as do we. Other cosmologies 
have sometimes been ruled out a priori because the stress-energy tensor does not behave in the correct manner, or 
for other reasons such as the lack of an FLRW limit - see Krasihski (1997) for a complete review. However, it has 
become increasingly difficult to avoid the conclusion that the expansion of the universe is accelerating, with the type la 
supernovae data being the latest and strongest evidence for this. It may be taken as evidence that there is some kind 
of 'negative pressure' driving the expansion of the universe. In standard FLRW models, this must correspond to an 
inflationary scenario, with A > 0, or a matter content of the universe which is mostly scalar field ('quintessence' - 
see Frampton 1999; Liddle 1999; Coble, Dodelson and Frieman 1997; Liddle and Scherrer 1998; see also Goliath and 
Ellis 1998 for a discussion of the dynamical effects associated with A). Either way, the real universe is behaving in a 
manner that is at odds with 'everyday' physics, so it is not appropriate to rule out Stephani models for exhibiting similar 
behaviour. 

We now briefly review FLRW models and then discuss the observational constraints primarily used in this thesis. 

1.1 The FLRW Models 

The standard model of modern cosmology are homogeneous and isotropic expanding FLRW models. They are param- 
eterised by three independent functions of time; {Hq, J7o, ^a}- The present day expansion rate is given by the Hubble 
constant, Hq; the density of matter today is given by Qq, which is normalised with Hq; and Q/^ represents a possible 
cosmological constant. 
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These models expand from a big ban^ into the universe's present state. The models may recoUapse after a finite 
time or expand indefinitely into the future: they are called 'closed' or 'open' respectively depending on this future fate, 
while the limiting case between the two possibilities is called 'fiat'. The different possibilities depend on whether the 
matter density and cosmological constant (in the combination fio + ^^a) are large enough to halt the expansion of the 
universe. The time since the big bang depends crucially on Hq. A significant problem facing cosmologists is an accurate 
measurement of these parameters. 

1.2 The Hubble constant 

The Hubble constant, Hq, represents the present day expansion rate of the universe. It is estimated, not surprisingly, 
by measuring the recession velocities of nearby galaxies, which move according to the Hubble law; 

V « Hod, (1.4) 

where v is the recession velocity and d is the distance to nearby galaxies. (This relation holds only in the 'local' universe.) 

Measuring Hubble's constant has proved to be extremely difficult in practice; despite 70 years searching it still 
would be difficult to claim anything better than 10% accuracy. In fact it is only within the past few years that the 
'factor of 2' uncertainty has been resolved. One of the biggest problems lies in calculating distance accurately: this 
becomes more difficult the further away galaxies lie. Distant galaxies are preferable because their random or peculiar 
velocities are substantially smaller than the Hubble expansion. The most common method of distance measurement is 
to use the observed magnitude. To infer distance from this one must obviously know the actual luminosity - and one 
must therefore use objects which have a narrow range of intrinsic luminosity which is independent of distance. A good 
example of such an object are Cepheid variables which pulsate regularly - the period of which is known to correlate 
to their luminosity. The most promising example of distant standard candles are the supernovae la (SNIa), which can 
outshine whole galaxies ~ and are thus observable to great distances. The maximum brightness, and subsequent light 
curve shape is know to correlate to luminosity. Other examples include spiral galaxies whose rotation velocity correlates 
with luminosity (the TuUy-Fisher relation). 

The main error in using these standard candles^ arises because they must be combined together to obtain Hq via 
the 'cosmic distance ladder'. One of the main sources of error in Hq lies in finding the period- luminosity relation for 
Cepheids, which is a consequence of errors in the distance to the Large Magellanic Cloud. 

These methods are finally reaching conclusion, albeit with different groups reaching slightly different mean values; 
finally, however, all the results seem to be consistent. At the 95% confidence level Freedman (1999b) quotes 57 < Hq < 
85 km s^^ Mpc^^. Despite this many would argue that iJo = 65 ± 5 km s^^ Mpc^^ is a safe bet; see Trimble and 
Aschwanden (1999) for an entertaining discussion of this. 

There are promising new methods for measuring Hq from observations of objects at large distances. These include 
time delay in gravitational lensing events, and use of the Sunyaev-Zel'dovich effect from the X-ray emitting gas in 
clusters affecting the cosmic microwave background (CMB) spectrum. 

Measurements of Hq from gravitational lensing make use of distant variable sources which are lensed, producing 
multiple images. The different paths the light travels will result in a time delay effect, which is measurable. Measurements 
of the angular separation of the two images then allow Hq to be determined. The problem with this is that the mass 
distribution of the lensing galaxy will not be known independently. This, combined with the difficulty of finding a 
suitable system makes gravitational lensing a test for the future (Freedman 1999a). 

The Sunyaev-Zel'dovich effect is the scattering of the CMB photons off electrons in the X-ray emitting gas of 
rich clusters (Sunyaev-Zel'dovich 1969). This results in a change in the CMB spectrum. The X-ray ffux is distance 
dependent, but the SZ temperature ffuctuations are not, therefore Hq may be determined. As with the gravitational 
lensing measurements, it is still early days, with large uncertainties, and only a few clusters. 

1.3 The Density Parameters Qq and ^Ia? and the Curvature 

In the standard model the energy density, curvature, and cosmological constant are all related by the normalised (with 
respect to Hq) parameters, {no,flk,^A} by 

no + nA + nk = i. (1.5) 

^Some models may 'bounce' indefinitely depending on the model parameters. 

^Actually, a standard candle is an object with a very narrow range of brightness, but there is no need to labour this point here. 
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We see that determining any two of the three will give the third. Determining these will determine the origin and fate 
of the universe. 

The density and curvature of the universe remain in question, and this uncertainty has led to the dark matter 
problem. The main problems here lie in the fact that the luminous matter density of galaxies nearby give a value of 
ilo ~ 0.01 (Big Bang Nucleosynthesis also gives similar, but slightly higher, results - see Wainwright and Ellis 1997, or 
Olive, 1999, for details), while dynamical studies of the same galaxies suggest that flo S 0.3 - ie, they are gravitationally 
bound by a far greater mass than we see. This means that most of the matter in the galaxies is 'dark'. There have 
been plenty of suggestions as to what it may be, see Carr (1994) for further details. In addition to this the latest SNIa 
measurements suggest that the deceleration parameter qo = ^ilo — < which necessarily implies that Q\ > - see 
figure 2.7. If the SNIa results are correct, then there is no way, in the standard model, to avoid invoking a cosmological 
constant or scalar fieldQ model. 

Measurements of these parameters from 'local' {z ^ 1) observations (eg, SNIa) , and from observations of the CMB 
(z ^ 1000) complement each other (Tegmark, Eisenetein and Hu, 1998a; Tegmark e t ai, 1998b; Eisenstein, Hu and 
Tegmark, 1998a,b; see also Tegmark, 1999; White, 1998; Efstathiou et at, 1 998; a nd § ^.8.2| ). The SNIa determine the 
quantity ^Iq — ^Ia , while CMB results give JIq + ^a (ie, the curvature) . See § 2.8.2 for details why. This means that the 
two data sets will give all of {fto, fife, ^^a}- 



Supernovae Results Riess et al. (1998), Schmidt et al. (1998) and Perlmutter et al. (1999). These all suggest that 
A is distinctly non-zero and positive. Accurate measurements of the SNIa constrain the deceleration parameter because 
they measure the change in the expansion rate with distance. They are observed at roughly half the age of the universe. 

The primary errors in using SNIa for the determination of the density parameters are similar as for iJo. In addition, 
the SNIa may evolve (as has been suggested by Drell, Loredo, and Wasserman 1999) ie, intrinsic brightness changes 
with distance. This simply amounts to not knowing precisely enough how intrinsically bright the objects in question 
are at the time of emission. The SNIa have caused such a stir initially, because they are thought not have this problem. 

Perlmutter et al. (1999) quote their best fit as fio = 0.73, Ha = 1.32; while their best fit results assuming a flat 
universe is = 0.28, f^A = 0.72 with errors of ±0.14. Regardless of the type of fit its claimed that A > at the 99% 
confidence level. The other group achieve similar results. 



CMB Measurements The CMB was discovered in 1965 by Penzias and Wilson, and interpreted as a 'relic' of the 
big bang by Dicke et al. (1965). 

The CMB is observed today to be a blackbody at a temperature of To = 2.734 ± O.OIK, with a dipole moment of 
Ti = 3.343 ±0.016 x lO-^K and quadrupole moment as big as T2 = 2.8 x lO'^K (see Mather et al, 1994; Partridge 1997 
for details and references). It was emitted at a time when the radiation was no longer hot enough to keep Hydrogen 
ionised, causing it to decouple from matter, which happens at T,^^ ~ 3000K. At this time, the small perturbations 
(created by inflation) which were present in the universe left their mark on the CMB surface. The baryons (electrons 
and protons) fell into the potential wells created by the small perturbations. Because the baryons were still coupled to 
the photons, the photon pressure acted as a restoring force against the motion of the baryons into the potential wells. 
This led to acoustic oscillations. These oscillations may be decomposed into their Fourier modes, with the first acoustic 



peak at / '-^ 200/v^Io+TIa - see figure 1.2 



Because the position of the first Doppler peak depends essentially only on the curvature, the CMB is a good test 



of + f^A- From figure 1.2, we can see that the curve favours a flat universe - the first peak is around I ^ 200. The 



estimated curvature is flo + ^a = 1 ± 0.2. See Turner (1999), Straumann (1999), or Rocha (1999). 

This is still a preliminary result, but the issue should be settled soon with some new satellites being launched in the 
next year or so (eg, MAP, Plank, etc.). These results are not particularly relevant for the work contained in this thesis. 

The SNIa and CMB measurements may be combined to give more accurate information on the density parameters - 



2.8.2. They can also be used to find the 'equation of state' of the universe; see Perlmutter et al. (1999), or 



Efstathiou (1999). 



1.4 Age 

There are two ways to determine the age of the universe. Determinations of the ages of 'old' objects such as globular 
clusters in the local universe provide a lower bound to the age of the universe. Alternatively, given a model of the 
universe, a measurement of Hq provides an estimate of the age of the universe. Both are plagued with uncertainty and 

^Also known as a quintessence model - see later. 
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Figure 1.2: The CMB multipoles as a function of temperature. 

From Max Tegmark's home page: http://www.sns.ias.edu/~max/foregrounds.html. 



there has been considerable disagreement until very recently, with globular clusters being significantly older than the 
age as implied by the expansion rate. 

From Hq we can obtain a range of ages from 8 Gyr for an Einstein-de Sitter model, all the way up to 16-17 Gyr for 
a low density flat model, both depending on the value of Hq. Given a set of parameters {Ho,flQ,nj\} the age can be 
calculated simply from 

^ dz 



^-1 ''-I lm(YT7y ^'-'^ 
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where H{z) is given by equation ( 2.13l| ). 

Measuring the ages of globular clusters is quite a complicated process; it is based on a number of complementary 
methods. It combines estimates from stellar models, with separate evaluations of the stars in the clusters turning off 
from the main sequence. As stellar models are not perfect uncertainties of about 7% arise. The largest uncertainties 
in the cluster ages arise from estimating their distance - and hence their intrinsic luminosities which are required to 
calibrate the stellar models. Distance estimates are made from parallax measurements, and from distance indicators 
such as RR Lyraes. As in estimating Hq, calibration of these stars relies on knowing the distance to the Large Magellanic 
Cloud, which again introduces considerable uncertainties. 

The ages of the globular clusters have traditionally been quite high - certainly too high for a flat A = model, and 
led to the 'age problem' as it was known. However recent recalibration of RR-Lyraes has led to a considerable reduction 
in the ages of these to about 11-14 Gyrs (Chaboyer et ai, 1998; Krauss, 1999); also recent parallax measurements using 
the Hipparcos satellite confirm these results (Reid, 1997; Gratton et ai, 1997). For a high value of Hq or a high density 
model, the margin of error between cosmological and globular cluster ages is still very small. See Freedman (1999b) for 
discussion and references. 

The age of the universe is important for this thesis because the Stephani models we consider here were initially put 
forward for their 'naturally' high age, and as a solution to the age crisis (D§,browski and Hendry, 1998). We constrain 
the Stephani models later using the age. 



Chapter 2 

Relativistic Cosmology 



The goal of cosmology is to find a model tliat best describes the universe on any particular scale, preferably on all 
scales. The framework which is usually used, and will be used throughout this thesis, is that of General Relativity 
(GR) as discussed in eg, Wald (1984); Stephani (1990); Misner, Thorne and Wheeler (1971); Schutz (1990); Hawking 
and Ellis (1973); that is to say spacetime is a manifold, which has a geometry described by a Lorentz metric gab, and 
associated connection F"^^, containing matter whose physical properties are described by the energy-momentum tensor. 
Tab- The curvature of the spacetime is felt via the Riemann tensor, Rabcd, which gives the lack of commutivity of 
derivatives when parallel transporting vectors in a curved manifold (Schutz 1980); 

(VaVb - VbVa) = i?abcde'', (2-1) 

for any vector field ^°; these are the Ricci identities. The curvature of spacetime and the matter content interact via 
Einstein's field equations, 

Gab = Rab — \Rgab = Tab — J^9ab, (2.2) 

where Gab is the Einstein tensor, Rab and R are the Ricci tensor and scalar, and A is the cosmological constant.]^ 
These equations give local conservation of energy from the geometry of the spacetime. The vector formed by taking the 
4-divergence of Tat, ^bT""^ represents the local creation or loss of energy. The Bianchi identities, 

^[aRbcf' = 0, (2.3) 



when twice contracted (over the ac and ed indices) imply VfcG" — 0, from (2.2) so that, 

VbT"'' = 0, (2.4) 

so that energy is locally conserved. 

For any model of the universe, we must bear in mind that we have to extract observational predictions from it; (2.2) 
are not particularly intuitive. Rather than just start from a metric and see what we get, it is desirable to break (2.2) 
down into a simpler (or at least more intuitive) set of equations while retaining their covariant character, and at the 
same time introducing quantities that are directly measurable. 



2.1 The Covariant 1+3 Formulation of Fluids in GR 



In this thesis, we don't make much use of the 1+3 formalism (Ehlers 1993), but it is necessary for the results of §2.3 



and §3.2, and also for a proper understanding of modern relativistic cosmology. 

For a complete cosmological model one must specify not only a metric defined on a manifold, but also a family of 
fundamental observers. These worldlines may be used to represent galaxies at late times (where the galaxies sit on the 
worldlines) or radiation at early times. It is usually assumed however that at any particular epoch there is, on average, 
one dominant or fundamental congruence of worldlines, eg, that defined by the CMB. The combinations of such a set 
of fluids will be described briefly later. These worldlines will have a velocity 

"" = i^' («) 

^In this chapter units where SttG = c = 1 will be used; in later chapters, where observational quantities are important, they will be put 
into the equations. 
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where r is the proper time measured along the worldlines, such that UaU°' = —1; ie, it is timelike. This velocity field is 
very important and its properties (together with a matter description) can be used to give invariant definitions of many 
cosmological models. 

This velocity field can be used to look at tensors along these worldlines, and orthogonal to them. This is the '1+3' 
splitting of the spacetime into 'temporal' and 'spatial' parts; it is also covariant because can be defined uniquely and 
without any coordinates which would be required for splitting the spacetime into space+time. 

Given u"' we define a projection tensor by 

hab " gab + UaUb, (2.6) 

which defines a '3-metric' orthogonal to the congruence, and satisfies 

ft/ = 3, habv'^Q; (2.7) 

ie, it is orthogonal to u", it therefore projects things into the instantaneous rest space of the congruence. 
We also make use of the projected alternating tensor (a generalisation of the Levi-Civita symbol) 

Tlabc = VabcdU'^, (2.8) 

where 



Vabcd '^~V\9\5\a5\5\5\^ (2.9) 

is the spacetime alternating tensor (the 4D Levi-Civita symbol). 

We use angle brackets to denote the projected, symmetric, trace-free part of a 2nd rank tensor, 

F^ab) M (h^^h.f - ^habh'") F,d, (2.10) 

(and for a 1-form or vector, K(^a} = ha Kb) so that any projected 2nd rank tensor has the irreducible covariant decom- 
position 

Fab = \Fhab + riabcF" + F(„6> , (2-11) 

where F = Fcd.h"'' is the spatial trace, and Fa = ^rjabcF'"' is the spatial dual vector of the antisymmetric part of Fab- 
In the 1+3 covariant formalism, all irreducible quantities are either scalars, projected vectors or projected, symmetric, 
trace- free tensors. 

With this projection tensor and velocity field, we can define two derivatives; firstly, differentiation along the fluid 
flow - a time derivative, denoted by a dot. 

Fab M U^ycFab-, (2.12) 

and also the derivative projected orthogonal to the flow lines - a spatial derivative, 

VcFab if h\h\h\VdF,f. (2.13) 

This new tensor is entirely orthogonal to u"; contraction of any index with is zero. Va is no longer a proper 
3-dimensional covariant derivative however, because in general, 

V[„V6]Q^0; (2.14) 

which means that the congruence is no longer orthogonal to the spacelike surfaces. 

The covariant derivative of any scalar can be split into orthogonal and parallel parts; 

VaQ = V„Q-Q«a; (2.15) 

while derivatives of vectors and tensors can be split into their irreducible covariant bits, see Maartens, Gebbie and 
EUis (1999) for the equations. 
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2.1.1 Kinematics of the Fundamental Congruence 

The derivative of the fundamental congruence VaUf, can be spht into irreducible quantities. These are: 
The expansion, 

0^'Va^i" = VaU^ (2.16) 

which is the trace of the derivative of u° and also the 3-divergence of the congruence. It represents the volume 
rate of expansion of the fluid elements. In general, one can associate a fundamental or average length scale by 

I'^J (2.17) 

which describes the volume change of the fluid. 
The acceleration, 

ii" Mu^'Vbu'', (2.18) 

is the time rate of change of u"; it describes motion of the flow moving under forces other than gravity alone. This 
comment can be understood by thinking about a Schwarzschild black hole: if someone sits at a constant proper 
distance from the black hole, then the force of gravity will pull them into the center. In order to stay at a constant 
proper distance from it, then there must be some non-gravitational force pushing outwards. Acceleration is zero 
if and only if the flow is geodesic (freely falling observers). 

— The shear, ^ 

CTab^Vi^aUb), (2.19) 

is the rate of shearing of the congruence; ie, the trace-free symmetric part, which describes how the congruence 
will distort in time. While it doesn't affect the volume change of the congruence, the relative distances of objects 
will change because of the presence of shear. Its magnitude is defined by cr^ ^ ^Cafof and cr = 4^ (Tab = 0. 

The rotation, 

iOab^y[aUb] (2.20) 

describes the rate of rotation of the congruence; it is the orthogonally projected anti symmetric part of the 
derivative of the flow. Objects will change their position in the sky because of rotation. It also has magnitude 

Lu^ icjaf,cj°^, with a; = <^ ujab = 0. We also define a rotation vector by = ^r]°-'"^LJbc- 

The covariant derivative can now be decomposed as 

yaUb = -UaUb+^Ohab + Crab+^ab, (2-21) 

giving a covariant representation of the changing congruence in terms of invariantly defined quantities - which also give a 
physical breakdown of different aspects to the kinematics of the congruence. This equation is of fundamental importance 
to cosmology - and relativity in general - and was first given by Ehlers (1961), translated recently as Ehlers (1993). See 
also Ellis (1971). 

We are now in a position to understand and expand ( ^.14| ), and a derivation is probably useful: 

V[,Vb]Q = -uJabQ (2.22) 

for any scalar field Q. This result was first presented in Bruni, Dunsby, and Ellis (1992). It shows explicitly that the 
fundamental congruence will define a spatial '3-metric' orthogonal to the flow line if and only if the rotation is zero, 
because that is when V defines a proper covariant derivative. 



which implies, using (2.20) 



CHAPTER 2. RELATIVISTIC COSMOLOGY 



11 



2.1.2 The Energy-Momentum Tensor 

Regardless of the particular matter present, any energy-momentum tensor can be decomposed with respect to the chosen 
fundamental congruence in the following manner: 

Tab = f^UaUt + phab + 2q(^aUb} + T^ab, (2.23) 

with each of the quantities having a physical interpretation: 
The relativistic energy density 

/X M u'^u'^Tab = u^u''Gab " A; (2.24) 

— The isotropic pressure 

p M Ih'^'^Tab = i/j^'Gab + A; (2.25) 

— The energy flux, or relativistic momentum density 

qa'^T(a)bu\ (2.26) 

which can be interpreted as the heat flow relative to it satisfles gait" = 0; 
The anisotropic pressure 

TTab'^T^ab). (2.27) 

which is trace-free, 11^°' = 0, and u°'Trab = 0. 

2.1.3 The Sphtting of the Weyl Tensor 

The trace-free part of the Riemann tensor is the Weyl tensoi]^, Cabcd, which describes the free gravitational field - tidal 
forces and gravitational waves - and to some extent describes the null structure of the spacetime. It is defined by 

Rabcd ~ Cabcd + 9a[cRd\b ~ \Rga[c9d]b- (2.28) 

In the 1+3 splitting of the spacetime, the Weyl tensor can be split too, into 'electric' and 'magnetic' parts: 

Eab ^ CabcdU'u'' (2.29) 
Hab ^ ^rjacdeC"'bju'uf. (2.30) 

Both of them are symmetric, trace-free, and orthogonal to w". 

2.2 Splitting Einstein's Equations 



If we take Einstein's equations (2^), the twice-contracted Bianchi identities, the Bianchi identities (2^), or using ( 2.28 ) 

VCahcd = V[a [ — Rb]c + ^R9b]c) (2.31) 



and the Ricci identities (2.1) for the congruence w", and separate out all the independent parts, we arrive at a set of 
evolution and constraint equations describing the structure of spacetime. Not all of these are used in this thesis, but 
they are included for completeness. 

Evolution: 

The energy conservation equation 

fi + e{^i +p)+ VaQ" = -^UaQ" - f^abTT'"'; (2.32) 

shows that, for a perfect fluid, the expansion correlates directly to the change in energy density along the flow 
lines; 



Also known as the conformal tensor, due to its invariant nature under conformal transformations; see §2.3. 
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— The Raychaudhuri equation (expansion evolution) 



i6»2 = VaM° + ilaU" " 2(7^ + 2uj^ - + 3p) + A, (2.33) 



which is the equation of gravitational attraction. This shows that a positive cosmological constant, or rotation, or 
acceleration will make the expansion rate increase - which is the accepted scenario at present, see later - whereas 
shear will slow the expansion rate. In a sense this is obvious; one would expect the 'rotation of the universe' 
to increase the expansion rate, and similarly a large negative pressure (or a positive cosmological constant - 
see ( 2.6lD ) will do the same; 

The momentum conservation equation 

q{-) + leq- + V> + VfcTT"^ = -a'^^b - (m + P)"" - utn''' - v'^'^^OJ^qc, (2.34) 
which shows that a perfect fluid can have acceleration only if there are spatial pressure gradients present; 

— The shear and rotation evolution equations 

Cj^-^^ + pUj"" = ^J]"''"^bUc + fT^^'cjfc, (2.35) 
Cr(°''> + |6l(7'^f' ^(''u'') + ui^ub) _ ^{a ^^b)c _ ^{a ^b) _ ^ab _^ l^ab. (2.36) 

the second equation shows that a fluid with oja = (Jab = Ua = will have the electric part of the Weyl tensor 
proportional to the anisotropic pressure. 

— The evolution of Eab and Hab, 

j^(ab) ^ g^ab _ ^cd{a^ ^jjb)^ ^ l^{ab) ^ _ 1 jy(a^6) _^ _^ ^yab] _^ ^^{a ^ ^^b)c _ l^b)c^ 

-ii^-q"^ - leT:-" + 7^-'^- {2urH\ + lo,[e\ + \^'\]) , (2.37) 

- r?^^^'^ {2u,e\ - loM\ - \o'\qa) . (2.38) 

Constraint: 

The shear and vorticity divergence equations 

VaW° = WaW" (2.39) 

Vfccr"'' = IV^e ~ r;"''" {VbOJc + 2ubOj)j - (2.40) 



The second equation is crucial in §3.2, where it is used to set energy flux equal to zero in a QCDM model with an 
isotropic radiation flcld. 

The 'curl a' equation 

Ty^'^^'^VcCr^^d - H"^ = 2u<"w''> + V^'^w^^ (2.41) 
The Eab and Hab divergence equations 

Vfei;'^'' - iijJbH''^ = ^[vv - eq''] + i[cr°Vfc - VfcTr"''] 

+ 77°^^ {cJbdH\ - Iwbqc) , (2.42) 

- V"""' Uybqc + (Jbd[E\ + . (2.43) 
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2.3 1+3 Splitting Under a Conformal Transformation 

A conformal transformation is an angle preserving transformation that changes lengths and volumes. The importance of 
these types of transformations lies in the fact that, under a conformal transformation, the null structure of the spacetime 
is preserved: indeed, it trivially follows that the causal structure is preserved. We also have the important property 
that the the Weyl tensor, C^^^'^) invariant (note that one index must be raised) so that a conformal transformation 
will introduce no tidal forces or gravitational waves; that is, a conformal transformation will only introduce 'non- 
gravitational' forces and matter into the new spacetime (by changing R^t and thus the matter tensor Tab via Einstein's 
equations) . 

We perform the conformal transformation 

gab = e^^gab, u" = e-^u\ Ua = e^iia] (2.44) 

where Q > is an arbitrary function, u° is a velocity vector with respect to gab'- gabu'^u^ — Uau'^ = —1; and u"" is 
the conformally related (parallel) velocity vector, and is normalised with respect to gab'- gab'u°"u^ ^ ^1-0 "^h*^ covariant 
derivative of any one-form field Va transforms as 

VaWfc = VaWfc - 2Qi^aVb) + gabQ^Vc, (2-45) 

where Qa = Q,a = ^ aQ — Qua- The expansion (9 = VaU°), acceleration (u° = u^Vhu"'), rotation {ujab ~ V[aUb])i and 
shear {uab — {a^b)) of the two velocity congruences are related by: 



e 


= e'^ie^SQ) 


Ua 


^ Ua - ^aQ 


l^ab 


= e^'^UJab 


Cfab 


= e~^aab- 



(2.46) 

The equation for the acceleration corrects equation (6.14) of Kramer et al. (1980). These show that a conformal 
transformation mat induce acceleration and expansion into the new spacetime, by not shear or rotation: in particular, 
a conformally flat model must have shear and rotation vanishing (as in eg, the Stephani models considered later). 
With respect to gab, a dot denotes differentiation along the fluid flow - a time derivative, Fab = u'^^cFab', and Vq is the 
derivative projected orthogonal to the flow lines - a spatial derivative, WcFab = h'^ Jf^ a^^"^ dFef , where hab — gab + UaUb 
is the usual projection tensor. 

The Einstein tensor transforms as (see Wald 1982) 

Gab - Gab - 2VaQb - 2QaQb + gab [2V,Q= - Q^] , (2.47) 

where Gab is the Einstein tensor of gab, and ~ QaQ'^- For clarity with the above, we decompose derivatives of Q 
into time and space derivatives: 

Qa = ^aQ ~ QUa, 
'^bQa = VaVbQ + UaUb [q - il^V cQ^ - Q {\6hab + CFab " ^^ab) 

+2u(a [-Vb)Q + i0Vfc)Q + (a^f + c^,)^) V,q] . (2.48) 
We also write Tab — Gab, and Tab ~ Gab as general fluids, both with respect to u°: 

Gab = I^UaUb + Phab + 2q(^aUb) + T^ab, (2.49) 
Gab = H^aUb + phab + ^q^a^b) + TTab', (2.50) 

where {p,,p, qa,'!Tab}, and {^J-,p,qa,^^ab} are the energy density, isotropic pressure, heat flux, and anisotropic pressure of 
Gab and Gab respectively. 



•^When performing conformal transformations, confusion can arise over which metric to use when performing contractions; this will be 
avoided here by always using g^b- ie, v'^Va = gab^°"u^- 
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manner: 



fi = u^u'Cab = e^'^fi - 3g g - -0 - 2VaV°g + v^gv^g, 



3 



We can decompose Gat given by ( ^.47 ) into the fluid variables in (2.50) by using ( 2.48] ) in the following covariant 

2„\ ^ ^2.51) 

(2.52) 

(2.53) 
(2.54) 



qa = -u''G^a)b = e'^Qa - 2V aQ + 2 



2g + ^VaV"g-^Vagv'^g 

+ 2M,v=g, 
'-g) v„g + 2 (a/ + c^/) Vfcg 



T^ab = G(^ab) = ^ab + '2'Qo ab - '2'^ (o^ b)Q - '^'^(aQ^b)Q 



(aV6)l 



2.4 Relativistic Thermodynamics 

Any realistic cosmological model must include some sort of thermodynamic scheme. This means that we expect the 
laws of thermodynamics to hold throughout the evolution of the universe. 



Together with the conservation equations (2.4) the energy-momentum tensor (2.23), which lead to the equations (2.32) 
and (2.34), there are a number of other conservation equations which one may or may not impose upon the fluid. See 
Ehlers (1971) for derivations and a discussion of kinetic theory in GR. 

Firstly, we normally assume that the comoving number density is conserved; ie. 



V,(A^'^) EE 0, 



(2.55) 



where iV° = nu"" and n is the number density of particles - see Krasinski (1997), or Maartens (1996). Now, ( 2.55| ) is 
equivalent to the condition 

n + 9n^0 ^ = const., (2.56) 

where (|1^ was used to demonstrate that the number of particles in a comoving volume is constant. 

However, in a more realistic model this may not be true at some stages of the evolution, see for example Gunzig 
et al. (1997) where a particle creation rate is introduced to model creation of radiation due to inflation. Other general- 
isations include adding a viscous pressure term; eg, Coley, van den Hoogen, and Maartens (1996), or Maartens (1995). 



The First Law of Thermodynamics 

The first law is the Gibbs equation which applies in equilibrium: 

TdS = d(^) +pd 
where S is the entropy density, and T is the temperature of the fluid; or 





da;" + pVa 


■f" 


-71. 




n 



dx"- 



(2.57) 



(2.58) 



which becomes, upon dividing by an increment of proper time along the congruence, and substituting from ( ^.56 ) 



(2.59) 



which shows that = for a perfect fluid (cf, ( ^.32 )): ie, entropy per particle along a particular flow line remains 
constant. 



The Second Law 

Entropy flux is generally defined as a vector; 



= Sr 



(2.60) 
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where i?" represents all dissipative processes, and S and T are related by (( |2.57 ) - ie, they are scalars defined in local 
equilibrium). In the simplest cases i?" is assumed zero, or equated with the energy flux. The second law states VaS'" > 0, 
where equality applies in equilibrium (where i?" — 0). 

Equations ( 2.55| ), and (2.57) comprise a thermodynamic scheme which implies that the right hand side of (2.57) has 
an integrating factor - see Krasiiiski, Quevedo and Sussman (1997). These are additional constraints on cosmological 
models. (This is used in §3.1.3.) 



2.5 Fluids in Cosmology 



Equation (2.23) is a general equation describing any fluid. Different types of fluid have different energy-momentum 
tensors, which we will now discuss in order of importance for cosmological models. 

The most important fluids used in cosmology are dust solutions, often with a cosmological constant. The simplest of 
these are the FLRW models, which are generally believed to represent the universe on a suitably large scale; perturbations 
of these models give a more realistic representation of the universe. However, simple does not mean correct and therefore 
other models have been studied to give a different view of the universe; among these are the dust Lemaitre-Tolman (LT) 
models and Bianchi models - the latter often used because they are 'close to FLRW for some suitable length of time, 
and therefore give a better understanding of FLRW models themselves. 

Generalisation of a simple dust solution to a perfect fluid is a necessary, but not always simple, task if one wants 
a model of the universe before recombination; most solutions of Einstein's equations which are used as cosmological 
models are perfect fliuds (in particular, a scalar field can always be written as a perfect fluid). 

Here it will be useful to note that a cosmological constant can always be absorbed into the energy-momentum tensor 
by redefining the density and pressure. 



Gab 



Tab + J^9ab 

flUaUb + Khab - ^UaUb 



{fl - A)UaUb ■ 



Ahab; 



(2.61) 



ie, a perfect fiuid with constant pressure. A cosmological constant does not need to be studied separately in the other 
perfect fiuid cases. 

2.5.1 Perfect Fluids and Dust Models 

A pe rfect fluid occurs when iTab = = 0; in the case of dustg ('CDM') we also have p — 0, which necessarily implies 
(by ( 2.34 )) that the fundamental congruence is geodesic; = 0. These models are implicitly in equilibrium; that is, 
the dynamics are completely reversible and no heat is generated from friction or anything else (5 = 0). 
From the energy conservation equation (2.32), we get the simple relation for the expansion. 



and similarly for the acceleration ( 2.34 ), 



fi+p' 



(2.62) 



(2.63) 



However, the Raychaudhuri equation remains unchanged - emphasising that it is an equation relating the important 
kinematic quantities of the flow, rather than the more complex anisotropic pressures, energy flux or Weyl curvature. 
We can substitute ( 2.62] ) and (2.63) into the Raychaudhuri equation ( p. 35 ) to get a second order differential equation 
for the matter quantities. 

From the Gibbs relation ( p.57[ ) (which defines the temperature) the entropy of a perfect fiuid cannot change along 
the congruence (cf, 2.59 ); 

5 = 0. (2.64) 



The adiabatic speed of sound is given by 



(2.65) 



■'Also known as an incoherent fluid. 
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which must be less than the speed of hght. 

Perfect fluids can be split into different types, according to their dependence of /i on p. 

A Barotropic Equation of State 

Generally a useful simplifying assumption to make upon the fluid is that of a barotropic equation of state (EOS). It is 
a fairly ad hoc assumption which simply requires that the pressure depends only on the density: p = p^ji). It is adopted 
usually for mathematical simplicity rather than for any physical reasons^. This functional dependence can in principle 
take any form. There are a few results which follow from this simple EOS. 
From ( 2.63D we get 

V[>fc] = ^ 77"'"V6?ic = 0, (2.66) 

so that rotation is conserved, 



leu" 



cf (HH) 



For these types of cosmological models the entropy is a global constant; (from 2.57), 
which implies V aS = 0. This is an isentropic fluid - an unrealistic situation. 



(2.67) 



(2.68) 



Dust: p = {). 

This is the most common fluid form used for a cosmological model, 
and thus eliminating many terms in Einstein's equations. 

In this case the energy density scales inversely with volume 
of a model by 



Its simplicity lies in making the motion geodesic. 
In general, one defines the characteristic length scale 

(2.69) 



hence, by (2.62) we have 



3£ 

7 



-T = ^ 



Dust is a particularly simple cosmological model to deal with 
dust moving under gravity alone. It requires that the random velocities of the galaxies be negligible 
is very low, so may not be used in the radiation dominated era. 



r'^. (2.70) 

It models the late universe of galaxies as 'particles' of 

ie, the temperature 



7-law equation of state 

This is a classic simple type of perfect fluid, which can describe a number of situations. We have an equation of state 
of the form 

P=(7-1)a*: 1<7<2, (2.71) 

where 7 is constant. 7=1 corresponds to dust, described above. Again, from (2.62), we have p. ~ . The other 
important case is that of radiation, 7 4/3. 

If 7 = 2, then we have a 'stiff fluid', in which the speed of sound equals that of light. This is, in a sense, the equation 
of state of the 'ether'. It is usually discounted as a physically reasonable form of matter as it fails the energy conditions, 
but some scalar field models have an effective equation of state of this form. A cosmological constant may be accounted 
for by 7 = 0, ^ = A. 

In the standard model (FLRW), we have a multi-fluid description: we have a mixture of (non- interacting) matter 
and radiation, with the radiation dominating at early times. If we have the characteristic length scale I in an expanding 
universe, then radiation density falls off as but the matter density of the dust decreases as t'^: hence, if we proceed 
forward in time from the initial singularity it — 0), then radiation will dominate early on, but matter will take on the 
dominant role later. 



^J. Ehlers, private communication. 
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2.5.2 Scalar Fields 

Scalar fields are fashionable at the moment for their use in describing inflationary scenarios in the early universe, but 
they may be pre sent today (quintessence). There is no need to go into any physical detail here, but the basics are 
required for §3^. A scalar field has an energy- momentum tensor of the form 

Tab ^ ^a^l^yb^b ~ gab (^Ve</)V=0 + 1/(0)) , (2.72) 

where V must satisfy the Klein-Gordon equation; 

r{cj)) = V^Va^. (2.73) 
If Va^ is timelike then we may define the velocity field 

(2.74) 



so that Tab has the form of a perfect fluid with effective energy density and pressure 

/x = -iv,(/.v'^</)-Hy(</)), p = -iv,(/.v"</)-y(</)). (2.75) 



2.5.3 Velocity Fields 

In general there is no canonical definition of the 'correct' fundamental velocity field to use and the energy-momentum 
tensor will have a different interpretation depending on the chosen congruence. That is, it is not obvious whether 
fundamental observers should be identified with the natural congruence in, say, a perfect fluid, or with some other 
velocity field (Coley and Tupper, 1983, 1984, 1985). This is fundamental to the 1+3 formalism: the decomposition 
above is velocity field dependent. 
If we take a perfect fluid 

Tab = pLUaUb + phab (2.76) 

and consider this with respect to another congruence defined such that 

u'^ MY^u" + V'') (2.77) 



with v'^Ua = (ie, f ' is the 3- velocity relative to u") and T = Xj^JX — - basically a Lorentz boost - then Tab will have 
the general form ( 2.23 ), with the relative energy densities etc., related by 

fi = fi + T'^v^{(i + p) 

p = p + lr^v^ifi + p) 

Qa = T'^{fl+p)Va 

TTab=T^{tl+p){vaVb-lv^hab}. (2.78) 

(See Wainwright and Ellis, 1997.) We can see that for a perfect fluid a relative energy flux will always be introduced by 
a change of velocity field. This is required for 

In the case of a general fluid with respect to another velocity field, the transformations are given in Maartens, Gebbie 
and Ellis (1999). 



2.5.4 More General Fluids 

The perfect fiuids considered up to now are clearly unphysical: the entropy is constant and there is no frictional heating. 
While this is not a problem if all we want is a rough and ready model to work with (eg, the standard model), we require 
something a little more involved for a decent model of the universe: thermodynamic processes in the real universe are 
(probably) not reversible. Unfortunately, there doesn't seem to be a well formulated theory of relativistic dissipative 
fluid mechanics, and there only seem to be a couple of cases occasionally used, such as bulk viscosity to describe 
inflation using the truncated Israel-Stewart theory of irreversible thermodynamics (Israel and Stewart, 1979, 1980); see 
Maartens (1996) for a review. See also Gariel and Le Denmat (1994). 
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2.6 The CMB Anisotropies 

In we discussed the observed characteristics of the anisotropies in the CMB radiation. Relativistic cosmology must 
explain the effects which cause these anisotropies before decoupling (ie, perturbations of the spacetime structure, density 
fluctuations etc.) and effects that occur after the last scattering surface which affect the size of the anisotropies as the 
photons travel to us (eg, lensing by structure, the Rees-Sciama effect - the changing redshift of the CMB photons as they 
pass through varying gravitational potentials; Rees and Sciama 1968 - interstellar dust - Finkbeiner and Schlegel 1999) 

The first attempt at a relativistic treatment was by Sachs and Wolfe (1967) by integrating the null geodesies in 
a perturbation of a flat A = FLRW model to examine the redshift function along the null rays. This has been 
done in a gauge-invariant and covariant 1+3 way by Challinor and Lasenby (1998, 1999) for scalar perturbations. These 
approaches have been generalised by Maartens, Gebbie and Ellis (1999) to include some more non-linear effects - see also 
Maartens (1999), Challinor (1999), Gebbie and Ellis (1999). The essential steps, after choosing a suitable velocity field, 
are to decompose the photon distribution function, the collision term in the Boltzmann equation and the temperature 
fluctuation into spherical harmonics to get the multipole moments for an inhomogeneous spacetime. 

The covariant and gauge invariant results are presented in the above references, and are a bit of a mess. In this 
thesis we will only require a small part of the calculations which give additional evolution equations for radiation other 
than those given by ( |2.32D -(|2.38|). 



For a radiation field, the energy-momentum tensor is given by 

Tf{x') = fi^u^u" + '^fi^h^' + 2q^^u'^ + - J pyfix\p^)d^p, (2.79) 

where 

pa ^ ^(^^a j^^ay^ -^p^Ua, (2.80) 

is the photon momentum, with a unit spacelike vector, and (Pp = EdEdfl is the volume element on the future null 
cone at a;', /(x^p"^) is the photon distribution function, which gives the number density of photons at a;* with momenta 
p", and must satisfy the Boltzmann equation, 

f=C(/), (2.81) 

with C being change of / along the geodesies (parameterised by v) from all types of collisions, absorptions, etc., and 
which is effectively zero after decoupling. The terms in ( 2.79| ) define the first three multipole moments: /i^ is the 



monopole, and represents the average temperature over the whole sky via the Stefan-Boltzmann law; 

(T^) (2.82) 



sky 



so if the distribution function is a Planck distribution, then the temperature is that of a black body. Any fluctuations 
in brightness across the sky are given by the higher order multipoles, tt"'', n°i°^'"°*. Now, if the observers' frame 
(ii°) is moving non-relativistically relative to the CMB frame (m°) then these evolve according to 

Ah + + Vaq^ + 2Uaql + O a^T^f = n ^ _ ^a^^ ) ^ ^[g] (2.83) 

if'^^mnc + n^^t - n^a^ {\ii,v'^ ~ + <^b) + 0[3] (2.84) 



which are the new forms of ( |23^ ), and ( p3^ for the case when the radiation field is not quite aligned with the (baryonic) 
observer. There are evolution equations for all the higher order multipoles as well, which are not given in a simple fluid 
description: the quadrupole evolves according to 



R '3 



- -J:^'^ - §n'"'=^c - h^vK'^ ) + 0[3] (2.85) 
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while the higher muhipoles {£ > 3) evolve according to 



tliAi) 



(2£ + 3) 



(^-l)(^+l)(£ + 2)_ 



5£ 



(2^ + 3) (2^ + 5) 



n-^* -n 



(Ae-i^at) 



{2e + 3) 

£+1 ■ 



2^ + 3 



n 



Aea„ 



0[3]. 



(2.86) 



For £ = 3, the n^'^^-if"*^ term on the right-hand side of equation (2.86) must be multiplied by |. is the free electron 
number density and is the Thompson scattering cross section. The series expansion terminates at 0[3] = 0{ev'^,v^) 
where e is a small parameter to represent how close the radiation and baryonic frames are: there is no neglect of physical 
and geometric quantities. All that is assumed is that the matter moves non-relativisticaly. 

These results are required for the EGS theorem later: without them (especially without (2.85)) Theorem ^ cannot 
be proved in the 1+3 formalism. 



2.7 Observations in the 1+3 Formalism 

Observations in cosmology are made by observing light^ which has travelled on our past light cone. They are made 
essentially from one spacetime point - 'here and now'. All we can hope to know directly from this is where the light is 
coming from and how it travelled here. The point of cosmology is to infer as much information about the universe as 
possible, hopefully with as few assumptions as possible. As a generic problem, the absolute limit to the amount we can 
know has been shown by Ellis et al. (1985) to be that part of the universe enclosed in our past lightcone, even assuming 



infinite precision in the observations (and assuming Einstein's equations (2.2)). 

The part relativistic cosmology has to play is to determine the cosmological model (or indeed the entire class of 
models) which fits the observations best: that is, what is the metric gab and matter content Tab = Gab (and necessarily 
the fundamental congruence, u"") of the universe? This is distinct from other areas of cosmology insofar as it need not, 
in the first approximation, describe the details of structure, or how the structure formed. It must, however, be able 
to describe gross features like the CMB; for example, the anisotropy of the CMB has been used to limit some global 
properties of m° such as the shear - see Ellis, Treciokas and Matravers (1983a,b); Stoeger, Maartens and Ellis (1995); 
Maartens, Ellis and Stoeger (1995); and Theorem ^ 

We must be able to relate the light that is observed to the metric gab and velocity field The most important 
methods are also the most direct: light observed from discrete sources, and light from the CMB. 



2.7.1 Observable Quantities From Discrete Sources 

The foundations of this subject come from the now 'famous in the right circles' paper by Kristian and Sachs (1966), who 
first treated the subject in a covariant manner (see also MacCallum and Ellis (1970)). Their method was to identify 
the various distance measures available in cosmology and relate them to redshift using the intrinsic and observed 
brightnesses (or equivalently, magnitudes) of the sources. 



Luminosity-Distance Relations 

If we observe a star to have some flux JF that has some intrinsic luminosity C then its luminosity distance is defined to 
be 




(2.87) 



This is a simple extension of the variation of the brightness of objects in non-relativistic physics because the flux that we 
observe not only decreases because of the decreased number-density of photons, but also because each photon is losing 

^Or gravitational waves, or neutrinos. 
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energy due to the expansion of the universe[[ This loss of energy of each photon is given by its redshift, z; 

-t , ^<,f '^observed ^emitted ^^o /r^ oo\ 

1 + zM- = = aT' ^ 

^emitted ^observed ^^e 

which shows redshift to be a time dilation effect. Although this last equality is true for photons,^ it is actually true for 
any null-connected points in any (smooth) spacetime. Redshift, therefore, reflects the stretching effect of expansion of 
spacetime; since it is a directly measurable quantity, it is of fundamental importance. It is therefore of some importance 
to write down other quantities in terms of redshift. There is, however, a problem, and one which is quite difficult in 
practice. If a galaxy is observed to have a redshift z then determining how much of this is cosmological, Zc, and how 
much is due to the relative motions of source (ze) and observer (zq) requires a careful study of the galaxy movement 
and a knowledge of our relative motion. The redshifts are related by 

1 + Z= (l + Ze)(l + Zo)(l + Ze), (2.89) 

although the Zq term can be inferred from a knowledge of the CMB dipole ~ the CMB frame. 
Redshift 

Redshift can be treated in a more rigorous way. By solving Maxwell's equations on a pseudo-Riemannian manifold in a 
charge and current free region, and assuming that the wavelength of light is small compared to the spacetime curvature 
(the geometrical optics approximation), we find that light is tangent to null surfaces of constant phase (f>, and therefore 
travels on null geodesies. So if a photon's velocity is given by a null (geodesic) vector fc° M V^'cj) (<^ '^[a^b] — 0): 

= Va0V°(/) = fc^fc^; fc''Vafc'' = 0, (2.90) 

(the first of these is the eikonal equation - see Ehlers and Newman 1999) then the angle between this vector and some 
velocity field u° is the relative (angular) frequency of the photon as measured by an observer traveling on 

UJ = -U^ka. (2.91) 

This frequency is clearly observer dependent; an observer moving in a different manner would measure a different 
frequency (doppler shift). If a photon travels between two points then the redshift is the relative change of frequency 
between the two points; 

We can see this more explicitly by calculating the rate of change of frequency along the photon path; 

k'^VaiO = k^'Vai-Ubk'') = -rfc''VaUh = {-^eKb + UaUb - (Tab) k'^kK (2.93) 

k"' is a null vector, and so can be written as a linear combination of a component parallel and orthogonal to viz; 

fc" = -Wbfc^(w'' + 6"). (2.94) 
where e° is a normalised spacelike vector ort hogo nal to c^Cq = 1; Uac"" = 0; it simply defines the direction of the 



photon relative to Substituting this into ( |2.93D we find 



rVat^ = -w^ {\e + e°-Ua + e°-e^aab) ■ (2.95) 

It is now clear that the relative direction of the photon is important when considering frequency change; important 
effects also come from the kinematics of m". For example expansion will decrease the photon frequency and thus 
increase the wavelength, making it redshifted. However, the change in frequency caused by shear and acceleration 
is direction dependent: acceleration can either increase or decrease the frequency, depending on the direction of the 
incoming photon relative to the direction of acceleration; similarly shear will decrease the frequency by varying amounts 
depending on the direction of the incoming photon and the principle directions of the shear (ie, its eigenvectors). 

^For ease of discussion I will assume throughout this thesis that the universe at some time (ie, now) is expanding - although models can 
be constructed which are static that give similar observed redshifts - but the analysis applies in any relativistic model. 
*This is true in the geometrical optics approximation which basically says that photons move on null geodesies. 
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The Redshift Structure of Conformally Related Spacetimes. 

If we have two conformally related metrics gab — e^^cjab'- Q > 0, then null geodesies, gabk"'k'' = 0, are conformally 
invariant: 

fc''Vfc/c" = ^ k'^^bk" cx k"", (2.96) 

where Va and Va are the covariant derivatives associated with gab and (jab respectively. Although the null geodesies 
associated with Va are non-affincly parameterised, the structure of the lightcones is identical for both spacetimes (see, 
for example Wald 1984). The afhne parameters A, A associated with the Va-geodesics and Va-geodesics are related by 

^=<ie"'- (2.97) 
aX 

with q constant. Now we associate a different null vector with each metric: fc° is tangent vector to a geodesic with affine 
parameter A, and k"^ is tangent to a geodesic with afhne parameter A. In some basis, 

since the basis vectors are linearly independent, we now have established the general result 

= qe^^k". (2.99) 

Hence, if the geodesies with respect to Vq are known explicitly then we can find them easily for Vq. 

Suppose we have a metric gab in comoving coordinates, = — |3oo|~^^^^o ' which is conformally related to another 
metric gab and there exists a coordinate transformation x" {x") : gap — > Oa'p' where ga'p' = exp[2Q{x'^ )]5q'/3' 
explicitly in the primed coordinate system. This means that the frequency of a photon travelling along a null geodesic 
becomes 

uak'^ ^ uap^k-' = -\g,,\-^/Hg^k'^' = -^2^L^^!^^fc"' (2.100) 
It is then straight forward to show that if ga' p' is in comoving coordinates the redshift becomes 

where z is the redshift associated with gab- 
Area and Luminosity Distance 

In addition to ( ^.87 ) as a distance measure, we can also measure the angular size of an object, dfl. If this object has 
an intrinsic proper area dS, then the area distance is defined by the ratio of these; 



This is related to the luminosity distance ( 2.87 ) by the reciprocity theorem, which states 

ri-r^(l + z)2, (2.103) 
as first proved by Etherington (1933). It is essentially a geometrical result, but can also be viewed as a time dilation 



effect, relating geodesies traveling up and down the null cone. See figure 2.1 



Area distance is measurable, and can be found by integrating the geodesic deviation equation - see MacCallum and 
Ellis (1970). Using ( plo^ ), and ( p^ we can relate the metric to directly measurable quantities: redshift, luminosity 
and area distance]^ 

T = ^ (2-104) 

47rr^(l + z)4 ^ ^ 



'Relating area distance to the metric and redshift is not simple in general, and results are model dependent. 
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dS„ is the intrinsic 

(jr 



area of the galaxj 



(Kl^is the angle subtended at the observer 




Sphere of coordinate 
radius r at cosmic 
time T 



Galax 
at time T 



Spatial section at time T 



Worldline of 
Observer, r=0 



Figure 2.1: An observer measuring the area distance of a galaxy. 



The Magnitude-Redshift Relation 



If we take the logarithm of (2.104) we get the magnitude-redshift relation; 

m — Af — 25 = 5 logjQ r^. 



(2.105) 



As it stands though, ( 2.105| ) is virtually useless until a cosmological model is given (or, at the very least, a metric). 
However, the cunning Kristian and Sachs (1966) managed to expand va in a power series in z, and thus making a 
generalised magnitude-redshift relation. The result of that expansion is, in the notation of MacCallum and Ellis (1970), 



M - 25 = 5 logio z-5 logio K^K^VaUb I 



logi 
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where 



u''kb\(j 



-K + e")lo. 



(2.106) 



(2.107) 



Obviously, (2.106) is extremely complicated for a general cosmological model - even for the lowest order terms. It 
is also questionable how accurate this series expansion will be when truncated at low order in redshift. Certainly it will 
be useless for z > 1 unless the exact function m(z) has a very special form, or all terms > 0{z) die off very rapidly - 



see figures 2.5 and 2.6 



Number Count-Distance Relations 

A potentially interesting area of observational cosmology comes from number counts as a function of either redshift 
or magnitude. It works by simply counting galaxies at a certain redshift in a solid angle of the sky dfl up to some 
limiting magnitude. If we know the (mean) density of galaxies in the spacetime then by counting how many we see in 
this volume will give us information about the spacetime geometry. 

We do not use number-counts in this thesis (although the relevant formulae may be derived) because there is large 
uncertainty in the source evolution function (Mustapha, Hellaby and Ellis, 1998 - although it is unclear if their result 
holds if multi-colour observations are taken into accountPl). 



2.8 FLRW Models 

The standard model of modern cosmology is based on a hypothesis: the Cosmological Principle (CP; Ellis 1975). One 
of the key issues in this thesis is the validity of the CP. The cosmological principle assumes perfect isotropy about 
our location and extrapolates this to every other location using the Copernican principle which necessarily implies 
homogeneity of the universe]^ Once the CP is in place then one must conclude that the universe has an FLRW form. 
The proof of this is intuitive, and can be found eg, in Wald (1984). 
The metric of FLRW models has the form (in isotropic coordinates) 

ds^ = -dt^ + , "^.^f, Ut' + r'dV?) (2.108) 
(1 -I- \kr^Y ^ ' 

where 

dVi^ ^ dd^ + sin^ ^d(/- (2.109) 
fe is a constant which characterises the spatial hypersurfaces of the model: k <, =, >, corresponds to hyperbolic, flat. 



or spherical geometryj^. Almost all the properties of the FLRW models follow from (2.108). The fundamental velocity 
field is given by 

w° = 1; = 0, (2.110) 

which has expansion 

m-i% (2.111) 

a{t) 

with the shear, rotation, and acceleration vanishing. The energy-momentum tensor automatically has the form of 
a perfect fluid spacetime with respect to this congruence. In fact, the FLRW models can be invariantly classed as 
perfect fluid spacetimes with a fundamental congruence that has Uab = ^ab = u"" = 0, as can be proven from the field 



equations (2.32-2.43) - see Krasihski (1997). 



^"B. Bassett, private communication. 

^^Any universe which is isotropic about three points will be homogeneous. 

^^These are often called open, fiat, or closed geometries on the assumption that only the spherically symmetric case has closed spatial 
sections. This is not true in general and an open model may not be infinite - a non-trivial topology may be imposed to give eg, an open 
model with closed spatial sections (eg, Luminet and Roukema 1999; Cornish and Spergel, 1999, show that such a model is favoured over an 
infinite open model, on the basis of COBE data. 
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The energy density and pressure are given by 



k 




3^ H 










k 


-2^ 




a 





A, 



■A, 



(2.112) 
(2.113) 



as can be foun d usi ng ( 2.24 ) and ( 2.25 ). Eq. ( ^.112 ) is the Friedmann equation. Combining these gives the Raychaudhiui 
equation (cf, ( 2.33 )) 

-i(/i + 3p) + iA. 



a 
a 



(2.114) 



The Friedmann equation ( 2.112 ) and ( 2.114 ) gives equation (|2.32|) . 

The function a(t) is free, but will have a specific form once the thermodynamics is settled upon; this usually takes the 
form of a barotropic equation of state, or multiple (non-)interacting fluids. As far as the standard model is concerned, 
the matter is assumed to be radiation until decoupling and dust thereafter. Once the matter specified ( 2.112 ) and ( 2.114 ) 
give a differential equation for a{t), which can be solved in principle - see eg, Ellis (1998). 

We can define various scalars to get an idea about the dynamics of any particular model. We define the Hubble 
scalar as the expansion rate 

£ a 



the deceleration parameter 



and the density parameter 



aa 

• 9 ' 



3ff2 



A 
3ff2 



1. 



(2.115) 
(2.116) 
(2.117) 



where the last equality (from ( 2.112| )) shows that the density parameter may be used to characterise the spatial surfaces 
in a A = model. We can also define 

^^A ^ (2.118) 



and 



k 



which implies that the deceleration parameter becomes, using (2.114), 

g = ijiZ 1 i + - I - r^A; 
for dust, we have the simpler, more familiar, form 



From (2.117) we have the normalisation condition 



n + nA + ^k = 1, 



(2.119) 



(2.120) 



(2.121) 



(2.122) 



which means, even though all the parameters are functions of time above, that we can always characterise the curvature 
of the spatial surfaces in the following way: 





h > 1 


closed spatial surfaces 


k>Q 




< 




h^A = 1 


flat spatial surfaces 


k = 


nk 


= 




h ^A < 1 


open spatial surfaces 


k<0 




> 
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2.8.1 The Magnitude-Redshift Relation 



The magnitude-redshift relation is quite easy to derive due to the symmetry of the solution, as can be found in, eg, 
Peebles (1993). The key steps are outlined below. (Factors of c, the speed of light, are included in this section.) We 
are aiming to find the funct ion m (z) - a relationship betwee n mea surable quantities. The informatio n we have about 
the spacetime is the metric ( ^.108 ), the Friedmann equation ( ^.112| ), and the Raychaudhuri equation ( 2.114 ). What we 
don't know is the function a{t) which is going to be needed if we want to compare the models to data. The usual route 
is to assume that the universe is dust after decoupling, and solve ( 2.113 ) for a{t). 

The magnitude-redshift relation in the case of A = was first solved by Mattig (1958). The general formula for 
vanishing pressure was first found by Kaufman (1971) after the initial generalization by Kaufman and Schucking (1971) 
to include A > in a closed model. 

We will derive it explicitly here to facilitate comparison with later derivations for the Stephani models. 



Redshift: For a general spherically symmetric spacetime it is necessary to solve the geodesic equation for the null ray 
connecting the galaxy to the observer in order to determine the galaxy's redshift. More specifically, for a null ray from 
the galaxy, G, at (r, r) to the observer, O, at (0, tq), the tangent vector along the ray fc" is obtained as a solution of the 
geodesic equation. Once we have /c° along the ray, though, the redshift is obtained immediately from (cf, Ellis 1998) 



1 



(2.123) 



where is the four- velocity of the perfect fluid, ie, of the galaxy (G) or the observer (O). Since u° is a four- velocity it 
is normalised by 

gabu'^u'' - -1, (2.124) 

and in comoving coordinates only m° is nonzero, so that ( 2.124| ) completely fixes However, the high symmetry of 
FLRW models means that the redshift can be obtained directly without having to integrate the geodesic equation - see 
eg, Wald (1984) for a derivation. However we can use a novel approach outlined in j |ll| , which can be implemented if 
we write the metric in the conformally static form 



ds^ = ait{ri)y [dr]^ + (1 + ^kr^)-^{dr^ + r^dfl)] 
where rj is the conformal time coordinate; 

f dt 



(2.125) 



(2.126) 



Using (2.101), we immediately have (because the part in square brackets in ( 2.125| ) is static and there are no gravitational 
redshifts) 

1 + ^ = ^; (2.127) 
a{t) 

where qq = a{T) is the scale factor today. 



Lookback time: Consider radial null ray (ie, dd — dip ^ ds^ = 0) from a source at (r, r) reaching the observer (r — 0) 
at time tq. The relationship between tq, t and r can be obtained directly from the metric symmetry guaranteeing that 
the raypath is purely radial so it is parameterised by a function T(r) by integrating 



dr 



dt 



(1 + ifcr2) a{t) 



between r and r = 0. Now we have, on using (2.127) 



dt 



da 1 



dz 



where H is given by (12.1151) . From (^J0|) we have 



a{tf 
H = fio — 3- 



(2.128) 



(2.129) 



(2.130) 
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so that (2.112) becomes, using (2.127) once again, 

H{zf = Hl [rio(H-z)3 + r!fe„(l + z)2 + r!Ao] ; 
where a subscript '0' means the present day value; viz; 



^0 - 



k 



3Hi 



A 



Together with ( ^.12^ ) this gives the function r{z) 



fc > 
fc = 
fc < 



^tan-i(iVfcr 



^tanh"^ f iVlfclr) 



dz' 



(2.131) 
(2.132) 



(2.133) 



We now have the coordinate distance as a function of redshift. We must now relate this to some measurable distance 
quantity. 



Angular diameter distance: Angular size distance, ta, is the ratio of an objects physical diameter to its apparent 
angular diameter. As is well known, in spherically symmetric spacetimes the angular size distance of an object as seen 
from the centre is given by (the square root of) the coefficient in front of the angular components of the metric (dil^) 
evaluated at the time the light was emitted. 

For any inward radial null rays, in the plane 9 — 7r/2, say, the angle coordinate cj) obviously does not change along 
the path. This means that the apparent diameter, of a galaxy is just determined by the difference between the 
coordinates at either edge of the galaxy, at the time the light we observe was emitted (ie, at t{r)): Sq — 5(f>. The physical 
diameter, D of the galaxy is given from the metric (with dt — dr — dd = 0): 



D = 



a{t{r))r 
1 + U"r2' 



a{t{r))r 
l + ifcr2' 



{a{t{r))r/{l + jkr'^) just being the coefficient in the angular part of the metric). The angular diameter distance, ta, is 
simply defined to satisfy D — taSg, so that 



a{t)r 



ao 



r{z) 



1 



\kr^ 



l + ifcr(z)2 1 + z 



(2.134) 



for FLRW models, where we have used ( 2.12g ) in the second equality. 

The area distance can behave in quite an odd way: in non-relativistic situations when things move further away 
they get smaller. This is not necessarily what happens in a cosmological model: as something moves away it may get 
smaller at first but then start to get bigger again. This is shown in Fig 2.2 where we can see that in a closed model 
this refocusing occurs at some redshift. This is because the light rays are effectively traversing an expanding sphere, 
so a galaxy situated at the antipode to a particular observer will be spread all over the observers sky. In (dust) FLRW 
models, this refocusing only occurs once as z — s- oo because the expansion is always just fast enough that the light rays 
never 'catch it up.' 



Luminosity distance: The key to deriving the luminosity distance, r^, is the reciprocity theorem (see Ellis and 
MacCallum, 1970, and references therein), which allows us to write 

rL = {l + zfrA- (2.135) 
This equation is exact, and completely general, applying in any spacetime. Combining this with ( ^.134 ) gives 

(1 + z)aor{z) 



rhiz) = 



1 + ifcr(z)2 



for FLRW models. Note that we can write this in the more familiar form 

f(z) — -j^ sin 2 tan^-'- \\fkr(2 



rL{z) = (1 + z)aof(z) 



r{z) — r{z) 



f(z) = ^= sinh 2 tanh ^ ^J\k\r(z) 



fc > 0, 
fc = 0, 
fc < 0. 



(2.136) 



(2.137) 
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IN? 




Figure 2.2: Area distance as a function of redsiiift for a variety of 'standard' models. Tlie closed model has {flo, f2Ao}={0.3, 1.5}, 
Einstein-de Sitter (1932) has {1, 0}, flat: {0.3, 0.7} and the open model has {0.1, 0}. All the models have Ho = 60 km s"^ Mpc"^ 
Note that all the curves match up for low enough redshift. 



The function r{z) is plotted in Fig. |2.3| . Hence we have the result 

1 



rL{z) = (1 + z)ao 




dz' 



j """^^ \/^o(l + Zf + (1 + z)2 + Oa„ 



fc > 
fc = 
A: < 0. 



(2.138) 



Magnitude-redshift relation: The magnitude-redshift relation for an object of absolute magnitude M then follows 

(2.139) 



from (2.13g) and 



m(z) — M — 25 = 51og;^Q (1 + z)ao x 



1 



1 



smr 



dz' 



H^A^\ J V^o(l + zf + f7fc„(l + Z)2 + fiA„ 



fc > 
fc = 
fc < 0. 



In Fig 2.4 is shown the m ~ z relation for the same models as Fig. 2.2. The closed model has the refocusing issue as 
before; the apparent magnitude becomes singular at that point. 

As the integral in ( 2.140 ) is elliptic, it cannot be evaluated explicitly. However, a series relation can be found 
relatively easily; 



m — M — 25 = 5 log;^o TT 
+ ^logioe|(l-go)z + 



-(3go + l)(9o-l)-2r!Ao 



z^ + 0{z^) 



(2.140) 



We can now consider quantitatively the comments made in § |l2| on the validity of the low o rder e xpansion. Consider 
Fig 2_^: the plot shows the difference between the exact m{z) function and the series form ( 2.140 ) with an increasing 
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-o.i-L 




Figure 2.5: The exact magnitude-redshift relation minus tlie same function as a power series relation to various orders. \m^ 
'^series I > 0-1 would give a significant detectable error. 



power of redshift z up to order 9. What we can see from figure 2^ is that all the series expansions diverge from the true 
function to more than O.lmag (which is roughly the errors in the magnitudes of the SNIa) at some redshift less than 
unity. While not massively important for the FLRW models at low redshift, for supernovae at high redshift z ~ 1 it 
will be increasingly important to use the exact relation, as even the ninth order expansion diverges from the true result 
at these redshift, by a measurable amount. 

If we jump ahead somewhat and consider the issue for a non-FLRW model (the Stephani model considered in 
substantial detail later) then we can see from Fig (2.6) that the situation is far worse: the series expansion can diverge 
very dramatically even at low redshift. 

We shall come back to the series form of the magnitude-redshift relation later in relation to general cosmological 
models. 



2.8.2 Observable Quantities in FLRW Models 

For a dust FLRW model there are only 3 independent parameters to fit from data: {HQ,^Q,^\g}. In practice, these 
have boiled down to the equivalent set {i?Oj90i^o}- This arises from the fact that most measurements are made at 
z < 1, or z ~ 1000 at the CMB. Although the measurement of all three parameters may be made in principle from 
any observations, it is essentially impossible to determine any more than two of these from observations in a narrow 
range of redshift. This has been demonstrated explicitly in a series of papers by Tegmark, Eisenstein and Hu (Tegmark, 
Eisenetein and Hu, 1998a; Tegmark et ai, 1998b; Eisenstein, Hu and Tegmark, 1998a, b; see also Tegmark, 1999; White, 
1998; Efstathiou et ai, 1998; see Ehlers and Rindler, 1989, who first introduced the 'phase plane' arguments used here). 
Essentially they have shown that the set of possible rio,f^Ao found from low-redshift observations will complement the 
set of values preferred by high redshift experiments (ie, CMB results). (It should be noted that they stress the need to 
fit all the data sets at the same time, instead of taking their conclusions separately.) 

For the low-redshift experiments (z S 1), we can see from the series expansion of m(z) ( 2.140| ) that they will be able 



to determine go to some accuracy, but the extra detail required for a measurement of and f2Ao separately would need 
knowledge of the 0{z'^) term. For example, the latest supernovae la results of Perlmutter et al. (Perlmutterei aL,1999) 



give the results of their data as an exclusion plot in the flo — flAo plane. Fig. 2.7, We can see in this figure that their 
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best fit regions do not provide any substantial constraint on D,q and Q\g , but they do provide quite good constraints in 

qo- 

At high redshift a very similar thing happens, but in a less obvious way. The series expansion of the rn(z) curve 
is clearly useless for any intuitive idea of what's going on. However, if we consider the CMB to have an intrinsic and 
apparent magnitude, then we can examine the exact m — z relation at z = 1000 and plot the contours of constant m 
on the r^o — ^^Aq phase plane. See Fig. 2.S, where m — M — 25 is plotted in the fig — fl^g plane, with contours of 
constant m — M — 25 marked on it. In practice we can calculate m — M — 25 for any object from observations (giving 
m) and assumptions about the objects intrinsic brightness (M). Thus we can select just one contour in (2.8), and get 
some limitations on fio and SIa^ . Just as in the SNIa results, the constraints will be weak, with any particular contour 
selecting a large range of allowed values for and (unless fio is small). It should be noted that this is a very 
simple analysis, because the degeneracy of the CMB in the Og — plane comes from the position of the first Doppler 
peak (cf, figure 1.2) rather than simply the m — z relation, but the principle is the same. 

Tegmark et al. have noticed that when the results from the CMB are taken together with the SNIa results, then it is 
possible to get very precise predictions of llo and separately. If we have a look at Fig. 2.8 we can see why. We have 
plotted contours, as described above, for redshifts of 0.5, 5, and 1000. Typically, a set of observations at low redshift 
will determine one of the blue contours (obviously this means a set of contours in a small neighborhood of some mean or 
'best fit' contour); similarly, the CMB measurements will determine a similar set of the yellow contours. Because these 
contours intersect at a fairly large angle, it is only the overlap of these curves where the real universe will lie (provided 
the universe is FLRW of course) . 

For completeness, we show in Fig, 2.9 the results of Tegmark, Eisenstein and Hu (1998a) who computed this overlap 
for future SNIa and CMB observations. 

We also show another stolen picture. Fig, 2.10 , which gives the same sort of results as Fig, 2.8, but for number-counts, 
age and growth as well. 

What is also interesting about figure 2.8 is something that none of the papers above have mentioned; a quick look 
at the figure shows the contours chang angle from ^ 45° to > 90°. This means that there will be some redshift {z ~ 3) 
for which the lines will be vertical; ie, JIq will be determined to very high accuracy from just one set of experiments 
independent of JIa- In fact it may be worth specifically looking for eg, quasars or SNIa at this redshift to give a very 
accurate determination of f2o (and hence f^Ao from other redshift surveys). 
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Supernova Cosmology Project 
Perlmutter et al. (1998) 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




12 3 



Ap.J. 

astro-ph/9812133 



Figure 2.7: The best fit to tiie m{z) function for tlie supernovae data of Perlmutter et al, (1998b). the confidence regions do 
nothing for determining Qo and Qaq- 
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Figure 2.8: The phase diagram for the m — z relation for z — 0.5, 5, 1000 presented in the Qo — ^Aq plane. We can see that the 
contours move round the diagram with increasing redshift. The faint dashed lines are constant qo, which shows that low redshift 
observations will give good results for go- 
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Tegmark, Eisenstein, Hu & Kron 1998 
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Figure 2.9: A futureistic diagram of what might be expected from SNIa and CMB complementing each other in the future. The 
errors on f2o and Q,Aa (small black ellipse) are extremely small (by cosmology standards) in comparison with the independent 
results. Plagiarised from Tegmark, Eisenstein and Hu (1998a). 



2.9 Observable Quantities in a General Spacetime 

In the FLRW models, there is a simple relation between the fundamental length scale a and the observable parameters 
HqjQq etc. It is interesting to determine the correct covariant definitions of Hq and qo for a general cosmological 
model and to ask whether the FLRW SIa generalises to other cosmological models. Almost invariably these parameters 
are measured fitting distance-redshift data to an FLRW model (essentially eq. (2.4)). What can this say about other 
inhomogeneous models; and what could be causing a negative qq other than a cosmological constant? 



The standard generalisations of Hq and go are based on the fundamental length scale defined by (2.17) simply by 
generalising the definitions for FLRW models, (2.115) and ( p. 116 ), viz; 



H 



I 3 



Using the Raychaudhuri equation ( ^.33| ) we find simply that 



(2.141) 
(2.142) 

(2.143) 



(See eq. (52) in Ellis, 1998. Note that his statement about the acceleration terms being small from CMB anisotropics 
is false - see Chapter |[) 



The definition of the expansion Q is the volume change along defining H via ( 2.141 ) is misleading. It gives the 
average distance change to neighbouring particles along u°; whereas the actual distance change to a particle in direction 
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Tegmark, Eisenstein, Hu & Kron 1998 




e° is given by 



m 

Si 



^^9 + aabe-'e' = H- 



(2.144) 



as one would expect, the shear will affect the distance to neighbouring particles, as some may be moving towards or 
away from you relative to the overall expansion - see Ehlers (1993). Now this is a good definition of the Hubble scalar, 
and has been used before. Often, the expansion tensor is defined as (eg, Wainwright and Ellis 1997) 



Qab = (Jab + \Ohab 



(2.145) 



so we can now define 



(2.146) 



to give a covariant, direction-dependent definition of the Hubble scalar. This has been used in Humphreys et al. (1997). 

If we want to compare a non-standard cosmological model with the real universe (or with an FLRW model, with 
values of Hq etc. derived from them) then we need a generalisation of Hq and qo suitable for a decent comparison with 
data. It has been suggested that equating the various order terms in the generalised m — z series relation with those 
of the FLRW series would give the proper generalised and (Humphreys et al. 1997; Paczyhski and Piran 1990; 
Dominguez-Tenreiro 1981a,b; also see Palle 1999). So, comparing (|l0§) with ( |2ll0| ) we can define 



H 



m K'^K^^aUbL, 



^ K''K''K^^aVbUc 
~ {K'^K'^VdUef 



(2.147) 
(2.148) 



Expanding H gives 



a A) . 



+ Mae" -I- (Tabe"■e^ 



(2.149) 



monopole ^ipole quadrupole 
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where the quantities on the right hand side are understood to be evaluated at the present time. This definition does 
not have the physical significance of ( 2.146| ); rather it is the zero point of the linear Hubble relation. To compare a 
non-standard model to an FLRW model at low redshift, then (2.149) is the definition which must be used. Eq, (2.149) 
shows that a general inhomogeneous cosmological model has the effect of inducing a dipole and quadrupole in the m — z 
relation as z ^ 0. However, this conclusion only concerns measurements made at very small redshift, say z ;S 0.1, 
because, as we have seen in § ^.8.1 , we cannot be sure that the low-order terms in the series expansion will give a 
function even approximating the real thing for larger redshift (say, z Z 0.1). In other words, at higher redshifts (eg, for 
SNIa) higher order terms may modify these moments to some degree. 

All of these low redshift surveys essentially use the linear Hubble Law which, for a general cosmological model reads 



z = K'^K'^aUblvA + 0{r\) = <Va + 0{r\). 
If we write this in the more suggestive form 



(2.150) 



Hnr 1 



Hn 



(2.151) 



where Hq is the expansion rate ( ^.14l| ) as in the FLRW models (basically the all sky average expansion rate) Equation 
(2.151) reflect the effects of local small-scale inhomogeneities, ie, perturbations of the smooth background model. If 
we and add in these local effects, (this is often done in terms of spherical harmonics, eg, Regos and Szalay 1989; or 
Branchini et al. (1999); Tadros et al. 1999 - but we shall just use cost? expansion) represented by constants 6i, 62, 63 • • •, 
we have 



i/oH 1 



Ha 



Ho 



hi cos -d + 62?" cos ■§ + h^r cos d -\- 



(2.152) 



where ta has been replaced by r because all definitions of r agree at very small redshift. hi corresponds to the bulk flow 
of the local group moving relative to the CMB frame, and 63 often corresponds to a shearing effect from our infall to the 
Great Attractor and Virgo cluster - see eg, Lilje, Yahil and Jones (1986) who first reported the effect. However, there 
is no physical local effect that will produce 62 7^ 0. It is certainly not clear at all that 62 = on observational grounds; 
this leaves acceleration as the only effect that w ill produce a dipole in the Hubble law which grows linearly with distance. 

If we try to determine from equation ( 2.152 ) from galaxy surveys then we must be able to determine how much of the 
observed dipole is bulk flow and how much may be attributed to any possible acceleration term; ie, whether any of the 
dipole grows with distance. In principle we may be able to limit the acceleration quite strongly in this way. However, 
preliminary results using the IRAS catalogue only constrain 



half 5ky 



Ho 



< 0.2. 



(2.153) 



(Clarkson, Rauzy and Barrett, in preparation.) This is not strict at all, but it is likely that a more complete survey 
such as the PSCz catalogue may provide stronger constraints. 

T he ge neralised Qq is not nearly so simple; expanding ( 2.148| ), and substituting for 9, a-(ab) and Vbcr'^'' from ( 2.33 ), ( 2.36| ) 
and ( 2.40| ) respectively we find 



obs 

Qo 



Ho 

+e'' 



6^2 2, ,2 

5^^ ~ 3^ 



"^^{aUb) + U{aUb) + E^ab) ~ {ab) + f CTc(a<^b) 



-e e 



+e''e°e'' 



3 c(a b) 



-^U{a'^bc) - '^{aO'bc) 



-^{a^b) - 

3e''e''e=e''cr(afcCTc<j) } • 



monopole 

dipole 

quadrupole 
higher multipoles 



(2.154) 



This is somewhat complicated for an observable quantity. Measurements of go are nowhere near accurate enough to 
detect any multipole moments; in fact it is only recently that the sign of the deceleration parameter has been decided 
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upon, let alone a variation in the sky. However, if we consider the mean value of (2.154) around the sky, then we hav( 



m 



ohs 2 ob5 

H q 



a , 6 „2 2, ,2 



(2.155) 



which, we can compare to the definition given by the length scale ( 2.143 ): 



(2.156) 



We find that the shear and rotation become less important in the observationally derived relation, while the acceleration 
terms become relatively more important. More interestingly though is the contribution from the acceleration vector u^, 
has disappeared: it will actually contribute nothing to the value of the deceleration. 



^''Obviously the multipoles of odd order are zero when averaged over the sky, and the even order multipoles have an average value of zero 
because they are contractions of traceless tensors. 



Chapter 3 



The Ehlers-Geren-Sachs Theorem and 
Some Generalisations 



The high isotropy of the CMB is usuaUy taken as strong evidence that the universe is homogeneous and isotropic, ie, is 
well described by an FLRW model. The principle justification for this is an important theorem of Ehlcrs, Geren and 
Sachs (1968) (based on earlier work by Tauber and Weinberg 1961), which states that if all observers in an expanding, 
dust universe measure an isotropic CMB then the universe is FLRW and the cosmological principle is valid. The 
importance of this theorem lies in the fact that it permits the homogeneity and isotropy of the universe to be deduced 
not from measurements of the actual isotropy of the universe about us, but from only measurements of the CMB, 
combined with the Copernican principle (that is, the assumption that all observers in the universe see the same degree 
of isotropy). The Copernican principle is often regarded as a powerful but untestable assumption in cosmology, although 
there are suggestions that it may be testable using the Sunyaev-Zeldovich effect, for example (Goodman 1995). Here 
we simply assume that the Copernican principle is valid and study the consequences of applying it to the observed high 
degree of isotropy of the CMB. That is, we examine spacetimes with an isotropic CMB for all observers. 

The EGS theorem has been generalised by Treciokas and Ellis (1971) to include an isotropic collision term. Ferrando, 
Morales, and Portilla (1992) find the general form of the energy-momentum tensor and Einstein's equations for spacetimes 
with an isotropic radiation field, and consider some special cases with anisotropic pressure. It has also been shown by 
Stoeger, Maartens and Ellis (1995) that the EGS theorem almost holds when applied to an almost isotropic radiation 
field. 

There are counterexamples to the spirit of the EGS theorem (that is, when some of the assumptions are relaxed the 
result fails to hold). In particular, Ellis, Maartens and Nel (1978) show that the result does not hold if the expansion 
is zero (which is obviously not relevant to cosmology), and Ferrando et al. (1992) emphasise that homogeneity does not 
follow if there is anisotropic pressure in the energy-momentum tensor. Nilsson et al. (1999) provide a counterexample 
to the almost EGS result when the Weyl curvature is not negligible. 

The basis of the EGS theorem is the Liouville equation for photons, which tells us that if a radiation field (ie, a 
solution of the Liouville equation) exists such that for every observer on some timelike congruence the radiation field is 
isotropic, then that congruence is (parallel to) a conformal Killing vector (CKV). This may be expressed more formally 
as follows (Ehlers et al. 1968; Ferrando et al. 1992): 

Theorem 1 A spacetime will admit an isotropic radiation field if and only if it is conformal to a stationary spacetime, 
which happens if and only if there is a velocity field satisfying 

(Jab - 0, (3.1) 

V[„ {u^ - \eu^) = 0, (3.2) 
where Uab, and 9 are the shear, acceleration and expansion of u"" , respectively. 

(Then m° is the velocity field relative to which the radiation is isotropic, and is parallel to the CKV.) 

The proof of Theorem]^ is very straightforward in the 1-1-3 formalism, and follows from the results of §2.6 and §2.3 
in Chapter ^ 
Proof 
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With reference to §2.6 a radiation field is isotropic (with respect to u°) iff = n'^ = — V £ > 2, and 
from ( 2.83| )-( 2!85| ) we find that this condition imphes 



4 

V„Ma + Va^„ = 

Oah = 0. (3.3) 



Now, if we write Q M — | In /i„ in the first two conditions of (|3.3|) we find that 

= 3Q, Ua = ^aQ: ^ V[, (ufe] - i^Ufc]) = 0; (3.4) 

which, together with ( |2.46| ) and Oab — 0, tells us that the spacetime is conformally stationary 9 = ii = dab = 0. □ 
In the absence of some statement about the matter content of a spacetime, or further assumptions about the 
congruence Theorem ^is all that can be said. In a cosmological context the simplest, and most common, assumption 
is that the matter is dust (implying that u°' is geodesic), which leads to: 

Theorem 2 (Ehlers, Geren, and Sachs) // the fundamental observers in a dust spacetime see an isotropic radiation 
field, then the spacetime is locally FLRW. 

Proof : 

^ This is a straight forward extension of the proof of Theorem |l|: we have for a dust congruence u° = which implies 
VaQ = 0. Now, by ( 2.22| ) we find for 9^0 that Qujab — 0: hence the spacetime is a perfect fluid with geodesic. 



irrotational, shear- free flow, and is therefore an FLRW model (see Krasihski 1997). □ 
Alternatively, we can simply assume that is geodesic. The existence of an isotropic radiation field then ensures 
(for non-zero expansion) that the energy flux relative to is zero. If the anisotropic stress tensor is zero at any instant 
(so that the energy-momentum tensor has perfect fluid form) then it will remain zero and the spacetime will be FLRW 
(Ferrando et al. 1992, Corollary 1; but note that their statement that the anisotropic stress is invariant along m° in 
general is misleading - from Eqs. (31) and (40) of Ellis 1998 we have 7r(^ab) oc 9TTab)- 

It is worth emphasising that in applications of the above results to cosmology the motion of the fundamental observers 
must be identified with the congruence For example, in § [3.1.2 all Stephani models are conformally flat, and therefore 



conformally stationary, but for most of these spacetimes the fluid congruence is not aligned with the timelike CKV. 

The matter content of the universe is not precisely known. Certainly, there is a large number of possible contributors, 
including hot and cold dark matter (in their various manifestations), electromagnetic fields etc., as well as the more 
obvious radiation and baryonic matter. In particular, the type la supernova results of Perlmutter et al. (1999) suggest 
that an important component may be a 'quintessential' scalar field. However, the forms of matter that are thought 
to contribute significantly to the energy-momentum tensor may be treated in general as perfect fluids. That is, their 
energy-momentum tensor may be written in the form 

Tab = HUaUb + phab, (3.5) 

where fj, and p are the energy density and pressure, is the timelike velocity congruence of the fluid, and hab is the 
spatial projection tensor associated with Scalar fields may also be written in this form, with m° parallel to the 
gradient of the scalar field, provided that this is timelike (see i |3.2| ). Note that if several such components are present 
there is no reason why their fundamental congruences (the u°'s) should be parallel. If they are, then they may be treated 
as effectively a single perfect fluid with the energy densities and pressures added together. If not, the decomposition 
of the energy-momentum with respect to u° for one such fluid (as in equation (14) of Ellis 1998) will contain energy 



flux and anisotropic stress terms from the other fluids (again, see §3.2). The fundamental observers will be associated 
with one such congruence. Usually the fundamental observers in standard models of the universe are associated with a 
dust-like {p = 0) component, with the result that the acceleration m° of the fundamental congruence is zero. However, 
we wish to study the consequences of relaxing this assumption and consider models with acceleration. This acceleration 
must be caused by some non-gravitational force (typically pressure gradients for perfect fluid spacetimes, but in principle 
it could be the result of a coupling between the fluid and some other component such as the electromagnetic field). 

With this in mind, in this paper we consider perhaps the two simplest generalisations of the dust hypothesis. Firstly 
we imagine that the dominant form of matter is a single irrotational perfect fluid. We do not specify what form of matter 
this corresponds to, but we allow pressure gradients that give rise to acceleration. Secondly we consider cosmological 
models in which more than one matter component makes a significant contribution to the energy density and dynamics 
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of the universe. Specifically we consider 'QCDM' models, containing a non-interacting mixture of radiation, dust 
(CDM) and a scalar field (or a cosmological constant). The observers are associated with the CDM component, and are 
therefore geodesic. The difference between this and other theorems assuming geodesic observers is that the the scalar 
field component can introduce effective energy flux and anisotropic stresses relative to the dust congruence, and so the 
matter need not behave as a perfect fluid. 

In the following section we find all irrotational perfect fluid solutions admitting an isotropic radiation field for the 
fundamental observers, showing that they form a subclass of the Stephani spacetimes and are FLRW if and only if the 
acceleration vanishes. Then in |3.2| we examine QCDM models and prove that such models must be homogeneous and 



isotropic if they admit an isotropic radiation field. Finally, in §3.3 we emphasise the importance of acceleration for these 
results and show that the acceleration of the fundamental congruence is, in principle, detectable, and measureable, in 
galaxy surveys. Two appendices contain results relating to §3.1.2, 



3.1 The Irrotational Perfect Fluid Solutions. 



We wish to consider the constraints imposed by the existence of an isotropic radiation field for the fundamental observers 
on perfect fiuid spacetimes in which the rotation of the fundamental congruence is zero. Since it follows from Theorem |l] 
that the shear of the fundamental congruence must also be zero we immediately know that all the acceptable solutions 
are members of the Stephani-Barnes family, which is the family of all shear-free, irrotational, expanding (or contracting) 
perfect fiuids (see Krasihski 1989, 1997). It only remains, then, to impose condition (3.2) of Theorem]^ and thus find 
the sub-class of Stephani-Barnes models which admit an isotropic radiation field for all fundamental observers. 

The Stephani-Barnes family contains the Barnes solutions, which are of Petrov type D, and the Stephani models, 
which are conformally fiat, although these two classes overlap where the Barnes solutions degenerate to type O (these 
solutions then become Stephani models with symmetry). The FLRW models are a subcase of these solutions. The 
Barnes spacetimes all possess symmetry (see below), whereas the Stephani spacetimes, in general, do not. In all cases 
the metric in comoving coordinates can be written in the form (we use the same notation as Krasihski 1997): 



dy^+dz^] 



(3.6) 



where F — F(t) and V — V{t, x, y, z). F{t) is arbitrary, but there are some restrictions on the form of V depending on 
the symmetries of the solution, and these will be discussed in due course (but the impatient reader may wish to note 
Eqs. |3.12| , |3.14| , p.l5| , and |3.17D . Expressions for the energy density and pressure can be found in Krasihski (1997). 
The fiiud velocity is given by (without loss of generality we can assume that y > 0) 



V 



FVj 



l" 0' 



with expansion 



and acceleration 



-sign(yf) — , 



(3.7) 



(3.8) 



uo = 0, 



■In 



Vt 



(3.9) 



d 

Ui = — : 

dx' V 

where i = 1, 2, 3, x'' = {x, y, z}. Note that ( |3.8| ) differs from the expression usually given for the expansion (in Krasihski 
1997, equations (4.1.4) and (4.9.6), for example) by the inclusion of the — sign(y() factor. Neglect of this factor is 
inconsistent since F enters the metric only quadratically, so the sign of the expansion cannot depend on the sign 
of F. The sign of 9 does depend on that of Vt, though: for the Friedmann subcase of the Stephani-Barnes models, 
for example, V is related to the scale factor R{t) by V = 1/R, and \FVt/V 
—3Vt/V — — 3sign(yi)/|i^|. This is important here because the constraint 



1 (see § |3.1.2D , so that = iR/R = 
contains the expansion. Note, too, 
that F is not a true degree of freedom parameterising distinct spacetimes, but rather represents a coordinate freedom, 
corresponding to different choices of the time coordinate. 



From (3.7), (B.S) and (3.9), the condition (3.2) leads to the constraint: 



92 



dx'dt 

which is satisfied if and only if the function V has the form 



hiVt =0, 



(3.10) 



V{t, X, y, z) = T{t)X{x, y, z) + Y{x, y, z). 



(3.11) 
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where T, X, and Y are arbitrary functions. This equation is the key additional constraint on the Stephani-Barnes 
solutions. 

It is worth noting that it follows from (B.IO) that the acceleration scalar is constant along the fluid flow for every 
observer, ie, ii^t = (where u'^ = Uau"'), as can be seen by calculating (m^),* from ( |3.9|) . In fact, it can be verified more 
generally that for any conformally stationary spacetime (ie, a spacetime satisfying (|j) and (p|)), even with rotation, 
the acceleration scalar evolves according to 



where V denotes the spatially projected gradient (sec Ellis 1998). It follows from equation (32) of Ellis (1998) that VbS 
whenever the rotation vanishes (for a perfect fluid), which is the case for the Stephani-Barnes models. 

Note that if V has the form (3.11) we can immediately write the metric in a manifestly conformally static form 



±1 

T/2 



where dr = T^tFdt, which shows that these models will indeed be conformally stationary, as required by Theorem |l]. 
We now discuss each subcase in turn. 



3.1.1 The Barnes Solutions. 



The Barnes solutions all have spherical, plane or hyperbolic symmetry (ie, they possess three-dimensional isometry 
groups acting on two-dimensional orbits, cf. § 3.1.2| ). The restrictions on the metric function V depend on which of these 
symmetries the spacetime possesses (see Krasihski 1997). For the solutions with spherical symmetry (the Kustaanheimo- 
Qvist solutions) or planar symmetry we introduce a new independent variable u{x,y, z) defined by u = = + + 
in the spherical case and u = z in the planar case. Then V is defined hy V = V{t, u), subject to the condition: 



5^ 



f{u)V\ 



(3.12) 



where / is an arbitrary function. Since V =^ V{t^u) we know from ( 3.11 ) that in order to admit an isotropic radiation 
field for all observers V must have the form 



V{t,u) = T{t)X{u) + Y{u). 



(3.13) 



Equation (3.12) then imposes a constraint on the functions X and Y , which will be outlined below. 

For the solutions with hyperbolic symmetry the constraint on V is very similar. This time we introduce the vari- 
able u — x/y. V can then be written 

V{t,x,y) = yW{t,u), (3.14) 

where W satisfies 

-^^f{u)W\ (3.15) 



with / once again a free function. The condition ( |3.11 ) now gives 

V{t, X, y) = T{t)X(x, y) + Y[x, y) = yW{t, u). 
Dividing by y and redefining X and Y in an obvious fashion we obtain 

W{t,u) =T{t)X{u)+Y{u), 



(3.16) 



in which X and Y are again constrained by ( 3.15| ). 

For all three symmetries of the Barnes solutions, then, the constraints on the metric function V essentially reduce 
to the second-order differential equations ( 3.12|) or (3.15). Imposing the condition ( ^.11 ) introduces the additional 
constraint ( 3.13| ) or ( 3.16| ). Substitution of (3.13) or (3.16) into (3.12) or (3.15) respectively and differentiating twice 
with respect to time, dividing by Tt each time (and recognising that Tt 7^ 0, X 7^ 0, so that 7^ in (3.6)), leads 
directly to the condition 

f(u) = 0. 

Barnes solutions with /(u) — are conformally fiat (Krasihski 1997, p. 142) and are therefore actually a subcase of the 
Stephani models. That is, proper Barn es sp acetimes can be ruled out: they do not admit an isotropic radiation field. 
It only remains to apply the condition ( 3.11 ) to the Stephani models, which we do in the next section. 
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3.1.2 The Conformally Flat Solutions 



The conformally flat sub-case is the entire class of conformally flat, expanding, perfect fluid solutions, and is the 
Stephani solution (Stephani 1967a,b). The function V is most often written in the form (see Krasihski 1997; Barrett 
and Clarkson 1999a, b): 



1 



(l + ifc(t)|x-xo(t)|^ 



(3.17) 



R(t) ' 4 

In general the five functions of time in V are free. For our purposes, however, it turns out to be more convenient to use 

V ^ a{t) + b{ty - 2c{t) ■ r, (3.18) 

as in Barnes (1998), again with five free functions (we adopt three-dimensional vector notation, so that c = (ci, C2, C3), 
for example). In fact, this form is slightly more general than ( 3.17 ) - see Barnes (1998). 

We must be able to write V in the form ( [3.11 ) for the spacetime to admit an isotropic radiation field for all 
fundamental observers. From (|T^) or ( |3.18D it is clear t hat t he functions X{x') and Y{x') in ( |3.1l| ) can be at most 
quadratic in the x*. Writing X, and Y as quadratics in ( [j.ll ) and equating in ( 3.18| ) all powers of a;% we obtain the 
following constraint equations: 



a(t) = aiT{t) + a2, 
bit) = &iT(t) + 62, 

C{t) =CiT{t)+C2, 



(3.19) 



where T{t) is a free function of time and the ai.2, bi^2 and Ci^2 are ten independent constants. Not all of ai, bi and Ci 
can be zero (in order that V,t 7^ i n (t3.6[) ). 

Equations ( 3.1£ ), along with (|3.6| ) and ( |3.1^ ), provide the complete set of irrotational perfect fluid spacetimes 
admitting an isotropic radiation field. Not all of the possible choices of parameters give rise to distinct spacetimes, 
tho ugh, and we outline in appendix ^ how coordinate transformations may be used to eliminate many of the parameters 
in ( p.lS| ), and determine when the models can be reduced to manifestly spherically symmetric form (c = 0). 



The Constraints on R, k and xq. 



To find the constraints on R, k and xq in (3.17) corresponding to (B.l£) first equate powers of in (3.17) and (3.18) to 
obtain 



Solving these equations for i?, k and xq gives 



c{t) = ^xo. 



(3.20) 
(3.21) 
(3.22) 



ab - 



1 w ^ 

-k{t) = — — - 
4 ab-U 



xo = 



(3.23) 

(3.24) 
(3.25) 



which are valid whenever ab— cp 7^ (otherwise V cannot be written in the form (3.17)). 

To impose the constraints (3.19) perform the transformations of appendix |^ so that c = cz as in A. 3) and ai 7^ 0. 
Then |cp — c^, and b and a are related by 



bi 
— ( 

ai 



ai 



7a - 



(3.26) 
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where 7 = &i/ai and S = b2 — 6102/01 are constants. Using this in (3.23) and ( ^3.24 ) leads, after some rearrangement, 
to a quadratic relationship between k and R: 



(7 + SR) 



2 r>2 



7C^i? 



0. 



In addition to this constraint relating k and R we can trivially rewrite ( 3.22 ) or ( 3.25| ) as 



4i? , 
xn = — cz- 



(3.27) 



(3.28) 



Equations (3.27) and (3.28) are the constraint equations on i?, k and Xq corresponding to the equations ( 3.19| ), or 
rather ([A.3|). If desired (3.27) can be solved to obtain 



7 + (5i? ± vF+^^P + 4^yc2R2 



From ( |3.27 ) it is clear that the spherically symmetric Stephani spacetimes admitting an isotropic radiation field, for 
which c = 0, satisfy 



(7 + SR) 



0, 



which has the solutions k — Q and R{t) free (flat Friedmann), or k linearly related to R, 

h{t)^j + SR{t), 

with R{t) again free (when S = these become Friedmann models with curvature k = 47). This fixes notation for the 
rest of this thesis. 

At this point we can say that perfect fluid spacetimes admitting an isotropic radiation field for all fundamental 
observers are FLRW if and only if the acceleration of the fundamental observers is zero. This follows because the 
Stephani models with zero acceleration are FLRW (see Krasihski 1997, although it can easily be seen from ( |3.9| ): if 
u = then V = T{t)X{x^) for some functions T and X, showing that a, b and c depend on the single free function T 
and so must be FLRW by the results of the next section). Thus we have proved: 

Theorem 3 The irrotational perfect fluid spacetimes admitting an isotropic radiation field for the fundamental observers 
are Stephani models with the free functions restricted by ( 3.1{^ (or ( [A.^ jj- These spacetimes are FLRW if and only if 
the acceleration of the fundamental observers is zero. 

Spacetimes satisfying Theorem || will be referred to as IRF models. It is worth noting that all of the models admitting 
an isotropic radiation field are manifestly conformal to (part of) the Einstein static spacetime (c/. §5.1) once they have 
been transformed so that c = cz as outlined in appendix ^ From ( |3.18 ) and (A. 3) we obtain 



Vt = aiTt 1 



bi 
— 'i 
ai 

since c t = when c — c z (a s noted in appendix ^ we can assume oi 7^ 0). Changing the time coordinate via 
dt I— > aiT^F dt, the metric (|3.6|) becomes 



(l + iAr2)2 
V{z,r, t)2 



-df" 



1 



(l + iAr2)2 



(dx^ + dy"^ + dz"^) 



(3.29) 



where A = 46i/ai. The factor in braces is the Einstein static metric. In the following chapter we use this conformal 
relationship to simplify the study of the observational characteristics of these models. 



3.1.3 Symmetry and Thermodynamic Schemes. 

Having obtained the conditions ( |3.19| ) for a Stephani model to admit an isotropic radiation field it is possible to say 
immediately that all such spacetimes possess (at least) a three-dimensional symmetry group acting on two-dimensional 
orbits (just as for the general Barnes models). This follows from the work of Barnes (1998), who showed that the 
dimension of the isometry group of any Stephani spacetime is determined by the dimension d of the linear space 
spanned by the five free functions a, 5, and c: 
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1. if d = 4 or 5 (ie, at least four of the free functions are linearly independent), then the spacetime has no Killing 
vectors; 

2. if d = 3 there is a one-dimensional isometry group; 

3. if d = 2 there is a three-dimensional isometry group acting on two-dimensional orbits; 

4. if d = 1 there are six Killing vectors and the spacetime is Robertson- Walker. 



It is clear from ( 3.19| ) that a, b, and c depend on (at most) only two functions of time: /i(t) = T{t) and the constant 



function /2(t) = 1 (since Vt^Owe must have Tt ^ 0, so that these are necessarily linearly independent). Thus, c? = 2 
and the solutions have three-dimensional isometry groups as claimed. If 02 = 62 = C2 = in ( 3.19| ) then d = 1 and the 
spacetime is FLRW (see appendix ^) . 

It follows further from this and the work of Bona and Coll (1988) (see also Krasihski, Quevedo and Sussman 1997) 
that all Stephani models that admit an isotropic radiation field for all fundamental observers also admit a strict ther- 
modynamic scheme (that is, entropy and temperature functions can be found that depend on the energy density and 
pressure and satisfy the second law of thermodynamics). The converse is not true, however, since there are Stephani 
spacetimes with d = 2 (which must admit a thermodynamic scheme) that cannot have an isotropic radiation field (these 
are models for which the second independent function is not restricted to be 1). So, we have the following corollary 
to Theorem ||: 

Corollary 1 An irrotational perfect fluid spacetime that admits an isotropic radiation field has spherical, planar or 
hyperbolic symmetry and admits a strict thermodynamic scheme. 

On the subject of thermodynamics, let us mention for completeness that the thermodynamic scheme occurring most 
often in the literature is that of a barotropic equation of state. It is known that the only Stephani models with a 
barotropic EOS are precisely the FLRW models (Bona and Coll 1988; Krasihski 1997). Thus, the only spacetimes with 
a barotropic EOS admitting an isotropic radiation field are FLRW models. This also follows (when jj, + p 0) from a 
theorem of Coley (1991). Sec also Collins and Wainwright (1983). 



3.2 QCDM Models. 

The type la supernova results of Riess et al. (1998), Schmidt et al. (1998) and Perlmutter et al. (1999), which suggest 
that the expansion of the universe is accelerating, have led to an increased interest in cosmological models in which a 
significant contribution to the energy density comes from either a cosmological constant or a scalar field (quintessence 
component), which is capable of driving the expansion (Peebles and Ratra 1988; Ratra and Peebles 1988; Caldwell 
et al. 1998; Zlatev et al. 1999). In QCDM models the matter is an admixture of non- interacting cold dark matter 
(CDM), ie, dust, and a scalar field. The quintessence component may be thought of as fairies pushing the galaxies 
apartQ The fundamental observers (galaxies) are implicitly identified with the geodesic congruence of the CDM. Note 
that there is no reason a priori why the scalar field gradient (which defines a natural 'velocity' field) should be aligned 
with the CDM congruence (although it will turn out that they are aligned if the fundamental observers see an isotropic 
radiation field) . 

It is interesting to ask whether the EGS theorem can be extended to this case. We demonstrate that it can by 
proving the following theorem: 

Theorem 4 Any solution to Einstein's equations in which the matter consists of non-interacting radiation, expanding 
dust (CDM), and a scalar field (or cosmological constant), and for which the dust sees an isotropic radiation field, must 
either be an FLRW model, or have the gradient of the scalar field orthogonal to the dust congruence. 

(Note that the latter possibility means that gradient of the scalar field is spacelike, and is usually rejected as unphysical - 

although see below.) 

Proof: 

We may divide this proof into two parts: first we demonstrate from Einstein's equations in the 1-1-3 formalism that 
any energy flux component with respect to the CDM frame must be zero if the CDM observers see isotropic radiation, 
then we show that the contribution to the energy flux (with respect to the CDM frame) from the scalar fleld is zero 
if and only if the gradient of the scalar field is parallel (or orthogonal) to the CDM velocity u'^^^^ , so we deduce that 

^G. F. R. Ellis, private communication 
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the velocity fields are parallel (or orthogonal). The case where the field gradient is orthogonal to w"^^^ is probably 
unphysical, and will be rejected. Thus, the mixture of radiation, CDM and scalar field can be written as a single 



perfect fluid with geodesic fundamental congruence — u"^^^, and it follows from the results of §3.1 that the model 
is necessarily FLRW. 

Since the radiation is isotropic for the dust observers the energy-momentum tensor for radiation may be written 
in the perfect fluid form (3.5) with p — and the total energy- momentum tensor is: 

Tab = ^J.UaUb + \l-i.hab + pUaUb + (f>,a<P,b ~ gab (^'/'x'/)''' + *(0)) , (3.30) 

Radiation CDM Scalar Field 

where $(</)) is the scalar field potential (often assumed to be zero, in which case the scalar field can be interpreted as a 
stiff perfect fluid). Note that the cosmological constant case can be included by setting cj) = A =constant, $(</)) = (p. 

1. The fundamental congruence u°' is geodesic (u" = 0) because the CDM component does not interact with the 
other matter. This, in fact, implies that the rotation of must also vanish: from the momentum conservation 
equation for the radiation, we can write 

Ua = --Valn^ 



so that (using ( ^.22| )) 



= V[a?ib] = 7VrbVa] In^ = -7^ab- = ^i^bad, 

4 ^ ' 4 /i 3 

and we see that LUab = when 6 =/= 0. 
When Ua and ujab are zero, ( |3.2| ) becomes 

V[a(0Ufc]) = U[{,Va]6' = 0. 

This implies (since WaO = Wad ~ 6ua) that 

Va6l = 0. (3.31) 
(ie, the expansion is homogeneous). From the constraint equation relating the divergence of the shear to other 



kinematical quantities (2.40) we see that any energy flux component with respect to the CDM velocity field must 
vanish: 

qa = \^ae = 0. (3.32) 

This is the key step in the proof. 
2. Decomposing ( 3.30| ) with respect to m° we find that the relative energy flux component is 

= ga = -ha'u'n^ ^ ~^Va^. (3.33) 

So (7o = if (/) = (the scalar field gradient is orthogonal to and therefore spacelike), or if Vaf/* = (the scalar 
field gradient is parallel to m°). We take the latter case to be most important since the gradient of a scalar field 
is usually assumed to be timelike. 



free velocity field; it is thus an FLRW model by §3.1 



Since Va0 = Va(/>— Ua4> = —Ua<t> it is possible to writ e ( 3.30 ) as a perfect fiuid with geodesic, shear-free, rotation 



□ 



It is easy to see from the above proof that the fact that the fundamental observers correspond to dust-like matter 
was not used, only that they followed geodesies. Hence the above result applies for more general perfect fluids in place 
of the CDM component, as long as the fundamental congruence is geodesic. 

The idea of a spacelike scalar field gradient seems physically unappealing. However, such a field can (depending 
on the potential <&) satisfy the weak, strong and dominant energy conditions. The strong energy condition will be 
satisfied if and only if <&(</)) < everywhere, whereas as the weak and dominant energy conditions will be satisfied 
if $(0) > 0, although not only so. Thus, a massless scalar field (stiff perfect fluid) with spacelike gradient satisfies 
all energy conditions. It should be borne in mind, though, that scalar fields arising in cosmological contexts often fail 
to satisfy the energy conditions. This case may deserve further consideration. As can easily be seen, the scalar field 
component gives rise to anisotropic stresses in the energy-momentum tensor, so such spacetimes are not FLRW. Note 
that this theorem also applies to any number of non-interacting perfect fluids - that is, a spacetime consisting of dust 
seeing isotropic radiation and some perfect fluids will be FLRW. 
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3.3 Conclusions. 



We have proved that the irrotational perfect fluid spacetimes admitting an isotropic radiation field are Stephani models 

|)), and are FLRW if and only if the acceleration of the fundamental 



restricted by (3.19) (see also equations 

congruence is zero (Theorem ||). It follows from the fact that the constraints (3.19) depend on only two independent 
functions of time that all of the acceptable models possess three-dimensional symmetry groups acting on two-dimensional 
orbits (ie, have spherical, planar, or hyperbolic symmetry) and therefore possess a thermodynamic interpretation. We 
have also shown that spacetimes containing a mixture of radiation, dust and scalar field (QCDM models) for which the 
dust observers see the radiation as isotropic must always be homogeneous and isotropic (Theorem ^) unless the scalar 
field gradient is spacelike and orthogonal to the CDM congruence - a possibility we reject as unphysical. This result 
also relies on the geodesic nature of the fundamental congruence. 

Crucial, therefore, to the proof of homogeneity and the verification of the cosmological principle is the non-acceleration 
of the fundamental observers. Despite the intuitive appeal of cosmological models in which the fundamental observers 
are associated with a dust-like matter component, it is unacceptable to simply assume that we are geodesic observers, 
especially when such an assumption is, in principle, testable. Acceleration leaves a characteristic dipole signature in 
the redshifts of nearby galaxies that may be detectable using galaxy surveys. The physical principle underlying this 
effect is easy to see. For a set of uniformly accelerated observers ('galaxies') in Minkowski space each observer will 
see other galaxies redshifted or blucshiftcd in a dipole pattern, with the blueshifted galaxies lying in the direction of 
the acceleration, because during the light-travel time between galaxy and observer the observer's velocity has increased 
relative the velocity at emission, so that the galaxies the observer is travelling towards are blueshifted, and those it is 
travelling away from redshifted. It also follows from this that the magnitude of the dipole increases with distance, simply 
because the light-travel time from more distant galaxies is larger. In a cosmological context this acceleration dipole 
must be added to other terms contributing to the redshift of nearby objects, in particular the expansion. The method 
of Kristian and Sachs (1966) and MacCallum and Ellis (1970) gives, for any cosmological model with fundamental 
congruence u", the lowest-order term in the redshift z as a function of distance r and direction e° (a spacelike unit 
vector orthogonal to u", denoting the direction of observation): 



Hor 1 



(3.34) 



where Hq = ^Oq is Hubble's constant and the last term in brackets indicates the quadrupole introduced by the presence 



of shear. In ( 3.34| ) r can be any cosmological distance measure (area distance, for example) because for small r all such 
measures agree to first order. Note that just as the monopole (expansion) term increases linearly with distance according 
to the Hubble law, so does the acceleration dipole. This is important, because it allows the acceleration dipole to be 
distinguished from any dipole resulting from the peculiar veloci ty of our galaxy with respect to the cosmological average 
rest frame (usually identified with the CMB frame). Equation ( 3.34 ) applies in this rest frame, and any peculiar motion 
results in a doppler shift for each galaxy, which introduces an additional dipole component into the galaxy redshifts. 
This dipole is just a constant depending only on the peculiar velocity of our galaxy. A boost to the 'correct' rest frame 
can eliminate this constant component, but cannot remove the acceleration dipole because it is distance dependent. It 
is important to note in this context that the acceleration referred to here is not the same as the 'acceleration dipole' 
resulting from the gravitational attraction by the Great Attractor overdensity, which is often calculated using galaxy 
surveys (see Schmoldt et al. 1999). 

Galaxy surveys are often used to measure a possible bulk flow of the local universe, that is, the difference, if any, 
between the rest frame of the local universe and the CMB frame, which in standard cosmological models should be 
the same (see Willick 1998). A simple extension of these techniques (Clarkson, Rauzy and Barrett, in preparation) 
permits the acceleration to be constrained by observations. However, preliminary results suggest that the constraints 
on u°- are quite weak: it appears not to be possible to conclude definitively that we are geodesic observers ( 2.153 ). 
The accuracy of determinations is limited both by uncertainties in the distance estimates to galaxies as well as the 
peculiar velocities of galaxies. 

Even if the acceleration was measured to be zero, it is still necessary to show that there are no anisotropic stresses 
(Ferrando et al. 1992) before the cosmological principle can be verified. It follows from equation (31) of Ellis (1998) 
that this is equivalent to determining that the electric part of the Weyl tensor is zero. It is not clear how this may be 
achieved using observations. 

Of course, the Copernican principle (which is a vital element of EGS-type theorems, allowing the high isotropy of 
the CMB here to be assumed for other points in the universe) remains a purely philosophical assumption. It has been 
suggested by a number of authors that the Sunyaev-Zeldovich effect might be used to place constraints on the anisotropy 
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of the CMB at distant positions, but it seems unlikely that such observations will provide a definitive verification of the 
Copernican principle. Nevertheless, the arguments in favour of the Copernican principle are quite powerful, and it is a 
much weaker assumption than the cosmological principle. Note that if the acceleration is measured to be zero here, the 
Copernican principle must also be applied to give geodesic observers everywhere for the results of this chapter to hold. 

Finally, one might expect that the 'almost' version of Theorem || would lead to spacetimes that are almost the 
Stephani models of ( 3.19| ). However, when the assumption of geodesic observers is relaxed it is no longer possible to 
constrain the rotation to be small, and the class of perfect fluid spacetimes with an almost isotropic CMB may well 
include examples with distinctly non-zero rotation, unless other constraints are brought to bear. It would be interesting 
to determine the class of all perfect fluid spacetimes (including those with rotation) admitting an isotropic radiation 
field, and if it turns out that they are indeed all irrotational then the Stephani spacetimes defined by (3.19), (3.6) 
and ( 3.18| ) are indeed the complete set. 

The rest of this thesis concerns the IRF models. We derive all the relevant observational relations, for all observer 
locations, and demonstrate that the models are acceptable on observational grounds, while retaining the Copernican 
principle. 



Chapter 4 

Derivations of the Observational Relations 



This chapter presents a derivation of all the important classical observational tests of cosmology for the IRF models. 
These were presented in chapter || for the FLRW models, but the method here is slightly different. This is principally 
because the models have less symmetry (although the conformal symmetry allows us to bypass this). In contrast to the 
FLRW relations, none of the relevant quantities in this chapter have an integral in them (eg, ( 2.140| )); all the relations 
here, though, are given in parametric form, which makes them more difficult computationally. 

All the results in this chapter are new; Dsibrowski (1995) has presented observational relations for the Stephani 
models before, but only as a series expansion to first order, using the method of Kristian and Sachs (1966). 



4.1 The Stephani Models Which Admit an Isotropic Radiation Field. 



Stephani models are the most general conformally flat (expanding) perfect fluid spacetimes. They have vanishing shear 
and rotation, but non-zero acceleration and expansion. Although the general Stephani model has no symmetry at all, 
we only consider here the class possessing a three dimensional isometry group acting on two dimensional orbits, derived 
in chapter ^ The metric in comoving (and isotropic) coordinates can then be written (see § 3.1.2| ): 



(1 



iAr2)2 



Ritf 



(l + iAr2)^ 



(4.1) 



where c is the speed of light ,[| dVl^ — dd'^ + sin^ -ddip^ is the usual angular part of the metric and the function V is given 
by (using 3.25 and Xq — C^i to avoid confusion) 



V = F(r,79,<) = 



R 



2Cr cos 1? + AC 



R 



(4.2) 



R{t) is the Stephani version of the FLRW scale factor (note that we are using a different time coordinate to (3.29) 



which uses conformal time), and V{r,'d,t) is a generalisation of the FLRW spatial curvature factor (which is 1 + jkr 
in isotropic coordinates, with k = 0, ±1, but is more familiar as 1 — kr^ in Friedmann coordinates). Since k is a function 
of t the spatial curvature can vary from one spatial section to the next. In fact, it is possible for a closed universe to 
evolve into an open universe, or vice versa, in stark contrast to FLRW models (see Krasihski 1983). The models which 
admit an isotropic radiation field have k restricted by (B.27); viz. 



K{t) ^A + 2SR±2J{\A + SR) + AC^R^. 



(4.3) 



It is worth noting that this is the most general form for k, for the metric ( [4.l| ) with conformal factor (4.2) to be a 
perfect fluid; any change to k will introduce energy flux into the energy-momentum tensor - and, incidentally, make the 
expansion inhomogeneous. Note that the FLRW limit is given by i5 = C = 0. 

The four- velocity of the fluid is given directly from ( [l.l| ) since it° cx 6q in comoving coordinates and u'^Ua — —c^, 
and is 



RV 



(1 



iAr2) 



(4.4) 



^Factors of the speed of light will be kept in this chapter, and throughout the rest of this thesis, to facilitate comparison with observations. 
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The expansion is homogeneous, depending only on time, 

^RVt 



e{t) 



~2R 

3R,t 
R 



1 



1 ± 



c 



3i? 

'a" 

SR 



(- 

\R 



{\A + SRy 
0, or if A = 



A ^ 0, C ^ 



AC2i?2 



0; 



(4.5) 



(note that k must have the form 4.3). 

The acceleration, iia — u^WbUa has radial and angular components: 



In 



1 



iAr2i 



1 



Vir,^,t) 



iAr2 



2 
V 



u-g = --[lnF(r,i?,t)]^ 



Ccosi? 
2c2Csini9 



(4.6) 



where F = Iffttl^'^^ is the square root of the time component of the metric ( |4.l| ). 

The density and pressure for t hese m odels become (using equations (5) and (6) of Krasihski (1983, 1997); Eq. (3) in 
Barnes (1998), and the results of § |3.1.2D 
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3k 9 



-^=R^ 
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3k 
R^ 



3If_ 

C2A2 
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\R 



V 

i^' Vt 



(4.7) 



We see that the density is homogeneous on surfaces of constant time but the pressure is not - see also ( ^.14 ), below. 
So, for any fixed t, the spatial curvature, K(t), the expansion and the density are homogeneous, but there are pressure 
gradients that lead to acceleration of the fluid. It is interesting to note once again that the spacetime is locally FLRW 
if and only if = 0, which happens if and only if k = 0, which in turn is equivalent to the existence of a barotropic 
equation of state, p = p{iJ.) Krasihski (1983). 



4.1.1 Geometry. 

The metric {^^) is manifestly conformal to a Robertson- Walker metric with curvature A. However, if we multiply 
through by the conformal factor we see that the actual spatial curvature is time- and position-dependent and is given 
by the curvature factor K{t) in V{r, t) in the spherically symmetric case, but is quite complex in general, with variations 
depending on angle and distance. When C 7^ the spatial sections have a hyperbolic geometry - and thus a center of 
symmetry. However the spatial sections foliate in such a way as to make this 'center' wander around. Indeed there is 
no reason to assume that this central worldline will be causal, and an observer who happens upon this worldline may 
be able to see it in their past. It is possible in principle to write the metric in coordinates to make this apparent, but 
we do n ot do this here. 

(4.1) and (5^) are the usual forms in which the Stephani metric is presented (see D§,browski 1993, 1995; Krasihski 
1983, 1997). However, for our purposes they are not the most advantageous forms. The conformal geometry of the 



models is most easily studied by changing from the stereographic coordinate, r, to the 'angle' coordinate, x (see figure 4.1 



and equation (5.15) of Hawking and Ellis, 1973), appropriate to the value of A. Furthermore, for models with closed 
spatial sections (which will be our principal concern here) it is more convenient to choose a radial coordinate that is 
better able to reflect the fact that light rays can circle the universe many times. In such models the spatial surfaces 
have two centers or symmetry, r = and r = 00. Physically, there is nothing extraordinary about the point r — 00: 
it is not infinitely far away from the centre, and it is quite possible for light rays to pass through it. This last point is 
particularly important for subsequent discussions, so we make the coordinate change 

2 



, X 
tan — 
2 



(4.8) 
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Center of symmetry 



Coordinate radius r 




Hypersurface.of 
constant time 



Projection of galaxy onto plane 



Observer looks to 
galaxy 



Figure 4.1: Stereographic projection of a sphere showing the relationship between the coordinates r and x- Note that one can 
envisage the angle x increasing steadily round the sphere ad infinitum. 



Then r — > oo as x i"- As a coordinate x is restricted to the range < x < i"- However, it will prove convenient to use 
X not just as a coordinate but as a parameter along light rays. In the latter role its value can increase without bound. 
Strictly speaking, we should distinguish these two uses, but it should not lead to confusion. The absolute value is taken 



in ( p~q ) so that, when x increases beyond tt, r remains positive. Using (4.8) in the metric (5.4) gives 

1 ( o ,.o R{tf 



where 



W{x,^,tf 



W{x,^,t) ^ cos^^V {r{x) J, t)R{t) 



-idx'+sin'xdn') 



(4.9) 



- 1 



cos^ - + -V^ sin^ - - 2^-7^ sin y cos -d 
2 A 2 y^A 



(4.10) 



in general; and 



W{x, t) = cos^ ^V{r{x), t)R{t) = cos^ ^ + 



X , '«(*) . 2 X 



2 ■ "~ 2 ■ A 2' 

in the spherically symmetric case. Singularities in the conformal factor {\/W) correspond to spatial and temporal 
infinity. 

We can easily calculate the acceleration in these coordinates. Again it has only radial (%) and angle {d) components: 



,2W^ 

w 



A simple calculation shows that the acceleration scalar, which we will need below, is just 



(4.11) 



(4.12) 



For completeness we mention that when the IRF models are conformal to an FLRW spacetime with hyperbolic 
geometry (A < 0) the coordinate transformation is obtained from that just given by replacing trigonometric functions 
with their hyperbolic equivalent. 
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4.1.2 Non-Central Observers. 

Since the purpose of this thesis is to provide insight into the observational characteristics of the IRF models for all 
observer positions, we must find expressions for the distance-redshift relations and other observable properties of the 
models from any point in the spacetime. For a general inhomogeneous metric this is far from trivial, but the IRF models 
have features that make this problem tractable. In particular, they are conformally flat. As has already been noted, in 



equation (O) the part of the metric in braces is exactly the form of the FLRW metric in 'angle' coordinates, so that 
the IRF models are manifestly conformal to the (homogeneous) FLRW spacetimes. This means that there is a group of 
transformations acting transitively on surfaces of constant time that preserve the form of the FLRW part of the metric 
(but not the conformal factor). If the conformal spatial sections are closed (3-spheres), flat or hyperbolic (according to 
the value of A), the transformations are rotations, translations or 'Lorentz transformations', respectively. After such 
a transformation the metric will have the form of an FLRW metric centred on the new point, multiplied by a modified 
conformal factor. 

As will be shown in § ^.l.l| , we will be dealing exclusively with closed models (A > 0), and so will concentrate on 
this case. 

( — — The Spherically Symmetric Case 

In the spherically symmetric case, to find the coordinate transformation to a non-central position, we perform a rotation 
of the spatial part of the metric, moving the origin (x = 0) to the point x — the observer's position in what 

follows). In appendix]^ we derive this transformation. The old x is given in terms of the new (primed) coordinates 
by (p3.6[), which we reproduce here: 



cosx = cos?/; cos x' — sin'f/'sinx' cosi?', (4-13) 

where 1}' denotes the direction of the new location. The conformal factor now becomes, dropping the primes on the 
coordinates ^ ^ 

W — > W{x,i},t) ^ 7: + -r) + 7: — -r) (cos cosx — sin?/)sinxcos??), (4-14) 



2 V A/ 2 V A. 
while the rest of the metric ( [4.9| ) stays the same, but now in terms of the new coordinates. 

This transformation makes the study of our inhomogeneous models significantly easier, and allows us to find exact 
observational relations valid for any observer. 

C^O - The General Case 



When C 7^ 0, things are not so simple, and the form of the conformal factor (4.10) is not general enough. This is because 




in (4.2) we took the direction of the center of symmetry to be aligned with the z axis. We desire the generality of the 
case above: we wish to be able to move the observers to any location in the spacetime, regardless of the direction to the 
center. In the spherically symmetric case we only had to worry about the radial coordinate; in the axially symmetric 
case we need to worry about radial and an angle coordinate. Instead of performing a general coordinate transformation, 
we instead define our direction to the center to lie somewhere in the y = plane, in the direction given by 

(4.15) 

y cos^ J 

Now we may write 

V = V(r,i},t) = - (l + -K(t)rA +4:C^--2((xsm£_+zcosn. (4.16) 

it \ 4 / K 

In order to find W, we must first note that in terms of the mixed coordinates Xi z we have 

W ^-[1 + — + 4C^— + - 1 - — + 4C^ — cosx - 2-^ ^ 1 , ^ \ 4.17 

In order to use the transformations in appendix |b| we note that the cartesian coordinates used there are not the same 
as the ones used here, but are related by 

^/Aa;* 
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(this amounts to the transformation between the FLRW and stereographic coordinates) and the rotation is in terms of 
the i% so that we have in the new primed coordinates 



1 + iAr2 



1 + iAr2 



— = sin y' sin -d' cos ip 

(sin ip cos x' + cos -0 sin x' cos {)'). 



(4.19) 



We now have 



If K 



1 ( K R 
+ - < 1 — — + 4^^ — J- (cos ijj cos X — sin i/) sin x cos i?) 



CR 

2—= (sin x[sin ^ sin i9 cos ip + cos ^ cos ■0 cos z9] + cos ^ sin ■0 cos x) ■ 
V A 



(4.20) 

Li 

Note that W is of the form 

W = A + Bsin'd cosifi + C COS'S ^ A + Bx + Cz, (4.21) 

where A, B, and C are functions of time only, then we can rotate in the x — z plane (by tan^^(_B/C)) so that W depends 
only on z, and we have in this new coordinate system 



W = A+ + cosz9. 



(4.22) 



This shows that, as in the spherically symmetric case, we have only a dipole moment in W. 

To summarise, we have changed the coordinate system in order to center our coordinates on any point in the 
spacetime. In these new coordinates the metric is of the form 



1 



Rit? 

A 



(4.23) 



where 



W = + $_ cos cos X 



2CR 



cos ^ sin '0 cos x 



/ C^L^ QR^ 

sin X cos-dJ sin^ tp + 4 ^ (1 — cos^ ^ sin^ -0) + 2 <^os <^ sin 20, 



(4.24) 



and 



1 , ^ 



(4.25) 



4.2 Observational Relations 



So far we have considered only some of the 'global' physical properties of IRF models, but to really assess their potential 
viability as cosmological models it is necessary to confront them with observations. In this section we derive the 
distance-redshift relations that form the basis of the classical cosmological tests which are used later to compare them 
with available observational constraints to see whether any regions of parameter space are capable of providing a fit. 

Deriving the observational relations (redshift, angular size or area distance, luminosity distance and number counts) 
means relating the coordinates and metric functions to observable quantities. This requires knowledge of the observer's 
motion (4- velocity), which can, strictly speaking, be specified independently of the background geometry. However, the 
IRF models contain perfect fluid, so we will identify the observer's motion with the fluid velocity. We are not obliged to 
do this, and, given the strange form of the matter, it might be thought advantageous to instead assume that observers 
(ie, galaxies) constitute a dust-like test fluid moving freely through the spacetime whose geometry is determined by the 
exotic matter. It should be clear from chapter ^ that if we were to make this assumption a large dipole anisotropy in 
the CMB would result (although the dipole in Hq would be eliminated - see (5.3C) with ?i = 0) because such a flow will. 



in general, have a significant velocity relative to the IRF fluid flow, which, it will turn out, is very nearly in the rest 
frame of the CMB everywhere. This should not be a surprise because chapter ^ was devoted to proving just this. 



CHAPTER 4. DERIVATIONS OF THE OBSERVATIONAL RELATIONS 



52 



4.2.1 Redshift 

Probably the most important distance measure in cosmology is redshift. In general, it is no simple task to find analytic 
expressions for the redshift in any cosmological model; derivations usually rely on symmetries of the spacetime or other 
simplifying factors to solve the equations of null geodesies. Here we can take advantage of the conformal flatness of 
Stephani models (strictly speaking, of the fact that the IRF models are manifestly conformal to FLRW cosmologies) 
using the results of although we can also derive the redshift formula as a time dilation effect which we do here. 
When the true spacetime is conformal to a spherically symmetric spacetime the radial null geodesies connecting any 
point with an observer at the centre (of the conformal part) are obviously purely radial (since their paths are not affected 
by the conformal factor). They are therefore given (in terms of coordinates r and t with respect to which the spherical 
symmetry is manifest) by some function toi^r^, t^) relating the time, to , that the light ray is received by the observer, 
to the time of emission, t^, for an object at radius r^. This is just the lookback-time relation. Redshift, as the ratio of 
proper time intervals at the observer to proper time intervals at the emitter, is then given by 



dTo dto dTo 
dto dTE dto 



dto drs 



dto dts 
dtp. dTf_ 



1 



dt. 



i% \ dr„ 



dt 
dt. 



(4.26) 



(When the coordinates r and t are comoving - = Q - the r-derivative term disappears.) This will provide an analytic 
expression for the redshift whenever the lookback-time equation can be integrated. For the IRF models: 



0, 



1 500 



1/2 



= W 



(4.27) 



and the lookback time can be derived directly from the metric: on the past null cone of the observer ds — Q — dd — dip, 
leading to an expression for dx/dt, which, when integrated, gives 



x{t) = cVA 



dt' 



where T is the coordinate age. Differentiating this with respect to t at fixed x then gives 

dto _dT _Ro 
dt^ dt Re ' 



(4.28) 



which, together with (4.26), (4.27), and (4.23) we find 



l + z(V;,x,^9) 



i?o W{^,x,^;t) 
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R{t) 



(4.29) 



where i?o — R{T), Wq = W{ip,T) and t and x are related by equation (4.28). 

A more elegant approach is to use the result of §^ For the IRF models gab is an FLRW metric, the conformal factor 



is — 1/W (see (4.23)) and is the usual FLRW comoving velocity field. The well-known expression for redshift in 



(4.30) 



FLRW spacetimes, 1 + z = Ro/Re, then gives 

Wo Rb' 

(Alternatively, we could change to conformal time ( 5.1(]| ), making the metric manifestly conformal to the Einstein static 
spacetime, (5.11), for which the redshift is zero. The conformal factor is then e*^ = R/ (A^^'^W), which also gives (4.30).) 
Now, if we take our large expression for W (4.24) and write it as 



(4.31) 



where 
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(4.32) 
(4.33) 
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then we see that 

all inhomogeneity simply results in a dipole variation in z around the sky. 
4.2.2 Distance-Redshift Relations 

For metrics with spherical symmetry about the observer the angular size (and area) distance is given directly from 
the coefficient in front of the angular part of the metric, because symmetry ensures that for radial rays t9 and ip are 
constant along the trajectory. For our models we do not have spherical symmetry about every observer, but the metric 
is everywhere conformal to a spherically symmetric metric, as can be seen from (|^). Since null rays are not affected 



by the conformal factor they also remain at fixed and Lp, so we can again obtain the angular size distance, ta, from 
the coefficient of the angular part of the metric: 

rA{^,X,0) = ^J^L^Isinxl (4.35) 



(again, x a-nd t are related by ( 1.2§| )). The absolute value of sinx is taken to keep rA positive - which we are perfectly 



entitled to do since only sin x appears in the metric (4.23). We can find now, for the first time, the exact angular 



size distance relation parametrically by combining equations (4.35) and (4.29), which bypasses the rather cumbersome 
method of Kristian and Sachs (1966), used in Dabrowski (1995) to find a series relation. This is valid for any null- 
connected points in the spacctime. As far as I am aware, nobody has done this for an inhomogencous model before. 

The other classical tests can now be written down. Luminosity distance, r^, is related to va by the reciprocity 
theorem: 

rL^il + z)\A, (4.36) 

see MacCallum and Ellis (1970) and Ellis (1998), or § pl] . This then allows the magnitude-redshift relation to be 
determined in the usual way: the apparent magnitude m of an object of absolute magnitude M is given in terms of the 
luminosity distance by 

m-M-25 = 51ogiorL. (4.37) 
4.2.3 Number- Count-Redshift Relation 

Finally, for completeness, we discuss the number count-redshift relation, although we do not use it here. The only 
important consideration is the identification of some sensible comoving number density distribution, ndx)-, for the 
observers (fluid particles). This is not necessary for homogeneous spacetimes because it is natural to simply take ric 
also to be homogeneous (and independent of time, assuming no evolution). Furthermore, when there is fluid pressure 
ric cannot be directly related to the energy density of the matter, as is possible with a dust, and, in general, there need 
be no obvious choice for - the only constraint it must satisfy (again assuming no number evolution, such as would 
be caused by mergers, for example) is particle number conservation: 



Since the metric (4.9) is in comoving coordinates we can satisfy this by simply assuming that Uc is independent of time t, 
being given by whatever distribution we specify a,t t — 0. Moreover, the fact that the IRF models become homogeneous 
(W ^ 1) as t — > means that it makes sense to define ric to be independent of spatial position too. We are then back to 
the same situation as with FLRW models, except that now observations of galaxy numbers constrain comoving volume, 
but not proper volume (note that the proper density of particles is not constant, because the conformal factor modifies 
the proper volume by different amounts at different positions and times), and so are less directly a constraint on the 
radial component of the metric (which determines proper distance). The number of particles within some radius x of 
any observer is simply proportional to the volume of the FLRW 3-sphere of curvature A out to that radius: 

^(^^ ^ ^3/2 J ^^"^ ^ ^ ^ -^j^[X~sinxcosx) 



exactly as in FLRW models. The difference from the FLRW N — z relation comes when we use (4.29) to relate number 
counts to redshift: z{x) will be different to that for FLRW spacetimes. Nevertheless, it is clear from this discussion 



CHAPTER 4. DERIVATIONS OF THE OBSERVATIONAL RELATIONS 



54 



that the fluid particles (observers) in the IRF spacetime are distributed uniformly with respect to the FLEW volume on 
the spatial sections. Since these hypersurfaces of constant time are just spheres in the FLRW metric we know that the 
Copernican principle can easily be applied: the probability that an observer lies in some region of space is exactly the 
volume of that region with respect to the FLRW volume (divided by the total volume of the sphere). 

4.3 Discussion 

In this chapter we have, for the first time, derived all the relevant classical cosmological distance-redshift relations exactly: 
that is, we have given relations connecting observable quantities, which can now be used to pit the models against real 
data. Moreover we have derived these relations in such a way as to be valid for any observer in the spacetime. We 
have not discussed, though, observations connected with a blackbody (ie, the CMB); this will be done in the following 
chapters (in fact, the temperature evolves simply as 1/(1 + z)). 

In the following chapters wc use these observational rc;lations (with a specific form for the scale factor in order 
to integrate the lookback time relation) to demonstrate that these models are quite suitable as a cosmological model 
from an observational point of view - regardless of observer position. This is important because previous studies of 
non-standard inhomogeneous cosmologies have always assumed a privileged observer position (at a center of symmetry), 
or in a small neighborhood of such a location. 



Chapter 5 

Observational Characteristics With C = 



In chapter |j we derived expressions for all the important observational relations used to test a cosmological model, 
and we took the time to derive them for any location in the spacetime. Throughout the rest of this thesis we wish to 
demonstrate that these observational aspects of the models are compatible with present observational constraints. In 
this chapter we restrict ourselves to the spherically symmetric class, C = 0. 

These models are far more general than we will need because there are no restrictions on the function R{t). In 
fact we need a specific form for R in order to find the lookback time relation (4.28). We want a family of models with 
a small number of free parameters to simplify the analysis and produce graphs etc., so we restrict attention to the 
spherically symmetric two-parameter family derived in §IV of D§browski (1993). These are just an extension of Model I 
of Dg,browski (1995) (in these papers the parameter b in equation (pTl|), below, was set to 1). They were chosen largely 
because of their simplicity, and because they proved to be the most useful models in Dabrowski and Hendry (1998). As 
described in Dabrowski (1993), we choose the metric functions in terms of this new t to be of the form: 



R{t) ^ ct{at + h) 

K(t) = A - —R{t) (= 1 - Ritfjc^), 



(5.1) 
(5.2) 



where a M —cS and b are the free parameters and 



l~b\ 



and C = 0. The metric is then, 



R{t)^V{r,ty 



R{tf 



(l + iAr2)2 



(5.3) 



(5.4) 



We will henceforth refer to these models as Dabrowski models. Note that in Dabrowski (1993, 1995) the quadratic 
dependence of R{t) also includes a constant term, d. We restrict attention to models which have a big bang (a density 
singularity in the language of D;ibrowski 1993), which means that R{t) must have a root. We can therefore choose the 
origin of time so that d — Q. Furthermore, we require that after the big bang (and before any big crunch) i? > 0, which 
forces 

6 > 



We make the coordinate changes used in chapter^; viz. 

1 

so that the metric becomes 



, X 
tan — 
2 



ds' = 



1 



wix.tf 



-c^dt^ 



A 



(dx'+sin^ Xdn')}; 



and we can use (5.2) to simplify the expression for W, resulting in 



^(X.)^l-^sin^f. 



(5.5) 
(5.6) 
(5.7) 

(5.8) 
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This has no dependence on t9. The conformal factor {1/W'^) is non-singular for aU x if a < 0. When a > singulari- 
ties W = correspond to spatial and temporal infinity, and indicate that the universe has 'opened up'. 



We transform to a non-central location, as discussed in §4.1.2, so that W transforms as 



W — s- W(x, 'ip,'&]t) = 1 — ^"-^(^) /-^ _ cos cos X + sin V' sin X cos I?) , (5.9) 

cA 

while t he r es t of the metric retains its original form (but now in terms of the new coordinates). 



In (5.1)- (|5.3|) we have retained factors of the speed of light, c, to facilitate comparison with the references given 
above and with observations. The units we will use are as follows: [c] = km s^^, r is dimensionless, R is in Mpc and 
[t] — Mpc s km~^= [1/Ho], so that [a] = km s^^ Mpc^^= [Hq] and b is dimensionless. Note that these units are slightly 
different to those used in Dabrowski (1995) because the parameters a and b in that paper contain a factor of c (so 
[b] = [c], etc.). This explains the appearance of c in ( ^.l|) . 

We will use T to denote the coordinate time of a specific epoch of observation along some observer's worldline (ie, the 
coordinate age of the universe), again in Mpc s km^^, and r to denote proper time along a particular fiow line. When 
we state ages they will generally be given in Gyr: TQy^. w 978T. 

We will impose constraints on the value of Hq, age, size (the meaning of which will be explained below) and the 
anisotropy of the microwave background|], leaving the wealth of data available from galaxy surveys and high-redshift 
supernovae for consideration in a later paper; the complexities involved in interpreting such data and applying it to 
inhomogeneous cosmological models require separate treatment. 

The task now is to limit a, b, and T using present observational constraints. A full discussion of each constraint 
is made in the following sections, and they are followed by exclusion diagrams showing the regions of parameter space 
for which a and b give a plausible cosmological model for all observer locations in these models. First off though we 
demonstrate that the choice of scale factor ( ^!l| ) gives an immediate solution to the horizon problem (cf. Rindler, 1956). 

In the case of the Dabrowski models, for any time t the spatial sections are homogeneous and isotropic and their 
geometry depends on the sign of K{t). If, at some point during the evolution of the universe, t = —{b± l)/2a, the 



curvature changes sign, as can easily be seen from (5.2) and ( |3.lD . For example, if \b\ < 1 and a > then the spatial 
sections are closed (k > 0) for t < {1 — b)/2a. The spheres increase in size until t passes through t = {1 — b)/2a, 
when the universe 'opens up' and acquires a hyperbolic geometry. This does not happen in FLRW models, where the 
spatial curvature, fc, is fixed. The distinction between the time-dependent true geometry (k) and the fixed conformal 
geometry (A) should be borne in mind throughout what follows. As the universe opens up and the sections become 
hyperbolic the coordinate x represents a conformal mapping from a hyperbolic surface onto a sphere: spatial infinity 
will then correspond to some finite value of x < For a more detailed explanation of this see theorems 4.1 - 4.4 of 
Krasihski (1983). 

When a > 0, there is a difficulty in applying the Copernican principle, which is central to this thesis. This is because 
once the universe has opened up it becomes infinite in extent; applying the Copernican principle becomes extremely 
difficult, because there is no way to consider 'all' observer locations. Even worse, 'most' observers cannot be defined 
in any obvious way. This is not a problem in the homogeneous open FLRW models simply because it is homogeneous; 
a priori every location is equivalent to every other. In the Stephani case this is not true: we are attempting to show, 
through what one can observe, that all locations are equivalent. Obviously this cannot be done for a > 0. This isn't 



a problem though, because in §5.1.1 we reject the case a > because it fails the energy conditions. However see the 
discussions in chapter |^. 

5.1 Lookback Time and the Horizon. 

If, in equation ( [4.9|) , we absorb a factor of R^(t)/ IS. into the conformal factor and change to conformal time, r/, defined 
by drj — c\^dt/R{t), so that 

t(r,) =&/(e-'"'/^-a) (5.10) 
then the resulting metric is clearly conformal to the completely homogeneous metric 

dsl = -dif + dx^ + sin^ xdrj^ (5.11) 



which is precisely the Einstein static spacetime (ESS - Hawking and Ellis 1973, §5.1-5.3). Equation (5.10) implies that 



the big bang (t = 0) happens ai r] = —oo. When a < the big crunch {t = —b/a) occurs ai rj — +oo, so the Dabrowski 

^This may come as a surprise, given that the models were derived specifically to have an isotropic radiation field. The CMB has physics 
behind it though, and consideration of this produces anisotropy. 
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models with A > are conformal to all of the ESS when a < 0. If a > the models are not conformal to all of the 
ESS because then rj will be bounded above and, what is more, the coordinate x will not take on all values once the 
universe has opened up (the no-horizon argument outlined below applies for a > too, though, because it only relies 
on the Dabrowski models being conformal to a region of the ESS that is unbounded below for all x)- Since the null 
structures of conformally related spacetimes are identical we can find the null geodesies of the D;ibrowski models directly 
from those of the ESS, which can be derived trivially from the metric ( ^.11 ): the past null rays from x — V — Vo 



satisfy t] = rjo — x- That is, they are straight lines at 45° in the (x, ?7)-plane. If we represent the ESS (and the Dabrowski 



models) on the Einstein cylinder (see figure 5.1), where x is shown as a circle (0 < x ^ 27r = 0), the null rays just 
circle around the cylinder indefinitely into the past, going round infinitely many times regardless of the initial rjQ (or to). 
Thus, the whole of the big-bang surface is contained within the causal past of every point of both spacetimes, and there 
is no horizon problem for the D§,browski models. 

We can see this in another way if we calculate the lookback time in our models (ie, the time, t, at which a galaxy at 



some p ositi on x e mits the light that the observer sees now at time T), which we can do using (5.10) (or directly from the 



metric (4^) - see 4.28 ). From ( 5.1C ) and the equations for null geodesies in conformal time, t] = ?yo(2^) ~ X (with rjoiT) 
given by the inverse of ( 5.10| )) we get 

, ^ bT , ^ 

t{x) = ^ . (5.12 

(aT + 5) exp (6x/\/A) - aT 

Note that this function is continu ous through x = tt. Again, t —^ ii and only if x ^ oo- We can visualise this using 



the 'Dabrowski cylinder' of figure 5^ (like the Einstein cylinder but with t instead of rj as the time coordinate). The 
null rays still asymptote to the big-bang surface, but this diagram reflects the way that in reality the spatial sections 
shrink as time decreases to zero: as the sections shrink (going back in time) light rays can travel all the way round the 
universe in a shorter and shorter time, and if the spatial sections shrink sufficiently slowly as we move back towards the 
big bang this travel time can approach zero without the light ray hitting the big-bang surface. 

This is all in sharp contrast to FLRW universes and de Sitter space, which are conformal only to parts of the Einstein 
cylinder (they are bounded below at 77 = 0): at early times the particle horizon is finite and contains only a small part of 
the big-bang surface, so that widely separated points can share no common influences. See figures 17 and 21 in Hawking 
and Ellis (1973). 

The existence of a horizon is directly related to the rate at which the (FLRW or Stephani) scale factor grows for 
small t. li R t" at early times, then conformal time goes as 

.-/^-^i-" (or logt for a. 1). 

So, when a < 1 (as it is for dust FLRW models) 77 = at t = 0, the spacetime is conformal to only part of the Einstein 
cylinder and a horizon exists. Otherwise, as t — > 77 — > —00 and the horizon is absent. 

Other methods which lead to the solution of the horizon problem include using the LTB models; Celerier and 
Schneider (1998). See also Rindler (1956) who first reported the horizon problem. 

5.1.1 Matter Content and Energy Conditions. 

The Stephani models do not have an equation of state in the strict sense, with the relationship between pressure and 
energy density being position dependent. Along each flow line, however, there is a relation of the form p = p{p.) - see 
Krasinski (1983). The particular models we are using have an equation of state at the centre of symmetry (or everywhere 
in the homogeneous limit a — + 0) of the (exotic) form p = — l/i. The matter content of these models with regard to the 
natural (comoving) velocity field is a perfect fluid with energy density 

SttG 3 /r 1 o\ 

^~Ai=^7TTi> (5-13) 



my 



and pressure given by 



2^2 y(r,t) 



.-,cMl + ^sin^|), (5.15) 
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Figure 5.1: Lookback time and tlie past null cone on tlie Einstein cylinder. This shows the null structure for any model conformal 
to the Einstein static spacetime, including the Dabrowski and FLRW models. In FLRW models, the big bang occurs at 77 = 0, 
whereas in Dabrowski models, the spacetime extends to 77 = —00. 
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Observer 





Big Bang 



Figure 5.2: As in figure 5.1, but now using the actual Dabrowski model time coordinate, t. 



where (5.13) has been used to express the t-dependence of pressure (which arises through the appearance of R{t) 
in V{r,t) - see (4.2) and (5.2)) in terms of the density. In other words, we have an 'equation of state' of the form 



(5.17) 



The appearance of —-^fJ- as the dominant contribution to the equation of state immediately suggests a quintessential 
or scalar field model (see Frampton 1999; Liddle 1999; Coble, Dodelson and Frieman 1997; Liddle and Scherrer 1998), 
although it has been shown (Vilenkin 1981; Dabrowski and Stelmach 1989) that cosmic strings also give rise to this 
EOS. The interpretation of the Stephani matter content as a perfect fluid is by no means required; other interpretations, 
such as a scalar field, are equally valid, and the fact that there is no true EOS might even suggest a two-component 
interpretation. It is certain from chapter |^ that they will admit a thermodynamic scheme. We will discuss this in detail 
in a future paper. For the moment, however, this is as far as we will go to provide a physical motivation for the matter in 
the Dabrowski models: in this paper we are only interested in the observational consequences of the geometry; although 
when we impose the energy conditions below we show that the matter is certainly not obviously unphysical. 

We note that the Bianchi identities lead to the following conservation equations (Ellis 1998; Wainwright and El- 
lis 1997) 



Vae = -Ua(e + f) 
fi = -3i/(e+ |)Ai 
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the second of which is exactly the form of the energy conservation equation for FLRW models with a '7' equation of 
state - note also that it shows that /t is not homogeneous. 

We can see that there are singularities of density and pressure as R{t) — > (ie, at t = 0, —b/a), which correspond to 
the big bang and crunch for these models (the metric becomes singular at these points too) - see Dabrowski (1993). We 
can also have a finite-density singularity, where only the pressure becomes singular. This happens when r 2/^/—A. 
Such infinite pressure is clearly not physical, so we can reject models with A < 0. For models with A = 0, it is difficult 
to compare them directly with A > models due to the different geometries of the spatial sections, so we will not 
consider them in this thesis, and we are left with A > 0, ie, 6 < 1 - see (5.3). As we explained in §4.1 the natural 
assumption that R{t) > after the big bang ensures that 6 > (A < 1). 

Having calculated the pressure and density of the fluid we can now investigate its physical viability through the 



energy conditions. It is more convenient to use the original stereographic coordinates for this (ie, expression (5.14) for 
the pressure), since we wish to consider all values of A, until we find reasons to the contrary. The weak energy condition 
states that /i > and p + fjic? > 0, whereas the strong energy condition is equivalent to 3p + n<? > (see, for example, 
Wald 1984 for a discussion). The weak energy condition does not constrain our models at all: it implies that > 0, 
but this is always true since V ^ only at spatial (and temporal) infinity (even though the coo rdinates themselves may 
be finite). The strong energy condition, however, implies that K{t) > A for all t. From (5^) we can see that this is 
equivalent to a < (since i? > 0), so the models must have a big crunch {R{t) is an 'upside-down' quadratic). 

The dominant energy condition is more interesting. It states that \p\ < fic^ , from which ( ^.13 ) and (^.16|) immediately 
give 

1 V{r,t) 



< 



3 1 + iAr2 



< 1. 



The left inequality requires only that A > 0, which rules out infinities in the pressure (finite-density singularities). The 
inequality on the right says that for all t and r 



{K{t) - 3A)r2 



< 



(5.18) 



must hold. This condition is always true for a > (see (5.2)), as long as A > 0. When a < 0, r is unbounded (for A > 0, 
because k is then positive - see (4.2)), so the left hand side of ( 5.18 ) must always be negative, ie, K{t) < 3 A. It is easy 
to show that R{t)/c < —P/Aa, so, from (5.2), 

K{t) < 1. 



Then k < 3A for all t provided 



& < J- « 0.82. 



(5.19) 



For b larger than this the dominant energy condition will be broken at some time, in regio ns of the unive rse at large r. We 
will not consider further the intricacies of this. A glance at the exclusion diagrams, figures 5.12 and 5.13 , shows t hat (5.19 ) 
does not eliminate a significant area of the allowed region. In light of this we will, for the moment, overlook ( ^.16| ) and 
investigate the properties of all models with < 5 < 1. 

To summarise: we have used basic physical requirements, such as the occurrence of a big bang and the avoidance 
of pressure singularities, and energy conditions to restrict the ranges that the two parameters a and b (or A) can take. 
The results are: 

a < 0, < 5 < 1 (ie, < A < 1) (5.20) 

(we reject 6=1 for simplicity, as explained above, and we refrain from invoking ( ^.19 ) until § ^.5| ). Of course, we are 
not forced to accept the energy conditions. Although they seem physically very reasonable conditions to impose on any 
form of matter there are many examples of cosmological models based on matter that does not satisfy them (quantum 
fields, for example, can exhibit negative energy density). However, we will assume their validity and in most of what 
follows we only investigate the properties of the models satisfying ( 5.2Cl| ). 



5.2 Constraining the Model Parameters Using Observations. 

So far we have considered only the 'global' physical properties of Dabrowski models, but to really assess their potential 
viability as cosmological models it is necessary to confront them with observations. In this section we derive the distance- 
redshift relations that form the basis of the classical cosmological tests and compare them with available observational 
constraints to see whether any regions of parameter space are capable of providing a fit. We will impose constraints 
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on the value of Hq, age, size (the meaning of which will be explained below) and the anisotropy of the microwave 
background, leaving the wealth of data available from galaxy surveys and high-redshift supernovae for consideration 
in a future paper; the complexities involved in interpreting such data and applying it to idealised cosmological models 

require separate treatment. 

To recap from chapter |4| when the change of coordinates is made to an arbitrary location, we have the metric (4.23) 
with W restricted to 

W W(x, -0, 1?; i) = 1 — ^^^^^1(1 _ cos cosx + sin V'sinxcosi^). (5-21) 



Using the expression for the lookback time (5.12) and the expression for the redshift (4.30 and our simplified 
expression for W (|5.9|) we find 



l + z(V',X,^) 



Wa R{t) 

Ro f 1 2a 2a . 

Wo \ Rjt) ^ ~ ^smVsmxcosi^ } ; (5.22) 



where i?o = RiT), Wq — W{iPtT), showing that, for objects at any fixed Xi the inhomogeneity of u nivers e manifests 
itself in the redshift as a pure dipole in angle around the sky (cos d term) . This will be important in § ^.2.4 . 



For metrics with spherical symmetry about the observer the angular size (and area) distance is given directly from 
the coefficient in front of the angular part of the metric, because symmetry ensures that for radial rays and ip are 
constant along the trajectory. For our models we do not have spherical symmetry about every observer, but the metric 



is everywhere conformal to a spherically symmetric metric, as can be seen from (4.9). Since null rays are not affected 
by the conformal factor they also remain at fixed "d and ip, so we can again obtain the angular size distance, r^, from 
the coefficient of the angular part of the metric: 

rAi^.X,^) = y=-^^^^—-\sinx\ (5.23) 

again, x and t are related by ( ^.12|)). W e can find, for the first time, the exact angular size distance relation para- 



metrically by combining equations ( [4.35| ) and (4.29), which bypasses the rather cumbersome method of Kristian and 
Sachs (1966). This is valid for any null-connected points in the spacetime. 
Luminosity distance, r^, is related to ta by the reciprocity theorem: 

ri = (l + z)2rA, (5.24) 

see MacCallum and Ellis (1970) and Ellis (1998). This then allows the magnitude-redshift relation to be determined in 
the usual way: the apparent magnitude m of an object of absolute magnitude M is given in terms of the luminosity 
distance by 

m-M-25 = 51ogiorL. (5.25) 

The task now is to limit a, b, and T using present observational constraints. A full discussion of each constraint is 
made in the following sections, and they are followed by exclusion diagrams showing the regions of parameter space for 
which a and b give a plausible cosmological model for all observer locations in these models. 

5.2.1 Hubble's Constant. 

The expansion rate of the universe has been measured with reasonable accuracy. Hubble's constant is believed to lie in 
the range 50 S Hq < 80 km s~^ Mpc~^, and we will u se th ese limits to constrain the D^browski models. The Hubble 
parameter for these models is independent of position ( [4.5| ): 

R(t)t Rt(T) 
H M- = —-^ =^ Hq^ '} ' (5.26) 
3 R{t) " R{T) ^ ' 

We can use this to place constraints on the time at which observations can be made at any position: for our models H 
decreases monotonically, so it will only lie in the observed range of Hq for some range of T. For any observer with 
coordinate age T, we require 
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Usually, for simplicity, we will choose a specific value for Hq (almost invariably that which produces the 'worst case'). 



Then we can solve (5.26|) for T 



However, when Hq is actually measured, it is not necessarily equal to the expansion rate. What is measured in 
practice is the lowest order term in the magnitude-redshift relation, which gives the measured Hubble's constant, -ff™: 



til 



(5.28) 



(where fc" denotes the wave-vector of the incoming photons); see MacCallum and Ellis (1970). Equivalently, we can 
consider the gradient of the redshift-area distance curve at the observer (c/. Humphreys et al. 1997). The covariant 
derivative of the velocity field that appears in ( 5.2^ ) can be decomposed in terms of the expansion, acceleration, shear 



and rotation of the flow. In the case of the Stephani models these last two are zero, so 

VbUa = Hhab - UaUb- (5.29) 

Therefore the measurement of depends upon the acceleration of the observer. If we measure the magnitude-redshift 
relation for objects in some direction, then, comparing equation ( ^.2^ ) with the expansion given by ( 5.26| ), we find 



11 

H"' (9) = Ho - - cose, (5.30) 
c 

where ii is the acceleration scalar and 9 is the angle between the acceleration vector and the direction of observation 
(which is opposite to the direction in which the photons are travelling) . Since the acceleration is non-zero in the Stephani 
models there will be a dipole moment in H™. The size of this in the Dsibrowski models is given directly from ( 4.12| ). 



If this is large in any model we can probably reject that model because a large dipole moment in Hq is not observed. 
However, nearby it is difficult to measure Hq accurately due to peculiar motions and the discreteness of galaxies. There 
is a dipole moment in observations of more distant objects, which is assumed to be due to the fact that the Local Group 
is falling into the potential well produced by Virgo and the Great Attractor. The question is: what upper bound can 
be placed on the acceleration by observations? This issue will be discussed in § 5.2.5| . For now we simply use (5.27) to 
constrain the epoch of observation, T. 

Hubble Normalised Scalars 

If we wish to define, in analogy with FLRW models, Hubble normalised scalars, then we can do so as follows. We define 

n = ^ {2at + b)-\ (5.31) 



We can get a fairly good idea how these models behave, in comparison with what is know from local observations. We 
recall that 

0.3 <nQ<l, < r^An < 1 (5.33) 

for the present-day values of the density parameter and cosmological constant. For the Dabrowski models, we note that 
in the limit t ^ we have 

^=52' (5.34) 

We know that b is positive in order that the scale factor is positive, which implies that fl^t is negative at the big bang 
only if a > 0. However we also see that as t — > 0, then approaches some value larger than 1 if 6 < 1, as required by 
neglecting models with a finite density singularity. This means that if we desire SI S 1 today, then we force a > 0, in 
contradiction with the strong and dominant energy conditions (but not the weak one). 

As far as the quantity Tip' is concerned, we may 'equate' it with the cosmological constant term of the standard 



model. If we consider the expression for q, the deceleration parameter, given by equation (2.143) (which is essentially 
the Raychaudhuri equation) then the standard model has 

q^^n-QA (5.35) 

while for the D§,browski models we may write 

q^^n + ^Up (5.36) 
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neglecting the acceleration terms (which are zero at the center). We may therefore identify 

Qp ~ -2f7A. (5.37) 

Thus we may equate a positive cosmological constant with a negative pressure parameter. If we require this for all 
locations today (which means that R is quite small) then the sign of a doesn't matter: however, there will be a time 
(at t = —b/2a + -\/36^ — l/2^/2a) when the sign of Qp changes if a is negative. Of course we should take into account 
acceleration in the 'definition' of Q\, but this is only a qualitative discussion. 

What is interesting is that while the strong and dominant energy conditions favour a < 0, this seems to be ruled out 
by current measurements of the matter density of the universe. 

This may give cause for concern for these models, but it does not really concern us here. We will carry on with the 
analysis of models with a < throught this thesis, because we are more concerned with the Copernican principle in a 



non-standard cosmology, rather than finding the 'best-fit' IRF model (although, it will be shown in §5.4 that SNIa data 
favour a > too). This will be discussed later. 

5.2.2 The Age of the Universe. 

The original inspiration for D§,browski and Hendry (1998) to study Stephani models was the potential resolution of the 
age problem that they provided, which at that time seemed to be virtually insurmountable within the framework of 
FLRW models (even when a non-zero cosmological constant was invoked): the high measured value of Hq suggested 
an age, tq, of at most about 11 Gyr for a FLRW cosmology, whereas globular cluster ages were thought to be at least 
12-13 Gyr (Hendry and Tayler, 1996; Chaboyer et al. 1998). D^browski and Hendry (1998) and Paper I showed that, 
for the particular Stephani models they considered, this apparent paradox disappears: the D§,browski models have ages 
that are consistently 1-4 Gyr older than their FLRW counterparts (for an observer at the centre of symmetry, at least). 
However, the age problem has recently been alleviated by a recalibration of the RR Lyrae distance scale and globular 
cluster ages in the light of Hipparcos (Chaboyer et al. 1997), which has reduced the globular cluster ages considerably, 
to ~ 10 Gyr. The fit is still marginal, but the new ages are generally accepted as they allow a flat FLRW model to fit 
the observations provided that < 67 km s~^ Mpc~^. 

We certainly require, then, that at the epoch of observation our models are older than 10 Gyr. However, we will also 
consider the stronger constraint tq > 12 Gyr, partly to be conservative, but also because the diagrams for the 12 Gyr 
constraint are often clearer. Since it was shown in Paper I that best fit Dabrowski models are significantly older than 
their FLRW counterparts, we do not expect problems from this constraint. There is a surprise in store, though, when 
we consider non-central observers, and age will turn out to be the dominant constraint on these models. 

The age of the universe according to an observer at position ip and at coordinate time T is simply the proper time 
elapsed from the big bang (i = 0): 



To 



dt 



4|a|sin|Jl-62 cos^ ± 



where 




^ — 2al sm — 
2 



(for a 7^ and ip ^ 0: otherwise W = I and tq —T), where we take the value of T given by the solution of equation (5.26) 



as our constraint on the coordinate time for any specific iJg. In figure 5.3 we show proper age for observers At tp = n 
with Hq — 50 km s~^ Mpc~^ (which shows the behaviour of the age function most effectively) as a and b vary. The 
contour at tq = 12 Gyr is also drawn. As can be seen, below this curve the values of a and b can be rejected - the 
proper age is too low. 



In figure |5.4| the tq = 12 Gyr contours of the proper age function are plotted for Hq = 60 km s~^ Mpc~^ and for 
several observer positions, "0, showing how the age of an observer varies with for different parameters a and b. The 
shaded regions contain models that are less than 12 Gyr old for at least one of the observer positions. This demonstrates 
that the proper age of an observer is smallest at the antipodal centre of symmetry, ■0 = tt, so the contour at tq = 12 Gyr 
in figure marks the limit of viable models for this age constraint. Consequently, we will always use proper age 
at = TT to constrain the model parameters. We could weaken this constraint by requiring only that most observers 




Figure 5.4: to = 12 Gyr exclusion plot for different observer positions, for Hq — 60 km s ^ Mpc ^. In the shaded areas the 
models are not old enough. ^ = tt is clearly the most restrictive case. 
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are old enough, which would allow us to consider instead the age of observers at ip < 7r/2 while still satisfying the 
Copernican principle (half of the observers would lie in this region). For simplicity, though, we will not do this here. 
Finally, in figure |5.5| we show the age exclusion plot for the models (based on proper age at ■(/; = tt). We use three 




Figure 5.5: The age exclusion diagram for various Ho and proper age t — 10 Gyr. The shaded region represents the prohibited 
area. Also shown as dashed lines are the age limits for r = 12 Gyr. (Note that the region excluded for Hq = 70 km s~^ Mpc~^ 
contains the excluded regions for lower Ho - the progressively darker shading indicates this.) 

values of Hq: 50, 60, and 70 km s^^ Mpc^^, and a proper age of 10 Gyr (although the limits for 12 Gyr are also 
indicated). The shaded regions are excluded. It can be seen that unless we require the universe to be particularly old 
or the expansion rate high there is still a significant region of parameter space that cannot be excluded on the basis of 
age. 

It should be noted here that these plots are meaningless for 6 = 1, because then A = and the model is conformal 
to an FLRW model with flat spatial sections, for which x is not a good coordinate - as can easily be seen from (4.1 



5.2.3 Size and the Distance-Redshift Relation. 

When the spatial sections of a cosmological model are 3-spheres it may be possible for light rays to circle the entire 
universe, perhaps a large number of times. This is indeed the case for the models we are considering here. What will 
the signature of this be in the various distance-redshift relations? For the Dabrowski models this is a particularly easy 
question to answer because the fact that they are manifestly conformal to FLRW models allows us to determine the 
paths of light rays directly from the null geodesies of the underlying FLRW space. From this it is clear that light 
rays from a point directly opposite the observer (ie, from the antipode, x = will spread out around the universe 
isotropically from the antipode until they pass the 'equator' (x = 7r/2), where they will begin to converge and be focused 
onto the observer. As a result, a point source positioned exactly at the antipode will fill the entire sky when seen by the 
observer, so that its angular size distance, — (physical length/apparent diameter), is zero. Similarly, the refocusing 
of light onto a point produces an infinite flux at the observer, and therefore the luminosity distance, r^, is also zero 
(m ^ logio = — oo). Precisely the same argument applies to light that leaves the observer's position, travels through 
the antipode (where it will be focused to a point), and returns to the observer (being focused a second time). In fact, it is 
obvious that whenever the light rays travel through a parameter distance x that is an exact multiple of tt, va = ?'l = 0: 
this is reflected by the factor of sinx in (4.35). 
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Figure 5.6: Area distance from the centre as a function of x for two values of 6. T — 15 Gyr, a = — 1. 



This effect can be seen clearly in figures 5.(:-5.8, where we show the two principle measures of distance as they vary 
both with coordinate distance x ^nd redshift. Viewed as a function of x, in figure the zeros of the angular diameter 
distance clearly occur at multiples of tt for all model parameters. Looked at in terms of redshift, though (figure 5.7), 
it is clear that for small b the zeros are much closer together than for larger &, with the first zero occurring at z ~ 1 
for b = 0.25. Figure 5.8 shows the luminosity distance-redshift relation, for comparison. These effects are not as unusual 
as they look, and can be found also in FLRW geometries for models with positive A - see §4.6.1 in Ellis (1998) and 
references therein. It is noteworthy that a blackbody situated at the opposite pole to the observer would look exactly 
like the CMB if it was at the right redshift. For example, a star of surface temperature 3000K at a redshift z = 1000 
would have an apparent temperature of 3K, and would cover the whole sky. 

Can we rule out such apparently aberrant behaviour? Theories of structure formation are fairly well developed (see 
Bertschinger et al. 1997 for a thorough discussion and references), and the evolution of galaxies and the star formation 
rate (SFR), while not accurately known, are at least qualitatively understood. In particular, the SFR, which is very 
important for determining the luminosity of distant, young galaxies, is believed to fall off beyond z ~ 3 (see, for example, 
Loeb, 1999). As a result, one could argue that there will be relatively few bright objects beyond some redshift zsf that 
corresponds to the epoch at which galaxies 'turned on' and the SFR began to increase significantly. This would mean 
that the zeros in the distance-redshift relations would be essentially unobservable if they occurred at redshifts larger 
than ZSF because there would be no luminous objects to be seen magnified in the sky, whereas if the zeros occurred at 
lower redshifts than zsf one could reasonably argue that there ought to be some signature of this in the observations. 
Since galaxies have only been observed (in the Hubble Deep Field, for example) with redshifts up to z ~ 5 we tak e 
Zsf = 5 (although the constraints imposed by larger zsf can easily be inferred from the exclusion diagram, figure 5.9). 
If these arguments are not completely convincing, then, at a simpler level, the fact that the observed magnitude-redshift 
relation is known accurately out to z ~ 1 from type la supernovae (Perlmutter 1999), and is certainly not dipping down, 
allows us to say that there is no zero of luminosity distance below 2 ~ 2, say. We therefore also consider the constraint 
that results from requiring that there are no zeros below z — 2. 

We wish, then, to constrain the parameters of our models by rejecting any models for which the first zero in the 
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Figure 5.7: Area distance from the centre as a function of redshift for the same parameters as figure 5.6. For small b the angular 
size distance oscillates far too rapidly. 



b=0.7 




Figure 5.8: Luminosity distance from the centre as a function of redshift for the same parameters as figures 5.6 and 5.7 
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distance-redshift relations occurs at z < z^, where = 2 or = zg/;- = 5. Using ( 4.29| ), this means 

^ RoW{^,TT,t^) 



'Ax 



R(U) 



> 1 + Zjr, 



(5.39) 



where Rq = R(T), Wo = W{^,T) (giving the conformal factor at the observer) and t,r denotes the l ookb ack time ( 5.12 ) 
at y = TT . Again we determine the epoch of observation (ie, the observer's coordinate time T) using (|5.26| ). The solution 
of (5.39) for a and b is shown in figure 5.9 as an exclusion diagram. The effect of this constraint is to rule out small 




Figure 5.9: Exclusion plot obtained by requiring that the first zeros of rA{z) occur at z > Zt^, for Ztt — 2 and = 5. The shaded 
regions are excluded. Curves are given for Ho ~ 50 and Ho — 70 km s"'^ Mpc~^. 



values of b, for any a. This is a reflection of the fact that, loosely speaking, b measures the 'size' of the universe: for 
small times the scale factor goes as bt, so that when b is small the spatial sections are small, light rays don't take long 
to travel from antipode to observer, the scale factor changes relatively little during this time and the redshift of the 
antipode (which is dominated by i?o/-R(^7r) as in FLRW models) is small. 

As a coda to this section we consider the effect of de manding that the first zero of ta is effectively unobservable as 
a result of being 'hidden' behind the CMB. Figure 5.1C| shows how the redshift of the first zero of r^(z) varies with b. 
If, instead of choosing z^r = zsf as our primary constraint, we want the first zero of r^(z) to happen at a redshift large 
enough for the universe to be opaque (ie, before decoupling), then figure 5.10 shows that b must be quite close to unity. 
This figure also allows the extent to which values of b are excluded for any z^r to be estimated. 



5.2.4 The Microwave Background Anisotropy. 

The CMB is observed today to be a blackbody at a temperature of Tq = 2.734 ± O.OIK, with a dipole moment of 
Ti = 3.343 ± 0.016 x lO'^K and quadrupole moment as big as Ta = 2.8 x lO^^K (see Partridge 1997 for details 
and references). It was emitted at a time when the radiation was no longer hot enough to keep hydrogen ionised, 
causing it to decouple from matter, which happens at Tdec ^ 3000K. Idealised cosmological models do not have realistic 
thermodynamics (that is, they do not, in general, describe the thermodynamic evolution of the gas and radiation 
mixture that fills the real universe) . In FLRW models the epoch at which decoupling occurs is simply defined to be that 
corresponding to the redshift necessary to shift the temperature at decoupling to the observed mean temperature of the 
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Figure 5.10: Logarithmic plot of the redshift at which the first zero of rA occurs as a function of b for Ho = 50 km Mpc~^ 
and a = —1. For the first zero to occur at z > 1000, so that recombination occurs 'nearer to us' than the antipode - ie, at x < ^ ~ 
requires 6 ~ 1. 



CMB, Tq. From the redshift relation applied to the temperature of a blackbody (T will be used to denote temperature 
in this section), 

Tobs = (5.40) 

we infer that the CMB is formed at a redshift z « 1000. This definition is fine for homogeneous models, leading to a 
consistent definition of the time of decoupling for every observer at the same cosmic time, but raises an interesting point 
for the inhomogeneous D§,browski models, because the redshift depends on both the obse rver's position, -0, and the 
angle around the sky, 9. If we simply define the redshift of the CMB at any point to satisfy ( 5.40| ) with Tobs = To then. 



by definition, we obtain a perfectly isotropic CMB for that observer, but we must choose a different emitting surface 
for each different observer. Such an observer-based definition of the CMB surface is clearly unsatisfactory. Instead 
we propose several alternative ways to define what is meant by the CMB surface (it is important in these definitions 
to distinguish between the dominant strange matter that is responsible for the geometry of the D^browski models - 
see § ^.l.l| - and the putative 'real' gas that decouples): 

1. On a surface of constant cosmic time, t = Icmb - since the D§,browski models possess a cosmic time coordinate 
with respect to which only the pressure is inhomogeneous this is a natural extension of the FLRW definition; 

2. In general inhomogeneous models we could choose hypersurfaces of constant density, which at least has some 
physical basis ~ for the D§,browski models this is equivalent to [|. because the density is homogeneous on cosmic 
time surfaces; 

3. On surfaces such that p/ p ^constant (which gives constant temperature for an ideal gas) - for our models this 
would mean that decoupling occurs at different times at different places in the universe (probably a good thing in 
an inhomogeneous universe); however, the Dsibrowski matter is not an ideal gas so the validity of this definition 
for these models is debatable; 

4. On a surface of constant proper time (based on the assumption of some common evolution for the ideal gas 
component at different positions) - see equation ( 5.38| ); 



Finally, we could avoid all consideration of the physics of decoupling and simply assume that it happens at such 
an early time that we can effectively define the CMB to be free-streaming radiation 'emitted at the big bang', as 
is implicitly assumed in the EGS theorem (Ehlers, Geren and Sachs 1968). This only really makes sense if there 
is some natural definition of the radiation field at early times. For example, if the model is homogeneous and 
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isotropic at early times, we can define a homogeneous and isotropic radiation field. Our models have exactly this 
property of homogeneity at early times (as can be seen from ( ^.8| ), if i ^ then W ^ 1). 

Although I) is interesting we will not consider it here because it does not reflect the physics of the CMB, and, on a more 
practical level, it prevents us from constraining the CMB anisotropy, because there isn't one. Also, the homogeneity of 
D§,browski models at early times means that for small t the proper age is virtually identical to the coordinate time t 
(equation ( 5.38 ) with W ^ 1). It turns out that for times of observation that reproduce the observed Hq any reasonable 
definition of the CMB surface puts it at an early time, which means definition ^) is virtually identical to |^) We will 
therefore define the CMB according to |]) in this section. 

It still remains, though, to decide exactly which surface of constant cosmic time the CMB originates from. Consider, 
for an observer at position ^, the temperature distribution on the sky that the CMB would have if it were emitted from 
the surface t = tcMs (related by the lookback time formula (5.12) to some distance Xcmb)' 



Tdc 



^ + Z{'4',XCMB,0)' 



(5.41) 



Equation (5.22) shows that we can write 



where 



1 + z{ip,XcA 



1 + 2o(V', Xcmb) = 



= 1 + Zq{'>P,Xcmb) + Zi{lp, Xcmb) cost 
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Rc 



2a 



(1 - cos V' cos Xcmb) 



. , ^ 2a Rq . u{ijj) Rq 

zi[iP,Xcmb) = TTTir^^'^V sm Xcmb = — ~T^777 ^^'^ ^^mi 

cA Wo V AWo 



(5.42) 



(using (4.12) in the last equality and assuming a < 0). The mean redshift of the CMB surface is zq; zi gives rise to an 
anisotropy in the CMB. We can therefore define the location of the CMB surface to be the Icmb (or Xcmb) that gives 
a mean redshift of 1000. That is, Xcmb is the solution of 



1000 



(5.43) 



for any observer position ip. 

Having found Xcmb we can evaluate the anisotropy in the temperature of the CMB. Since Tobs depends on the 
reciprocal of 1 + z the dipole moment in z will give rise to higher multipoles when expanded as a binomial series: 
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that is. 
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The dipole moment of the CMB temperature is then (using ( 5.43D ) 

^1 - -,^-3 



cos^ 9 + 0(cos^ 



<5i 



1 + zo 



IQ-^zi. 



(5.44) 



(5.45) 



(5.46) 



Measurements of the CMB can now be used to constrain the model parameters. We at least require that the dipole 
moment should be no larger than the observed dipole anisotropy, \5i\ < Ti/Tq « 10~^ (ie, \zi\ < 1). If this is satisfied 
for any model then it is clear that the quadrupole and higher multipole moments will all be S 10~^ - certainly no larger 
than their observed values. In fact, such a constraint on zi is very weak, leaving vast tracts of parameter space entirely 
untouched. Moreover, there are very good reasons for believing that there is a significant contribution to the observed 
dipole moment from the peculiar velocity of the Local Group as a result of infall towards the Great Attractor (Lynden- 
Bell et al. 1988), which can be measured with moderate accuracy using galaxy surveys, and seems to be consistent with 
the motion of the Local Group with respect to the CMB (Schmoldt et al. 1999). In an inhomogeneous cosmological 
model, though, one expects local anisotropics that could be interpreted as bulk flows, and it is not beyond the bounds 
of possibility that all of the CMB dipole and the local anisotropics could be explained purely as cosmological effects in. 
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say, a Stephani universe without any need to invoke peculiar motions and local inhomogeneities (ie, perturbations of 
the perfect background cosmology). Such perturbations must exist, of course, and will make some contribution to the 
observed anisotropics, but it has not been shown that they are the only, or even the dominant, contributions. However, 
we choose here to reject iconoclasm in favour of the more conservative viewpoint that most of the observed dipole is 
due to the peculiar motions induced by local inhomogeneities, but that there remains some leeway - up to 10% of the 
observed dipole - due to observational uncertainties, for there to be a purely cosmological contribution to the CMB 
dipole. Then the largest dipole moment that we can accept from our models is \Si \ < 10^"*, or 



|zi| <0.1. 



(5.47) 



Given any model parameters and some observer position we adopt the following procedure. First we use ( 5.26| ) to 
determine the epoch of observation for some Hq, as usual, then we solve for Xcmb using ( ^.43 ). Having found all the 
parameters we need to determine zi we simply che ck (5.47 ) to see whether the model, or at least that observer position. 



must be rejected. In practice we can simply solve ( 5.47) as a n equality to obtain a as a function of b at the boundary of 
the allowed region, and this is what is shown in figure 5.11 (for tJj = 7r/2), where it can be seen that a low value of Hq 



— H=50 
-- H=70 




Figure 5.11: The exclusion diagram from CMB anisotropies for an observer at ip — -k 12. The curves for Ho = 50, 70 km s'^ Mpc" 
are shown. As Hq, increases, the 'fingers' move down to more negative a. The excluded region lies within the fingers. 



constrains our models most - in contrast to the age constraint. This is because Hq decreases monotonically with time, 
so small i?o corresponds to a later time of observation and therefore a later time for the CMB surface, which means 
that W has evolved to become more inhomogeneous. We choose ij) — ti/2 because, as is clear from (3.9) and ( 5.42| ), 
the anisotropy is generally worst there, so if a model is rejected at = 7r/2 it will be unacceptable everywhere, in 
order to satisfy the Copernican principle. Again, therefore, for models not excluded in figure ^.11 detection of a CMB 
anisotropy of the magnitude that we observe would be typical, and the Copernican principle need not be abandoned for 
these models to be viable. 

takes on different 



The finger-like excluded regions in figure 5.11 appear because for different model parameters Xc 



values, and for some parameters this value is very close to tt, so that the CMB surface is almost at the antipode of the 
observer. This means that the entire CMB is emitted from virtually a single point. Since redshift depends only on the 
relative conformal factors at emitter and observer for our conformally flat models, the CMB must be almost exactly 
uniform however inhomogeneous the model (ie, whatever the value of a). 
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5.2.5 The Local Dipole Anisotropy 



Although wc arc not in a position to use real observations to constrain the dipoles that would be detected in observations 
of the 'local' universe (in galaxy surveys, for example, where z S 0.01, or with type la supernova data, for which z S 1), 
we can at least consi der th ese effects qualitatively. 

From ( 5.22|) and ( 4.12| ) the redshift dipole for objects at any radius x from an observer at tp is seen to be 



u{ip) Ro 



AWo 



smx 



{cf. equation ( ^.42 )). It is clear that any constraint on the anisotropy in the properties of objects at some distance from 
the observer will amount to a constraint on the acceleration, zi, of the fundamental observers, and vice versa. In the 
same way, constraints on the acceleration place limits on the possible dipole moment in the measured if™ in ( 5.30| ). 
Defining 5H = u/c (according to ( 5.30| )), the dipole in redshift becomes 



i?o sinx 5H 
Al^n c 



(5.48) 



If we assume that the time of observation, T, is fairly close to t = (as is generally the case for the values of Hq we 
allow), then Wq ~ 1 and aT ^ &, which means that Rq/c « hT and Rq/c k, b. Equation (5.48) then becomes, with the 



help of (5.26) 



1 Rq sin X 
cHo 



SH'. 



SH 



(5.49) 



Note that dependence on distance from the observer only arises through the sinx factor, so for observations of objects 
at any distance from the observer 

1.1.-^^^. (.50) 

Local measurements in principle, therefore, constrain the acceleration and the anisotropy at all redshifts. For example, 
applying this to the CMB dipole, ( |5.50| ) shows that the constraint ( 5.42 ) imposed in § 5.2.4 will always be satisfied 
provided that 



< 0.1 



A 



(5.51) 



To ensure that the variation in Hq is less than 20% requires VA/b < 2, or b > l/\/5 ~ 0.45. For a variation of less 
than 10%, we must have b > l/V^ ~ 0.71. Since we are interested in all values of b > 0.5 (see figure 5.12| ) we must be 
prepared to countenance variations in Hq around the sky of up to 20% if the local anisotropy is not to impose tighter 
constraints on the model parameters than those already derived from the CMB. Note, though, that ( 5.51 ) is considerably 
stronger than is really required for most model parameters, owing to the fact that the sin^ factor in (5.49) was neglected. 
This amounts to adopting as the CMB constraint the envelope of the fingers in figures 5.11 , 5.12 and ^.13 , which would 
obviously overlook large areas of parameter space that should really be allowed. 

Although it would seem that the easiest way to constrain the acceleration in the Dabrowski models would be to 
measure the local Hq dipole, it is not really possible to measure SH/Hq accurately, because that would require accurate 
measurements of galaxies in different directions at very small redshifts; there are relatively few such galaxies, and those 
that there are possess significant random peculiar motions that make a large contribution to the errors on any estimate 
of SH/Hq. Moreover, there is known to be a dipole in observations of galaxies at somewhat larger redshifts (z ~ 0.01), 
which is usually interpreted as the effect of infall of the Local Group towards the Great Attractor (Lynden-Bell et al. 1988; 
Schmoldt et al. 1999). This infall manifests itself as a systematic relative motion of the Local Group with respect to 
distant galaxies, at a velocity v ~ 600 km s~^, corresponding exactly to the Local Group motion relative to the CMB 
frame. In order not to confiict with these observations we at least require that at redshifts zq ~ 0.01 the dipole moment 
due to the Dabrowski acceleration is no larger than the observed dipole: 



CZl 



= 600 km s" 



Using (4.30) and (4.35) it is possible to rewrite (5.48) in the form 



SH 

X) = (1 + za)rA{^, x) , 
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where we implicitly identify zg and as the mean redshift and angular size distance of objects at coordinate distance %. 
Then, since at low redshift the Hubble law is valid, czq = v = i?o''prop (^prop denotes proper distance), and = = 



3, we have 



Zl 



That is, 



CZi 

czo 



czq SH 

600 km s^ 
0.01c 



= 0.2: 



(5.52) 



at most a 20% variation in Hq around the sky, consistent with the CMB result ( 5.51 ). 

Finally, we consider the magnitude-redshift relation in the low-redshift limit and apply it to type la supernovae (for 
example, to the relatively low-redshift supernovae of Hamuy et al. 1996 - although it is not unreasonable to assume that 
this gives at least an order-of-magnitude estimate of the size of the anisotropy for observations at higher z characteristic 
of the supernova data of Perlmutter et al. 1999). If 2 ^ 1, then = czq/Ho- Equations (4.37) and (5.30) give 
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(assuming that SH/Hq is small). That is. 



Sm = 



5 SH 



InlOi^o- 

The dispersion on the magnitudes of type la supernovae (due to both observational errors and intrinsic dispersion) is 
estimated to be Sm ~ 0.3 (Hamuy et al. 1996; Perlmutter et al. 1999). For the dipole resulting from the acceleration to 
be undetectable via type la supernovae we must have 



SH_ 
Ho 



0.14. 



(5.54) 



This is somewhat stronger than the other local constraints derived above, and would suggest that the su pern ova data 
will indeed restrict the model parameters more than the CMB anisotropy. Comparison with ( 5.51 ) and ( ^.52 ) su ggest 
that b ma y only be half or a quarter as large as would be permitted by the CMB constraint (see equations ( |5.30[) 
and (4.12)). However, these are only very approximate results, and require more detailed analysis. In particular, 
there is observed to be a dipole in the supernova data (explained in exactly the same way as the other dipoles in the 
conventional interpretation), and this means that the constraint (5.54) must be re-evaluated: the real constraint is likely 
to be somewhat weaker than (5.54). 

It should be borne in mind that throughout the preceding discussion we have only considered the modulus of the 
dipole, not its direction. It is clear from ( 5.42 ) that the local dipoles may be in the same direction as the CMB dipole 
or in the opposite direction, depending on the sign of sin Xcmb- What is more, the variation of the dipole with distance 
is controlled entirely by sinx, so the dipole will change sign whenever x is a multiple of tt. 

In this section we have examined the local anisotropics resulting from acceleration in a rough and qualitative way. 
Although certainly not conclusive, the results indicate that it is important to give full consideration to the constraints 
imposed by local observations on these anisotropics. 



5.2.6 The Combined Exclusion Diagrams. 



When we combine all of the constraints derived in this section (figures 5.12 and 5.13) we can see that for Hq 
50 km s^i Mpc~ 



^ the strongest constraint comes from the CMB, with age placing somewhat weaker limits on the 

The 'size' restriction of S 5.2.3 



5.3) 



allowable degree of inhomogeneity (which is measured largely by the size of a - see 
eliminates quite a large region of parameter space for small 6, but this is not really a constraint on the inhomogeneity, 
which is our principle concern. 

Perhaps rather surprisingly, given the results of D§,browski and Hendry (1998) and Paper I, the strongest constraint 
for larger Hq really comes from the age. As can be seen in figure 5.13| , the exclusion plot for Hq = 70 km s~^ Mpc~^, the 
CMB constraint pokes out in places to eliminate certain regions, and the size constraint cuts off low values of b, but age 
does most of the dirty work. It can also be seen that for Hq — 80 km s~^ Mpc""'^ age imposes a very strong constraint 
on the models (dashed line in figure ^.13 ): the models must be very nearly homogeneous. However, if we relaxed the 
age constraint to (a probably quite acceptable) 10 Gyr the CMB anisotropy would be the dominant limitation for most 
values of Hq in the currently fashionable range (50 S Ho :£ 80 km s^^ Mpc^^). 
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Figure 5.12: The complete exclusion diagram for all the observational constraints studied (age, size and the CMB anisotropy), 
for Hq = 50 km s~^ Mpc~^. We have taken the 12 Gyr age constraint, to be conse rvative. The dominant energy condition should 
be added to these constraints: it eliminates models with h > 0.82 (equation (p.l9|)). 




Figure 5.13: As in figure 5.12, but for Ho — 70 km s ^ Mpc ^. The age constraint for Ho — 80 km s ^ Mpc ^ is also shown as 
a dashed line (close to the 6-axis) ; high Ho means a very low age, just as for the FLRW models. 



We should not forget, at this point, to reintroduce the restriction (5T_£) from the dominant energy condition, which 
rules out high b. This is not shown on the diagrams, in order to avoid clutter. Most of the models eliminated by 
this constraint have already been ruled out by the age or CMB constraints, and models that are rejected solely by the 
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dominant energy condition are not hugely inhomogeneous (see the next section). 

5.3 The Size of the Inhomogeneity. 

While we have considered many different aspects of the D;ibrowski models, what we have not done is to assess the 
extent to which the models that are not excluded are inhomogeneous. It is obvious from the exclusion plots, figures 
and 5.12 , that the homogeneous Djibrowski models (those with a — 0) are the 'most acceptable', in that all the constraints 
favour small a. This should not be surprising, as far as anisotropy constraints are concerned, at least. What is not clear 
is whether the allowed region only contains models that are very nearly homogeneous. We will show that it does not. 
The most natural way to assess the degree of inhomogeneity of the models is to examine the variation of the pressure 



over surfaces of constant cosmic time. It can be seen from (5.15) that the extremes of pressure occur at x and x = ""i 
so we define the inhomogeneity factor 11 to be the relative pressure difference between the two poles: 
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(5.55) 



If n > 1 then it is reasonable to say that the models are truly inhomogeneous, whereas if 11 <C 1 they are obviously 
nearly FLRW. Note, though, that 11 depends on the cosmic time surface under consideration: for small t, 11 w 0, and 11 
reaches its maximum at t = —b/2a (see (^.l[)), at which time 



n = 



2P 



1-62 



= 2- 



1- A 



so that the models are significantly inhomogeneous (11 S 1) when 



b> 



V3 



0.58. 



(5.56) 



Most of the allowed models in figures 5.12 and 5.13 are unaffected by ( 5.56D - models with smaller b have already been 



eliminated by the size constraint in §5.2.3 



This is not really a fair reflection of the inhomogeneity of the models at the times of observation that are relevant 
here, though, because the Hq constraint ( 5.26D generally ensures that the epoch of observation is quite early on in the 
evolution of the universe when the scale factor is somewhat smaller than its maximum size. To evaluate the impact 
of this, consider two specific examples. From the allowed region of (5.12) choose the model at a = —7, b — 0.75. 
For Hq = 50 km s"^ Mpc"^ the solution of ( ^.26[ ) is T « 0.016 Mpc s km"\ which gives 11 = 1.33. For the model at 
a — —8, b — 0.64 we get T = (3 — \/5)/50 0.015 Mpc s km"-'^ and 11 = 0.86, which is close enough to 1. These models 
are certainly not 'close to FLRW'. 

Nevertheless, it could be said that the models are not massively inhomogeneous, and that their degree of inhomogene- 
ity only reflects the looseness of the constraints applied. This is not so. Firstly, we have at every stage chosen stronger 
limits than were strictly necessary, especially as far as the age is concerned, where we could have adopted a 10 Gyr limit, 
rather than 12 Gyr, which would have increased the allowed region considerably. Most importantly of all, though, is 
the fact that even for models that are only inhomogeneous at the 10% level (11 — 0.1), the CMB anisotropy, 5i < 10~^, 
is at least three orders of magnitude smaller than the inhomogeneity it permits. This certainly conflicts with the spirit 
of the almost EGS theorem of Stoeger, Maartens, and Ellis (1995) (although not the actuality, since our models have 
acceleration, whereas the theorem deals with geodesic observers), which says that small CMB anisotropics give rise to 
small perturbations from homogeneity. 



5.4 Fitting to FLRW Models 

So far we have demonstrated that the D§,browski models do not have any problem with regard to the 'global' constraints 
we have considered. As far as 'local' constraints are concerned we have only considered the size of the dipole, and 
ensured it was not too large. We still have to show that the Dabrowski models are capable of providing an acceptable 



fit to the observed magnitude-redshift relation. A glance at figure 6.14 (with reference to figure p.4| for the FLRW case) 
shows that these IRF models are quite capable of producing very incompatable m — z relations. In order to demonstrate 
that the D§,browski models are acceptable in this regard we 'match up' the D§,browski m — z relation to the FLRW 
one. Because of the local dipole considered above we will simply fit the two curves for an observer located at the center. 
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(Obviously a best fit from any location will pick out the center as the location most able to fit the FLRW models.) This 
should provide confidence, if not conclusive evidence, that the Dabrowski models are acceptable in terms of this test. 

The complete solution of fitting the two m — z relations is quite complicated. Given an actual data set it is very 
easy to decide whether either model is consistent with it: form a statistic and check whether it is small enough for 
the model to be consistent at the desired confidence level. To distinguish two theoretical models, however, it is, strictly 
speaking, necessary to adopt a sort of 'two-level' approach and ask: assuming that the Stephani model is the true model 
of the universe, what is the probability, pi, that we would, for any realistic data set, reject the FLRW model at some 
confidence level, p2? To answer this rigorously requires consideration of all possible data sets drawn from both the 
FLRW and the Stephani distributions. In fact, we will use a much simpler and more intuitive resolution criterion based 
on the mean-square difference between the theoretical models. This will be more than adequate for our purposes. For 
brevity we give a detailed derivation for the m — z relations only, any other case being similar. 

The mean-square difference between ms{z) and mp{z) over some z-range < z < Zmax, for which we expect to have 
observations, is 

D2 = / [ms{z) - m,F{z)f dz (5.57) 



^max J 

We assume (as is justifiable for the case we consider here) that the difference ms{z) — mp{z) is smooth and slowly 
varying, and also that the integral does not need weighted - ie, that the number density is uniform in redshift. Then 
we imagine that the difference is just a constant; ms{z) — mi?(z) — \fD2- We can accept or reject this difference on the 
basis of the usual analysis. 

We adopt the null hypothesis that the FLRW model is correct so that the data set consists of n data points {wi, Zi\ 
such that the errors on the mi are normally distributed with variance cP' say; ie, mi — rapizi) -I- ti where ~ A^(0, cr^). 
Then, clearly 

ms(z)-mF(z) = /d^< (5.58) 

(since ms{z) — mi?(z) =const. implies the best fit value of ms{z) — mp{z) is (m^ — mp{z)) — (e^) and the variance on 
this estimate goes as jn). 

Now, for the SNIa of Perlmutter et al. (1999) n — A2 and their error on the intrinsic magnitude dispersion is cr « 0.15. 
This implies 



2 



Da < — w 5 X 10 



-4 



for these data. We take z^^x — 1- When we fit the Stephani to the FLRW model we require Z?2 ~ 10"'' to demonstrate 
that the Stephani models are capable of providing at least as good a fit to the SNIa data as an FLRW model. 

The procedure we use is to fix iJg and Og and find the 'best fit' parameters a, h and T for a range of (effectively 



JIa). For simplicity we consider two values of fig = 0.3, 1.0. In figures 5.14 to 5.16 we show the best fit values of a, h and 
T against go • We can see that a negative go requires a > as expected from the strong energy condition. In figure |3.17| 
we show the goodness of fit: provided D2 < 10"** then these models will fit the SNIa data acceptably. Finally, in 



figure 5.18 we demonstrate that a positive value of a will enable an old enough universe. 
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Figure 5.17: The mean-squared difference as a function of go- Provided this function satisfies -D2 < 10 * then the SNIa data will 
not be able to distinguish between the FLRW and Stephani models. This is clearly satisfied in these cases. 
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Figure 5.18: The best fit values of a vs. T. Tliis clearly shows that a positive value of a will provide a sufficiently old universe. 



5.5 Conclusions. 



We have discussed the observational relations and other properties of a family of inhomogeneous perfect fluid cosmological 
models which, in contrast to FLRW models, have pressure gradients and consequently acceleration of the fundamental 
observers. It was demonstrated that these models do not suffer from particle horizons (although they also do not contain 
a radiation field, which, were it to be present, would dominate at early times, leading to a more rapid evolution of the 
scale factor at early times and the reintroduction of horizons) . More importantly, our studies have shown that there is 
a significant subset of this family that arc markedly inhomogeneous but cannot be excluded on the basis of t he te sts 
considered here. It is possible, for every observer in each of the models in the allowed regions of figures ^.12 and 5.13 , to 
choose the epoch of observation so that the observed value of Hq is reproduced, the age is greater than the measured age 
of the universe and there are no obviously unacceptable features at low-redshift {z S 5) in the observational relations. 
It has also been shown that it is possible for these models to reproduce the FLRW magnitude-redshift relation to high 
accuracy. Most importantly, though, the dipole in the cosmic background radiation would be considerably smaller than 
the observed CMB dipole: despite the inhomogeneity of the models the anisotropy they produce is very small. The fact 
that this is true for every observer means that it is not possible to reject these models by appealing to the Copernican 
principle. As a result, the standard assumption that the observed high degree of isotropy about us combined with 
the Copernican principle necessarily forces the universe to be homogeneous (the cosmological principle) is seriously 
undermined. 

It is clear from 5.12] and 5.13 that the constraints on the inhomogeneity are much more severe for values of Hq at the 



upper end of the currently accepted range {Hq Z 70 km s^^ Mpc^^), and from the results of §5.2.2 it can be see that a 
higher age also results in stricter limits. It should be noted, though, that the constraints we have adopted are invariably 
stronger than is strictly required by observations, so that the results are conservative. Moreover, the principal interest in 
this thesis is the relationship between anisotropy and inhomogeneity, particularly with regard to the CMB. The fact that 
age places a strong constraint on the models for large Hq is less interesting than the fact that the anisotropy induced in 
the CMB is smaller for larger Hq: it seems quite reasonable to believe that it should be possible to find models that are 
similarly isotropic, though inhomogeneous, but that have much more efficacious age characteristics. In fact when a > 0, 
which we require for qo according to the best fits, are much older than their a < counterparts. 

In the following chapter we investigate the properties of the more general C 7^ models. We will demonstrate that a 
substantial deviation from spherical symmetry does not cause major problems for the 'global' tests we have considered 
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so far. However, the local dipole is large enough to possibly preclude some of these models. 



Chapter 6 

Observational Characteristics When C 



In chapter |5| we looked into the general observational characteristics of the spherically symmetric Stephani solutions 
{( — 0) which admit an isotropic radiation field, with the scale factor R restricted to be a quadratic, and we demonstrated 
that the models have no significant problems with fitting the broader aspects of our universe. It was also shown that 
some of the allowed models are significantly inhomogeneous. There is no particular reason for choosing a scale factor 
of this form, and it is highly likely that the freedom in R{t) could mean that there are even better forms for the scale 
factor to fit observations. This is a degree of freedom which we do not exploit here. 

This chapter is devoted to giving the results for the more general case of C 7^ (the 'isotropic radiation field' or IRF 
models). The idea is exactly the same as before, but the algebra is a bit more messy. In fact none of the observational 
relations can be given in analytic form, so graphical methods are used to show that there is a distinctly non-zero 'volume' 
of parameter space that cannot be rejected by current observations. 



6.1 Hyperbolic Stephani Models: Overview 



We have commented before on the conformal nature of the IRF models, which have allowed us to perform rotations to 
arbitrary locations in the spacetime. The same conformal nature has allowed us to find all the relevant observational 
relations for these models. The inhomogeneity of the models is 'contained' in the conformal factor W. In fact, because 
the light rays are unaffected by the conformal factor, we see inhomogeneities as they were in the past, which allows the 
CMB to be so isotropic (the universe is homoge neou s at early times. 
We have the conformal factor of the metric (4.23) 
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We have k restricted by (|4.3|); 
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where 



S± sign ( -A 



2a 
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(6.4) 



We take the sign in front of the square root term to be the sign of ^A — 2aR/c in ( [l.3| ) to give the Dabrowski models 
in the limit C whatever the sign of ^A — 2aR/c. We take R to be 



R{t) = ct{at + b). 



(6.5) 
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Note that A M 1 — 6^ and we replaced 6 — a/c. 

The units we wih use are as foUows: [c] = km s^^, Xj'V'jC ^-^e dimensionless, R is in Mpc and [t] — Mpc s kni^^ = 
[1/i/o], so that [a] ~ km s~^ Mpc^^= [Hq] and b is dimensionless. Note that we have redefined C '-^ C/c to give the 
units of [C] = [Hq]. 

We will now use all the constraints of chapter 0. 



6.2 Energy Conditions 



As before, we don't have a specific form of thermodynamic scheme for the IRF models, although chapter ^ proved the 
existence of one. Knowing the specific form of such a scheme would be desirable, but rather difficult (Sussman, 1999). 
Rather we wish to limit the models just using the energy conditions: 



Ai > 
f p > 
/ic"" + 3p > 
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\v\ > 0. 



(6.6) 
(6.7) 
(6.8) 
(6.9) 



The first must always be satisfied; the second is the weak energy condition, while the other two are the strong and 
dominant energy conditions respectively. They say that for all reasonable forms of matter the energy density is positive 
and the pressure cannot be very large and negative; the dominant conditions goes even further and says that positive 
pressures can't be large either. For the IRF models we have 
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(6.10) 



The conditions are pretty complicated and will not submit to an analytic treatment. We can turn to graphical methods 
to give a rough guide. It turns out that, as in the Dgibrow ski ca se the strong condition requires a < 0, and the dominant 
energy condition fails at certain times and places. As in § |5.1.l| , things are worst when the universe is half way through 
its evolution, t = — &/2a; similarly on one timelike slice the condition will hold in some places (near ■0 = 0) and fail at 
others (ip = n). We test the condition at this spacetime point {t — — 6/2a, tp = tt), on the assumption that if the models 
are unacceptable there then they will be unacceptable everywhere. 

In figure 6.1, we show that a high value of ^ combined with a small value for a fails the dominant energy condition. 
As a becomes large and negative, C becomes less and less important, with b being most important for determining the 
acceptability of the models. As in the case ^ — 0, any value of 6 > ^2/3 is always disallowed. 

To try to get a feeling for this, consider figure 6.2, where we show the same as figure 6.1, but now in the a — 6 plane. 
As C becomes large, a must be large and negative or the dominant energy condition will be broken. 

The Stephani models have enough problems with their matter content without chastising them for failing this 
condition. We will banish this problem from our minds, in the knowledge that it actually won't really matter because, 
as we will see below, the dominant energy condition doesn't rule out anything that isn't more or less ruled out by other 
conditions. The matter is 'reasonable enough' to pass both the weak and strong conditions. 



6.3 Ho 

We must have our observers living at a coordinate time which is representative of us. We do this by ensuring that the 
observers measure a (mean) Hubble's constant (ie, 9/3) the same as we measure. It will not change over a surface of 
constant time. 

As before, we require the expansion rate to lie somewhere in the range 50 < Hq < 80 km s^^ Mpc^^. This means 
that for any given choice of parameters the age, T, of the universe is given by the solution of 
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For simplicity, in this chapter we take a fixed value of Hq — 60 km s^^ Mpc~^, because the added complication of 
considering a range of Hq does not warrant the effort: in chapter ^ the differences were only small. So, for any given 
values of the parameters C, a, b we have a value of the coordinate age T given by the solution of 

R / K\ 
"A \R),t 

which is pretty complicated (it's a sixth order polynomial in T). 



= 60, 



(6.12) 



6.4 Age 

It is known that the universe is older than 10 Gyr old, and probably older than 12 Gyr (cf, i ^5.2.2| ). The standard model 
has gone through various crises over this problem, but at the moment things seem to be acceptable, provided that Hq S 
67 km Mpc~^. We require that any observer in the spacetime, with a local expansion rate Hq = 60 km s~^ Mpc~^ 
be older than 10 or 12 Gyr. The proper age of any observer is given by 



where T is the solution of ( 3.12| ). This integral may have an analytical solution, but we just use a simple numerical 



analysis, as for all the other tests. To gauge what is going on, in figure 6.3, we show a surface plot of proper age verses 
b and C for a = —5. contours at 10 and 12 Gyr are marked on, demonstrating that for high the age is far too low, for 
the expansion rate we measure. This means that if we introduce a large C the expansion rate will be too low when the 



universe is old enough to accommodate globular clusters. This can also be seen in figures |6.5| and 
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Figure 6.2: Regions of the a — b plane rejected by the dominant energy condition, at t 
again a small value of a (when C 7^ 0) rules out more of the parameter space. 



-b/2a, = n for 3 values of ^. Once 



6.5 Area Distance 



In § 5.2.3| we considered the area distance-redshift relation, with respect to the 'size' of the observed spatial sections. 
This is because the zero in the distance-redshift relation at small redshift, for certain model parameters (namely small 
b), is probably ruled out by present observations. Since the spacetime is conformal to a spherically symmetric spacetime, 
a zero in rA{z) would correspond to a star situated opposite to an observer (ie, at x = for all ip) would become spread 
all over the observers sky and have an infinite apparent magnitude - it would look as bright as if the star were right next 
to the observer. Similarly, an object close to a pole would become distorted in a similar manner to normal gravitational 
lensing - see MacCallum and Ellis (1970) for the distortion formula. As in §5.2.3, we make the assumption that such 

an effect would have been detected|^. As before we consider the cases 2:^ > 2 and z^^ > 5. 

To proceed, then, we simply need look at the z{x) function at x = tt and at a coordinate age given by ( 6.12| ); ie, we 
must find the range of parameter space for which 



z{x = tt) > 2, or 5. 



(6.14) 



We are dealing with a 3-dimensional parameter space so we must proceed with care. In figure |6.7] we show a 3-dimensional 
plot of z for a variation of the parameters a and C,, for b — 0.7. We see from this that for large C the redshift at which 
the first zero in the rA{z) function occurs at a very low redshift. The second thing to notice is that, even though a 
varies between —40 and the curves don't change much. In fact, if we look at figure we see that a doesn't do much 
damage at all: b and C seem to be most important. In chapter ^ this was also the case; most of the problems came 
from small b (cf, fig p^ ). In figure 6.£, we demonstrate that this is the case: small b is ruled out, while large C is only 
permitted if b is s uffi ciently close to 1. 

Also in figure |6.9| , we show that for the first zero to occur 'behind the CMB' requires b ^ 1. 

^Onc might argue that if such a phenomenon were to be detected it would require a systematic comparison of close-by galaxies with ones 
much further away, and at much earlier stages of their evolution. This would involve deciding what a particular galaxy would look like at 
say 1/4 or 1/2 its age. 
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Figure 6.3: Age exclusion plot of ^ vs 6 for a = —5 and Hq = 60 km s ^ Mpc ^ at i/) = tt. The contours mark the limits of the 
allowed ages. 



In order to u nderstand the effects that ( has on the spacetime, we need to study the a — b plane for various values of 
C. In figure 6.10 we show this exclusion plot, from which we can immediately see that a large C has the effect of moving 
the excluded region towards larger b, while at the same time making a more and more redundant as a parameter for 
this test. Physically the effect of a large C is to decrease the 'size' of the spatial sections, as one looks into the past; the 
first zero in the r^(z) function will occur closer to the observer. 



6.6 The CMB Anisotropy 

We assume for simplicity that the CMB is a blackbody which occurs at T = 3000K, which is roughly the temperature 
of decoupling. As it is a blackbody, the temperature evolves according to (see MacCallum and Ellis 1970) 



Tl 



'today I + Z 

Our observer at some location given by ip will measure the temperature of the CMB to be given by 



1 + XcMB,-*?) 



Now from (4.29) and (4.24) we may write 

1 + z{lp,XcMB,'d) = 1 + Zo(V',Xcmb) + Zl{ll^, Xcmb) COS 



(6.15) 



(6.16) 



(6.17) 
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Figure 6.4: Age exclusion plot of vs 6 for various a with Hq — 60 km s ^ Mpc ^ at 1/' = ^- The contours mark the limits of 
the allowed ages. 
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where ^Xcmb^= 'X(t{xcMB)y ■ The multipole moments in the CMB temperature fluctuations will then become 
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The strategy is then as follows: 

1. Given a set of parameters, a, 6, C, ^, ^/^ we solve ( |6.12 ) for the coordinate age T. 



(6.18) 



(6.19) 



(6.20) 



(6.21) 



2. Given our complete parameter set a, 6, C, ^, V', 2^ we then find the mean distance to the CMB surface Xcmb by 
solving 

zoUcmb) = 1000 (6.22) 
for Xcmb- This must be done numerically because the equations are such a mess. 
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Figure 6.5: Age exclusion plot of ^ vs a for a — ~5 and Ho = 60 km s ^ Mpc ^ at -i/j = tt. The contours mark the limits of the 
allowed ages. As in the case ( — a, large negative a is not allowed by this test. 



3. We can then calculate the temperature anisotropy by c alcul ating zi using the value of Xcme 
dipole moment of the CMB temperature is then (using ( 3.17 )) 



Si = 



Zl 



1 



zo 



found above. The 



(6.23) 



We do not wish that the IRF models account for the entire dipole moment of the CMB, because there is good evidence 
that this is due to our peculiar velocity of the local group. However, there is no reason to say that a small part of the 
dipole cannot be of a cosmological origin: a variation of 10% of the dipole could in principle be of a cosmological nature. 
We take this constraint here; that is 61 ~ 10^"*, or zi ^ 0.1. 



So, that's the plan. We are principally interested in the effects a non-zero will have. In figure 6.11 we show a 
surface plot of \zi\ against b and ^. The contour zi =0.1 is shown, with the surface below this line representing the 
acceptable parameter range. This is for an observer a,t ip — tt/A. We wish to constrain the parameters from the most 
restrictive position; it is not obvious, when C 7^ where the 'worst' position is. Figure 6.12 demonstrates that the most 
restrictive position for this test is an observer at ^ = 7r/2, and we can use that in what follows. 

The other thing to notice in figure [6.12 is that ( behaves very much as a does in §5.2.4 (cf. Fig. 5.11), that is the 
excluded regions are 'fingers', with the position depending on 6, and the size depending on ^. 



In figure 3.13, we show the exclusion plot in the a — b plane, as in figure 5.11, but for different values of C,. As can 
be seen, a high value of C is not exactly 'ruled out,' but makes things more 'unlikely' to fit. For any value of a, the 
range of allowed b decreases its total size, but never (except for very high C) rules out the models completely. This is 
in contrast to the other constraints, such as the size constraint, which completely rules out small values of b; or the 
age constraint, which rules out large negative values of a. Note that all the plots don't quite get close to 6 = 1, as the 
calculation becomes very (computationally) difficult there; it should be clear from the previous chapter that this region 
is ruled out. 
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Figure 6.6: Age exclusion plot of 6 vs a for various values of and Ho — 60 km s ^ Mpc ^ at t/j = tt. These curves are all for 
an age of 10 Gyr. Increasing restricts high values of b. 



6.7 The Magnitude- Redshift Relation 

Due to the complexity of the equations - basically the conformal factor - a proper analysis of the local inhomogeneity 
is really beyond the scope of this thesis, although not impossible in principle. Instead we will consider some examples 
to get an idea of the effect a large value of ^ will have. It is perhaps appropriate to use the magnitude-redshift relation 
for this purpose, because it is easier to compare with observations. 
To recap, we have 



m(z) — M — 25 = 5 log^g = 5 log 



Rq 



10 



51ogio{(l + z)|sinx(^:)|}, 



(6.24) 



where xi^) is given by the solution of 

i?o 



1 + z{x) 



WoRix) 
sin X cos I? 



^+{x) + ix) COS ?A cos X ]=- cos ^ sin -tp cos x 



^-ixf sin^ V + 4 ^^^y V - cos2 ^ sin^ ^) 



c2A 



-2 ;= cos 4 sm 2yj 



(6.25) 



with R{x) = X{t{x)), etc. 

We note first of all in figure [6.14 , we can see that a large C can cause all sorts of problems, with the situation 
illustrated in figure |6.14 being completely ridiculous. It is not all the fault of C though, as figure 3.15, demonstrates. 
A much higher value of b is used, and the problem is slightly relieved. However, the situation is still unacceptable as 
a cosmological model. It is clear that a blueshift in one direction in the sky is not reasonable - the results are similar 
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Figure 6.7: Exclusion plot of vs. a for b = 0.7 and Ho = 60 km s ^ Mpc ^ for the first zero of the rA{z) function. 



whatever the observers location in the spacetime, ip. However the results do depend on ^; with ^ = 7r/2 (as all the 
graphs here have) we get a large dipolc moment, which virtually disappears when ^ = 0. This is because when C is large, 
the dominant part of the dipolc in 1 + z comes from the s econd term in the square root. 

A far more reasonable situation is shown in figure 6.16, with a much smaller C. The observers location is at ^ = 7r/2, 
with the direction to the center given by ^ = 7r/2, which seems to be the 'worst case' scenario. 

The moral seems to be to have C and \a\ quite small, and at the same time b should be 'close to' 1. It is interesting 
though, that when C is distinctly non-zero, the dipole in the magnitude-redshift relation is more or less the same 
regardless of the observers spatial location, but does depend heavily on their direction to the center. 



6.8 Summary 

It this chapter we have considered the effects of a non-zero (■ We have shown that the age may be high enough, the first 
zero of rA{z) does not occur too close to the observer, and that the CMB anisotropies are not too large, for non-zero 
values of C- However, we have shown that a large value of C will introduce a significant dipole into the magnitude-redshift 
relation. This is probably enough to rule out these models, unless ^ is close to 0: this would violate the Copernican 
principle. This deserves further consideration. 
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Figure 6.8: Exclusion plot of ( vs. a for 6 = 0.7,0.5 and Ho = 60 km s ^ Mpc ^ for the first zero of the rA(z) function. This 
clearly shows that only b affects this function significantly. 
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Figure 6.9: Exclusion plot of ( vs. 6 for a = 



— 10 and Ho ~ 60 km s ^ Mpc ^ for the first zero of the rA{z) function. 
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Figure 6.10: Exclusion plot of a vs. b for various values of for the first zero of the rA{z) function. 
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Figure 6.12: Exclusion plot of ( ys b for two values of tp for the CMB anisotropy. a — —5. 
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Figure 6.13: Exclusion plot of a vs 6 for three values of for the CMB anisotropy. As increases, the 'fingers' move in to exclude 
more and more of the a — b plane, leaving a few narrow strips and isolated pockets. 
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Figure 6.14: The magnitude-redshift relation in opposite directions in the sky for an observer at 1/) = 0. We have b — 0.5, a = —10, 
and C = 100. 
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Figure 6.15: The magnitude-redshift relation in opposite directions in the sky for an observer at 
and C = 100. 



= 0. We have b = 0.95, a = 



-10, 
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Figure 6.16: The magnitude-redshift relation in opposite directions in the sky for an observer aX iji = 7r/2. We have b = 0.9, a ■ 
— 5, and C = 5. 



Chapter 7 

Conclusions and Future Work 



This thesis has examined the apphcation of the Copernican principle in inhomogeneous spacetimes. The class of models 
we considered were derived in chapter ^ under the condition that they admit an isotropic radiation field, and were 
shown to be a subclass of the Stephani solutions with symmetry. A cosmological model which satisfies the Copernican 
principle must satisfy Theorem ^ if one wants a spacetime with observers who see an (exactly) isotropic CMB. The IRF 
models we derive are generalisations of the FLRW models, that satisfy Theorem ^ of Ehlers, Geren and Sachs (1966). 
It was also demonstrated that the models we find all admit a thermodynamic scheme. The rest of the thesis was spent 
developing the relevant observational relations, and 'testing' various aspects of them. 

In chapter ^, we derived all the necessary cosmological tests. After the conformal symmetries were recognised, 
and the relevant coordinate transformations were made, this became a tractable problem. Indeed this transformation 
enabled us to 'rotate' the origin of coordinates to an arbitrary location, which gave the opportunity to derive all the 
observational relations valid for all observers in the spacetime. This is clearly crucial for the subsequent observational 
consideration of the Copernican principle. 

It chapters ^, and ^, we discussed observational constraints. We limited the coordinate time of observation using the 
present constraints on Hq by constraining the expansion rate with (5.27) and ( 6.12| ). For any choice of parameters this 
gave a value of the coordinate T. We then proceeded to examine the models on a variety of tests; the most important 
being age and the anisotropy of the CMB. Obviously, any plausible cosmological model must have reached the required 
age by the time the expansion has slowed to the rate we measure today. This is position dependent; on a given surface 
of constant coordinate time (determined by the expansion rate) the observers at the center will generally be older than 
their counterparts away from the center (cf, figure |]^). As with all the tests, we consider only the 'worst case' observer 
position; ie, that which excludes most of the parameter space. 

The CMB is also a crucial test. Despite the fact that the IRF models admit an isotropic radiation field by definition, 
they do not exhibit an exactly isotropic CMB. This is because the models are homogeneous at the big bang t = 0, 
or = 1 (and big crunch, incidentally), which is when the isotropic radiation field in question is 'emitted'. The CMB is 
not emitted then, but a short time after when the universe was inhomogeneous, but not drastically so. The anisotropics 
in the observed temperature are oc 1/(1 + z) ~ 1/W: that is all anisotropics arise from the inhomogeneities in W at the 
time of emission, and have nothing to do with the subsequent evolution and motion of observers. This means that, since 
t is small, W ^ 1, and the inhomogeneities are very small too; the conformal flatness of the spacetime allows the light 
from the CMB surface to travel unmolested by any subsequent deviation from homogeneity, and we can see the CMB 
as a reflection of the universe's inhomogeneity as it was at that time. The dipole moment of the CMB is the largest 
moment, which we restrict to be smaller than 10% of the observed dipole, on the assumption that we are unable to say, 
at present, that the dipole is entirely due to the motion of the local group. 

The other tests we considered were that the energy conditions are satisfled (see below) and that the flrst zero of the 
m — z relation does not occur too close to the observer. Both these tests are independent of observer location^. As we 
have discussed, these do not impose strong constraints on the available parameter space. 

The IRF models we chose to look at had the scale factor (arbitrarily) restricted to be a quadratic function of time. 
We must conclude, therefore that the allowed degree of freedom in the IRF models, which manifests itself in the freedom 
in choosing the scale factor, R{t), gives us an enormous range of models which would satisfy the constraints chosen here 
to a perfectly acceptable level. In fact chapters |[ and || may be considered as an existence 'proof that the IRF models 
are acceptable models of the universe on observational grounds; that their matter content is unappealing does not enter 

^Tho dominant energy condition actually depends on observer position at certain times. 
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into it. 



7.1 A Note On 'Best Fits' and a > 

The discussion in chapters |^ and || is obviously not complete. We should really provide a fit of the m — z relation to 
real data. However, this would not be simple in general, because any result would be position dependent. Moreover, 
because of the dipole mentioned above, the data would have to include the directional coordinates of the SNIa. As the 
complete SNIa data is yet to be made public, we are not able to do this as yet. 

However, in the spherically symmetric case, an observer at the center would see no dipole variation of the m — z 
relation, so it is of value to fit the models from this special location, if only to show a proper fit (ie, as good a fit as may 
be obtained with an FLRW model) is possible. In |5.4| we demonstrated that this was possible. In fact it is possible to 
show that any FLRW m — z relation may be 'matched' by a suitable choice of parameters in these IRF models (Barrett 
and Clarkson, in preparation), for an observer at the center. 

It should be clear from the discussion in the introduction on the active gravitational mass, that a fit to an FLRW 
model with go < will violate the strong energy co nditio np]. This means that local observations at the center will force 
a > (cf, § [3.1.l| ). Moreover, we have seen from § 5.2.l| , that requiring the normalised energy density to be less than 



unity today also requires a > 0. This implies that this thesis is incomplete. As we have noted before, a > implies that 
the universe will open up at <= (1 — b)/2a (in the spherically symmetric case), after which the universe will become 
infinite in spatial extent. Spatial infinity, in terms of this will happen at = (cf, |5.8| ); 

sin't^ = ^<l for t>'—^. (7.1) 

Applying the Copernican principle in this case becomes difficult - at least in terms of producing exclusion diagrams. It 
may not be possible to find a 'worst case' observer position, as we have done up to now. If we choose any -ip < V'maxi 
then the number of observers between the center and ■0 will be finite; whereas, between ip and V'max the number will be 
infinite. Thus it is impossible to say that for any we have considered all or most observers. 

However, it may be possible to take a constraint, say the CMB dipole, and consider the limit ^ — s- V'max- That is we 
may say that if 

lim |zi|<0.1 (7.2) 

then the model will pass the test, and satisfy the Copernican principle. 

This however, is a computational nightmare - as we take the limit, at any given ip, we must solve (6.18) for Xcmb- 
This must also be done over the whole 3-dimensional parameter space. This is a chore for the future. 



7.2 The Local Dipole 

The most interesting question this thesis has raised is the possibility that we may be living in a distinctly inhomogeneous 
universe (that is, not a slightly inhomogeneous perturbed FLRW universe). 

Acceleration makes the IRF models inhomogeneous, but will leave its mark in the local universe. This is testable 
(j j3.3| ) by considering the local dipole moment of Hq - ie, the local dipole in the linear Hubble law. The measured dipole 
is usually interpreted as our peculiar velocity with respect to the CMB frame; that the directions of our peculiar motion 
from local studies and that of the CMB frame line up to within about 20° (Saunders et at, 1999; Scaramella et ai, 
1991; Riess, 1999) lends support to this but it cannot be conclusive. It should be noted that in the IRF models, we also 
require the local and CMB dipoles to line up. 

If we wish to prove homogeneity then we must test this local dipole to determine if it grows linearly with distance. If it 
does then we may use the slope of this linear relationship to either (cf, 5.30| ): 1. determine the maximum inhomogeneity 



of the universe, assuming we are located at or near the 'equator' {ip — 7r/2); or, 2. determine our distance from the 
center, in an arbitrarily inhomogeneous model. Both would require violation of the Copernican principle to some degree. 
The first would require us to be situated at the point of maximal inhomogeneity (maximum acceleration) ~ which is a 
'non-typical' location; while the second may violate it depending on how inhomogeneous we may decide the universe to 

^This may not be quite true in general (the results are only preliminary); from the discussion in § ^.e| , it should be clear that things are 
not that simple when comparing observationally derived values of go to the 'length scale' derived version. It is true at the center, however, 
where the acceleration is zero and the negative pressure generates a negative go , regardless of definition. 
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be. That is, the more inhomogeneous the model (the greater the acceleration at the equator) the closer we must be to 
the center, and we would no longer be representative of 'most' observers. 

Preliminary investigations reveal that constraints on the dipole in Hq are weak - around 20% may be attributed to 
an acceleration term (Clarkson, Rauzy, and Barrett, in preparation). This has been investigated using the IRAS galaxy 
catalogue; however, the recently released PSCz catalogue with around 15000 galaxies may provide stronger constraints. 

7.3 Almost EGS Considerations 

One may expect that the IRF models derived in § |3.1| would extend to an 'almost' theorem, in the same manner as the 
almost EGS theorem of Stoeger, Maartens and Ellis (1995). If we take an almost isotropic radiation field characterised 

by 

aab^Oie), (7.3) 
V[, {ii,] - leu,]) = 0{e), (7.4) 



where epsilon is small in some sense then we may write, from equations ( |3.3| ) 

Audita + ^af^n ^ 0{e), (7.5) 

(7.6) 

or ^ 

Wa = -iV,ln^^ + 0(e) (7.7) 

which gives 

V[a?ib] - V[aVb]//fl - iUbaO ^0{e): (7.8) 

there doesn't seem to be any way to make the rotation small, unless one artificially introduces additional assumptions. 
Therefore perfect fluid solutions which admit an almost isotropic radiation field may not (at first sight at least) be 
perturbations of the IRF models. However, there may be ways to constrain rotation on the basis of local observations: 
rotation leads to anisotropic number-counts - see Fennelly (1976). It would seem unlikely that such observations would 
be anywhere near as accurate as those of the CMB anisotropics; ie, we would be unable to claim that ujab = 0(e) from 
number-counts alone. 



7.4 Further Considerations 

The IRF models, by definition, admit an isotropic radiation field. As all observers see such a radiation field to be 
isotropic, the isotropy of the CMB cannot be used to verify the cosmological principle. The usual suggestions of testing 
the cosmological principle directly, necessarily involve making measurements within our past lightcone. The Sunyaev- 
Zel'dovich (SZ) effect is often quoted as a suitable physical effect to produce such a test (Goodman, 1995). Such a test 
will involve measuring the anisotropy of the CMB around other observers. It is normally assumed that if the anisotropy 
is small around other observers then one may conclude that we are not at a special location (eg, at a center of symmetry) , 
and that the universe is FLRW - hence the cosmological principle is verified. 

This, of course, assumes that the acceleration is zero. The IRF models would also give the same effect; measurement 
of the CMB anisotropics as seen by another observer would be similarly small, but these would not imply that the 
universe is homogeneous, by Theorem |^. Hence, contrary to popular belief, the SZ effect cannot possibly show that the 
cosmological principle is a valid assumption. 

Clearly, the SZ effect cannot distinguish between different IRF models - thus the acceleration may be arbitrarily 
large. This means that the Copernican principle may be violated, because for our local dipole to be so small (at most 
20% of Hq) would necessarily require us to be close to the center of such an inhomogeneous universe. Therefore the SZ 
effect will not even be able to provide any proof of the Copernican principle. 

It should be noted that acceleration would leave a small signature in such CMB measurements: small dipoles around 
other observers would be correlated, rather than a random distribution of dipoles from their peculiar motions. The 
distribution of the dipoles would be a mixture of some pointing in one direction, and others pointing in the opposite 



CHAPTER 7. CONCLUSIONS AND FUTURE WORK 



101 



direction. In the spherically symmetric case, the direction of the dipole would depend on the number of times the 
scattered light had passed through the point opposite to us. This will depend on the redshift of the galaxy. It is not 
clear what pattern the more general IRF models would give. 

If we used the SZ effect and found such a pattern of dipoles, then this would imply that acceleration is present in 
the universe. Such accuracy for these experiments is a long way off, as it would be next to impossible to disentangle 
any dipole from peculiar velocities and that from acceleration we have enough difficulty in determining our own. 

This poses some interesting questions. It would be interesting to carry on this work along the lines of Stoeger, 
Maartens and Ellis (1995), and Maartens, Ellis and Stoeger (1995), which is to place limits on all the components which 
characterise a cosmological model: {u°',9,aab,0Jab,Eab,^Hab}, together with the matter content. In general this is not 
a simple task and, as this thesis makes clear, one cannot rely on the CMB alone. 



Appendix A 



Coordinate transformations and the 
Stephani Spacetimes. 



On the face of it the Stephani models admitting an isotropic radiation field in ( |3.19D depend on one free function and 
ten free parameters. However, it is possible to use coordinate transformations on the spacetime to eliminate many of 
these parameters, resulting in a considerable simplification. As is shown in Barnes (1998), conformal transformations of 
the coordinates on the hypersurfaces of constant time preserve the form of the metric but change the free functions a, 
b, and c. These transformations can be thought of as acting on the five-dimensional space spanned by a, b and c and 
constitute the Lorentz group in five dimensions, 5*0(4, 1): they leave —ab+ |cp invariant (a and b are 'null coordinates'). 
It will be convenient here to let a = a + /3 and b = a — (3 (so that the transformations preserve — + + |cp), and to 
adopt five- vector notation: q^^ = {a,(3,c). We will use the terms 'rotation' and 'boost' to refer to the transformations 
on q, and will call q timelike, spacelike or null if — + f3'^ + |cp is negative, positive or zero, as usual. Then it is easy 
to visualise the transformations on the free functions by imagining the 'mass hyperboliods' of representations of the 
Lorentz group in the usual way: a timelike vector can always be boosted so that it has the form (a, 0,0), whereas a 
spacelike vector can be boosted and rotated into (0, /3, 0), for example. 

In addition to the Lorentz transformations we also have the freedom to change basis in the function space spanned 



by the free functions. For the spacetimes of interest here, described by (3.19), it is desireable to preserve f2{t) = 1, so 
that the basis change is 

+ (A.l) 



In five-vector notation the equations (3.19) become 

q^^{t) = {a,l3,c)^q'iT{t)+q^ (A.2) 
(with qi and q2 constant vectors), and the goal is to reduce as many of the components of q^ and ^2 to zero as possible 



using Lorentz transformations in the five-dimensional space containing qi and 92, and the basis change (A.l). Note that 
the FLRW {d = 1) subcase of ( |3.19 ) is characterised by the hnear dependence of qi and 92- 



It is easy to see that c (which breaks the spherically symmetry of the metric ( jS.q )) may always be reduced to the 
form c = cz (z = (0,0,1)), with c a constant: perform a spatial 4-rotation to reduce to ~ (q!i,/3i,0), then a 
spatial rotation amongst the c-components (which obviously leaves qi unaffected) to give ql^ — {0.2, P2, cz). 

It is possible in general to make further simplifications, but precisely how qi and 52 are simplified depends on 
whether they are spacelike, timelike or null. For example, if either qi or (72 is timelike (or may be made timelike by a 
transformation ( |A.l[ )) it is possible to reduce the model to manifestly spherically symmetric form (c = 0): boost so that 
the timelike vector, say qi, becomes qi — (ai, 0, 0) and rotate spatially so that the c-components of the other vector are 



also zero, 52 = (q;2,/?2,0) (we could then use (A.l) to eliminate more of these constants). 



To summarise, we have demonstrated that it is always possible to reduce the Stephani models of (3.19) to the form 

a{t) = aiT{t) + a2, 

b{t) = biT{t) + b2, (A.3) 
c(t) — cz, 

(where we have transformed back from a and (3 to a and b), and when either of the qi or 52 is (or may be made) timelike 
we can set c — 0. 
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Finally, note that we may always assume that ai ^ in (A.2), because if ai = then bi ^ 0, otherwise Vt = 0, 
and it is possible to perform a coordinate inversion x i— > x/r^ that interchanges a and &). This does not exhaust the 
possibilities for simplification: we could, for example, use (A.l) to set a-i — 0. 



Appendix B 

Transformation to a Non-Central Position. 

We want to transform from the (x, ^, ^) coordinate system, whose origin is at the centre, to coordinates centred instead 
on some observer at % = while preserving the form of the FLRW part of the metric (O). It is therefore necessary 
to identify the transformations of the (homogeneous) FLRW spatial sections that leave the FLRW metric invariant, 
ie, the isometries of the spatial sections. This is simple. Since the spatial sections of an FLRW model with positive 
curvature constant (A > 0) are 3-spheres, the isometries we require are 4-dimensional rotations (ie, elements of SO{A), 
the isometry group of the 3-sphere). 

A sphere of radius R in 4-dimensional space with cartesian coordinates {x, y, z, u) is defined by 

a;^ + + + u'^ = 

We have three coordinates on this sphere: x ^^'^ the two spherical polar angles 9 and (j). These are related to the 
cartesian coordinates by 

X = i? sin X sin cos (B.l) 

y = i? sin X sin sin (/) (B.2) 

z = i?sinxcos6' (B.3) 

u = Rcosx (B.4) 

The origin, x = Oi is then a,t x = y — z — 0, u — R. We are only interested in rotations that move the origin, and, 
as the initial metric is spherically symmetric (really spherically symmetric, not just conformally: even the conformal 
factor is spherically symmetric about the centre), we need only consider moving the observer in one direction, which 
we choose to be the 'z direction' (ie, to a position with non-zero z, but x = y = 0). Clearly, then, we are looking for a 
rotation in the u — z plane. Since we have the conformal factor as a function of x we want to find x as a function of 
the new coordinates. Starting with coordinates x': ^' and </>', centred on some position x — "01 along with their primed 



cartesian counterparts a:', y', z' and u' (which are related in the same way as the unprimed coordinates in (B.1)-(B.4)), 
a rotation back to the original coordinates is given, in cartesian coordinates, hy x — x' , y ~ y' and 

z — cos-i/jz' sin-f/'u', 

u = — s\n^z' -\- cosipu . (B.5) 
(Note that at the origin of the primed coordinates, where z' = and u' — R, we have u — Rcosip, showing that x = "0 



there, as required.) Equation (B.5), along with the primed versions of (P3.3D and (B.4), then immediately gives 



cosx = cos cos x' — sin sin x' cos 6*', (B.6) 



and this is all we will need, since the only spatial coordinate that enters into the original metric (4.9) is X: and that 
enters only as cos x (2 sin^ ^ = 1 — cos x) ■ 
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